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ABSTRACT 

Open  form  expressions  for  the  velocity  distribution 
and  volumetric  flow  rates  were  derived  for  the  co- current 
stratified  fully  developed  laminar  flow  of  two  immiscible 
liquids  in  a  rectangular  conduit.  The  effect  of  the  presence 
of  the  less  viscous  layer  was  to  reduce  the  pressure  gradient 
necessary  for  the  transport  of  the  more  viscous  phase  at  a 
given  rate  by  a  factor  which  was  a  function  of  the  input 
ratio,  the  viscosity  ratio  and  the  conduit  aspect  ratio. 
Pressure  gradient  reduction  factors  and  holdup  ratios  were 
evaluated  with  the  aid  of  an  IBM  1620  computer  for  viscosity 
ratios  from  1  to  1000  and  aspect  ratios  from  7.896:1  to  1:3. 

The  pressure  gradient  reduction  factor  increased  concomitantly 
with  viscosity  ratio  and  aspect  ratio.  For  a  viscosity  ratio 
of  1000  the  maximum  pressure  gradient  reduction  factor  was 
3.38  for  an  aspect  ratio  of  7.896:1  and  only  1.13  for  an 
aspect  ratio  of  1:3. 

Experiments  were  performed  with  oil  and  water  having 
a  viscosity  ratio  of  5.326  and  a  rectangular  conduit  having 
an  aspect  ratio  of  7.95:1  and  an  overall  length  of  37.5  feet. 
For  given  input  conditions  measurements  were  obtained  of  the 
insitu  ratio,  the  pressure  gradient  and  the  interface  velocity. 


» 
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In  the  laminar- laminar  regime  no  waves  were  present  at  the 
interface  and  reasonably  good  agreement  with  the  theory  w as 
obtained.  Coincident  with  the  transition  to  turbulence 
in  the  less  viscous  water  phase  the  interface  became  disturbed 
by  small  two-dimensional  waves,  the  velocities  of  which  were 
essentially  the  same  as  the  actual  average  water  velocities. 

These  waves  appeared  to  have  their  origin  in  the  pressure 
fluctuations  in  the  water  phase.  At  high  water- to-oil  input 
ratios  large  two-dimensional  and  roll  waves  were  observed 
which  were  relatively  slow  moving  and  were  apparently  caused 
by  an  instability  of  the  mean  flow.  Wavelengths  and  amplitudes 
of  these  waves  were  estimated  using  a  stereophotogrammetric 
technique.  In  the  laminar- turbulent  regime  pronounced  devia¬ 
tions  from  the  laminar- laminar  theory  were  evident;  in  particular, 
the  pressure  gradient  reduction  factors  were  always  much 
less  than  the  corresponding  laminar- laminar  values.  Lockhart- 
Martinelli  type  parameters  were  found  useful  in  correlating 
the  pressure  gradient  data  in  the  laminar- turbulent  regime. 

The  transition  to  turbulence  in  the  oil  phase  caused  the  two- 
dimensional  characteristics  of  the  interfacial  waves  to  be 
destroyed  and  at  high  flow  rates  of  both  phases  three-dimensional 


waves  we re  observed. 
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I  INTRODUCTION 


One  of  the  more  interesting  and  potentially  useful 
phenomena  associated  with  the  simultaneous,  horizontal  pipeline 
flow  of  two  incompressible  fluids  is  the  fact  that  the  pressure 
gradient  and  power  requirement  necessary  for  the  flow  of  the 
more  viscous  phase  at  a  given  rate  may  be  substantially  reduced 
by  the  presence  of  the  less  viscous  phase.  The  most  common 
practical  example  in  which  this  effect  is  evident  is  the  pipe¬ 
line  flow  of  petroleum  and  water. 

The  advantage  to  be  derived  by  adding  water  to  petroleum 
flowing  in  a  pipeline  was  first  suggested  by  Isaacs  and  Speed 
in  a  patent^7)  obtained  in  1904.  However,  it  is  only  in  recent 
years  that  theoretical  and  experimental  studies  have  been  under¬ 
taken  with  the  objectives  of  obtaining  fundamental  information 
concerning  the  mechanics  of  this  two-phase  system  and  quantitative 
predictions  of  the  extent  to  which  the  pressure  gradient  may 
be  reduced. 

One  important  result  of  these  studies  which  will  be 
subsequently  discussed  in  detail  is  that  concentric  flow,  in 
which  the  less  viscous  liquid  completely  surrounds  the  more 
viscous  liquid,  is  much  more  effective  in  reducing  the  pressure 
gradient  than  stratified  flow.  Concentric  flow  can  be  established 


2 


if  the  liquids  are  introduced  concentrically  into  the  pipeline 
and  the  density  differential  between  the  liquids  is  almost 
non-existent.  However,  in  the  presence  of  even  a  slight 
density  differential  the  liquids  tend  to  stratify  and,  unless 
the  pipe  is  strongly  wet  by  the  less  viscous  phase,  stratified 
flow  results  with  part  of  the  pipe  surface  wet  by  the  more 
viscous  phase.  The  fact  that,  in  general,  pairs  of  immiscible 
liquids  have  different  densities  makes  stratified  flow  the 
much  more  widespread  and,  therefore,  important  configuration. 

Studies  of  the  stratified  flow  of  two  incompressible 
fluids  between  infinitely  wide  parallel  plates  and  in  circular 
pipes  have  contributed  to  the  knowledge  of  this  flow  system 
but  have,  at  the  same  time,  exposed  the  extremely  complex 
nature  of  the  interaction  between  the  phases,  especially  when 
waves  are  present  at  the  interface. 

This  introduction  consists  basically  of  a  critical 
review  of  the  literature  which  is  relevant  to  the  co-current 
horizontal,  stratified  flow  of  two  immiscible,  incompressible 
fluids  with  a  view  to  establishing  the  objectives  of  the 
present  investigation.  Certain  concepts  and  characteristics 
of  single  phase  flow  will  be  discussed  before  considering  two- 


phase  systems. 
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A.  THE  CONCEPTS  OF  LAMINAR,  TURBULENT  AND  TRANSITIONAL  FLOW 

The  concepts  of  laminar,  transitional  and  turbulent  fluid 
flow  are  basic  to  any  theoretical  or  experimental  study  in 
fluid  mechanics.  Consider,  as  an  example,  the  flow  of  a  fluid 
through  a  straight  circular  pipe  with  smooth  walls  and  constant 
cross  section.  At  low  flow  rates  every  fluid  particle  moves 
with  uniform  velocity  parallel  to  the  pipe  boundary.  Viscous 
forces  slow  down  the  particles  near  the  wall  relative  to  those 
in  the  central  core.  This  well-ordered  type  of  flow  is  termed 
laminar  or  viscous  flow.  However,  at  higher  flow  rates  laminar 
flow  ceases  to  exist  and  a  mixing  of  the  particles  occurs,  for 
superimposed  on  the  axial  motion  are  subsidiary  motions  normal 
to  it.  Such  fluid  motion  is  described  as  turbulent  and  has 
been  defined  by  Hinze^^^  as  "...an  irregular  condition  of  flow 
in  which  the  various  quantities  show  a  random  variation  with 
time  and  space  co-ordinates,  so  that  statistically  distinct 
average  values  can  be  discerned". 

The  two  types  of  flow  were  first  observed  by  Reynolds 
who  injected  a  thin  filament  of  dye  into  the  flowing  stream. 
Under  laminar  conditions  the  dye  filament  maintained  its 
identity,  but  under  turbulent  conditions  eddies  caused  the 
filament  to  diffuse  throughout  the  fluid.  Reynolds  found 
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that  the  quantity  PV--  ,  now  known  as  the  Reynolds  number, 
in  which  p  is  the  density  of  the  fluid,  V  is  the  average  flow 
velocity,  D  is  the  pipe  diameter  and  JJ,  is  the  viscosity  of  the 
fluid,  is  important  in  determining  the  conditions  under  which 
the  laminar  motion  is  unstable  and  becomes  turbulent.  The 
Reynolds  number  may  be  shown  to  represent  the  ratio  of 

inertia  forces  to  viscous  forces.  Inertia  forces  tend  to 
promote  instability  while  viscous  forces  tend  to  damp  out 
disturbances  and  maintain  laminar  flow.  The  critical  value 
of  the  Reynolds  number,  i.e.  the  value  at  which  the  transition 
from  laminar  to  turbulent  flow  occurs,  was  determined  by 
Reynolds  to  be  about  2300.  It  is  now  established,  as  was 
suspected  by  Reynolds,  that  the  numerical  value  of  the  critical 
Reynolds  number  depends  upon  the  size  of  disturbances  in 

(27) 

the  entrance  to  the  pipe.  This  has  been  demonstrated  by  Ekman 
who  maintained  laminar  flow  up  to  a  Reynolds  number  of  40,000 
by  the  use  of  an  inlet  section  exceptionally  free  from  distur¬ 
bances.  More  recently  Leite ,  who  also  used  a  specially 
designed  entrance  section  to  stabilize  the  flow,  found  laminar 
flow  to  be  stable  up  to  a  Reynolds  number  of  20,000.  In 
addition  he  artificially  created  disturbances  of  known  ampli¬ 
tude  and  found  that  the  transition  to  turbulent  flow  occurred 
when  the  amplitude  of  the  disturbance  exceeded  a  threshold 


. 
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value.  The  value  of  the  threshold  amplitude  decreased  with 
an  increase  of  Reynolds  number. 

While  there  is  doubt  about  the  upper  value  of  the  critical 
Reynolds  number,  little  doubt  remains  that  below  a  Reynolds 
number  of  approximately  2000  laminar  flow  is  stable  despite 
the  presence  of  large  disturbances. 

The  transition  to  turbulent  flow  causes  profound  changes 
in  the  velocity  distribution  and  the  law  of  resistance.  Whereas 
in  laminar  pipe  flow  the  velocity  distribution  is  parabolic, 
in  turbulent  flow  the  higher  lateral  transport  of  momentum 
causes  the  velocity  distribution  to  be  much  more  uniform.  Under 
laminar  flow  conditions  the  pressure  gradient  required  to  main¬ 
tain  the  flow  is  directly  proportional  to  the  average  flow 
velocity,  while  in  the  turbulent  regime  the  pressure  gradient 
is  proportional  the  average  velocity  raised  to  about  the  1.75 
power,  reflecting  the  additional  energy  losses  associated  with 
the  eddies. 

The  transition  from  laminar  to  turbulent  flow  in  pipes 
has  been  described  by  Schlichting  who  used  the  recent 

experimental  work  of  Rotta^7^)  as  a  basis  for  his  discussion. 

In  the  region  of  the  critical  Reynolds  number  the  flow  becomes 
intemittent  in  that  it  alternates  in  time  between  the  laminar 


and  turbulent  states.  In  particular,  the  velocity  distribution 
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alternates  between  a  laminar  distribution  and  a  corresponding 
turbulent  distribution.  The  nature  of  the  flow  can  be  described 
by  an  intermittency  factor,  ,  which  is  defined  as  that  fraction 
of  time  during  which  the  flow  at  a  given  position  remains 
turbulent.  For  continuous  turbulent  flow  If  =  1  and  for  con¬ 
tinuous  laminar  flow  f -  0.  Rotta  found  that  the  transition 
was  completed  in  the  range  of  Reynolds  number  from  2300  to 
2600,  f  showing  a  steady  increase  with  Reynolds  number  through 
this  range.  The  intermittency  factor  was  a  function  of  the 
axial  distance  from  the  entry  section  and  increased  with  the 
distance,  indicating  that  the  length  of  the  turbulent  patches 
increased  as  the  fluid  flowed  along  the  pipe. 

Much  effort,  especially  in  recent  years,  has  been 
devoted  to  the  theoretical  prediction  of  the  critical  Reynolds 
number  for  various  flow  geometries.  The  basic  principles  of 
stability  theory  have  been  reviewed  by  Lin  and  Schlichting  . 

Briefly,  the  theory  involves  the  superposition  of  a  small 
disturbance  on  the  basic  flow  under  investigation.  Solution 
of  the  differential  equation  of  the  disturbance,  an  Orr- Sommer f eld 
type  of  equation,  with  the  appropriate  boundary  conditions 
indicates  whether  the  initial  disturbance  is  amplified,  in 
which  case  the  flow  is  unstable,  or  whether  the  disturbance  is 
damped  out  and  the  flow  is  stable.  The  mathematics  is  always 


7 


exceedingly  complex. 

The  theory  has  been  successful  in  predicting  boundary 
layer  instabilities  for  flow  over  a  flat  plate  which  agree 
well  with  experimental  observations  but  as  far  as  fully 

developed  flows  are  concerned  contradictions  exist.  Several 
theoretical  studies,  notably  those  of  Wasow^^),  Zondek  and 
Thomas (^7) ^  Gallagher  and  Mercer  (^4)  and  Deardorf  f  (^1)  ,  all 
based  on  small  disturbance  principles,  appear  to  indicate  that 
plane  Couette  flow,  where  the  velocity  distribution  is  linear, 
is  always  stable.  However,  turbulent  plane  Couette  flow  has 
been  observed  by  Robertson (75)  at  high  Reynolds  numbers. 

Lin  (60)  reports  that  flow  between  infinite  parallel  plates  or 
plane  Poiseulle  flow  has  been  shown  to  be  both  stable  and 
unstable  and  suggests  that  the  confusion  may  have  resulted 
from  the  exceedingly  complex  role  played  by  viscosity.  Viscous 
forces  can  serve  as  a  source  of  instability  by  creating 
stresses  which  transfer  energy  from  the  basic  motion  into  the 
disturbance.  No  evidence  of  the  instability  of  fully  developed 
flow  in  a  circular  pipe  has  been  found  theoretically (^1)  while, 
of  course,  turbulent  flow  is  observed  experimentally.  Lin  has 
suggested  that  the  origin  of  turbulence  for  flow  in  a  circular 
pipe  may  be  in  the  entry  section  where  a  boundary  layer  exists 
close  to  the  pipe  surface  and  the  flow  is  essentially  uniform 
in  the  region  of  the  pipe  axis. 
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It  would  appear,  therefore,  that  fully  developed  laminar 
flows,  except  perhaps  plane  Poiseuille  flow,  are  stable  as 
far  as  very  small  disturbances  are  concerned.  Observed  tur¬ 
bulence  would  appear  to  originate  either  from  boundary  layer 
instabilities  in  the  entrance  section  or  from  disturbances 
in  the  main  flow  which  have  a  sufficiently  large  amplitude, 
as  demonstrated  by  Leite(^). 

Recently  Benjamin^)  investigated  theoretically  the 
effect  of  a  flexible  boundary  on  the  stability  of  boundary 
layer  flow  and  found  that,  under  certain  conditions,  the  flexible 
nature  of  the  boundary  could  tend  to  stabilize  the  flow. 

B.  THE  GENERATION  OF  WATER  WAVES  BY  WIND 

The  generation  of  waves  by  wind  blowing  over  a  water 
surface  has  intrigued  philosophers  for  many  years.  As  Ursell 
has  pointed  out  the  problem  is  merely  one  of  mathematics 
provided  it  is  accepted  that  the  fluid  motions  are  subject  to 
the  laws  of  classical  mechanics  with  the  appropriate  boundary 
conditions.  Unfortunately  the  resources  of  mathematics  are 
insufficient  for  a  complete  solution,  although  progress  can  be 
made  when  suitable  physical  assumptions  are  made.  The  real 
problem  is  to  determine  the  proper  assumptions,  and  this  can 
only  be  done  through  a  combination  of  theory  and  experiment. 


« 
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Scott  Russell ' s observations  of  the  initiation  of 

(55) 

waves  on  water  have  been  quoted  by  Lamb'  '  .  Scott  Russell 
observed  that  the  water  surface  was  perfectly  smooth  and  reflected 
the  image  of  surrounding  objects  up  to  a  wind  velocity  of  about 
0.5  miles/hr.  At  a  velocity  of  about  1  mile/hr.  the  water 
no  longer  gave  a  clear  reflection  due  to  the  presence  of 
minute  capillary  waves.  These  waves  ceased  almost  simultan¬ 
eously  with  a  lull  in  the  wind  and  were  incapable  of  travelling 
any  distance  in  the  absence  of  the  wind.  On  a  further  increase 
in  wind  velocity  to  about  2  miles/hr.  small  gravity  waves 
arose  uniformly  over  the  whole  water  surface,  while  the 
capillary  waves  remained  superimposed  upon  them. 

Similarly,  Benjamin ^  observed  the  flow  of  air  over 
thin  layers  of  water  and  highly  viscous  syrup.  In  the  case  of 
water,  initially  some  apparently  random  agitation  of  the  water 
surface  was  observed,  while  long- crested,  slow  moving  waves 
subsequently  appeared.  In  the  case  of  the  highly  viscous  liquid 
the  random  disturbances  were  not  visible  but  long- crested 
waves  appeared  which  were  very  similar  to  those  on  the  water 
surface. 

Benjamin  suggested,  therefore,  that  waves  can  be  initiated 
in  two  basic  ways.  The  random,  small  disturbances  initially 
observed  on  the  water  surface  are  presumably  caused  by  turbulent 
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pressure  fluctuations  in  the  air  stream,  while  the  regular 
long-crested  waves  would  appear  to  have  their  origin  in  the 
instability  of  the  mean  flows  of  air  and  water  adjacent  to 
the  interface.  Two  entirely  different  theoretical  approaches 
are  therefore  possible  -  the  first  concerned  with  the  generation 
of  waves  by  random,  turbulent  pressure  fluctuations  and  the 
second  based  on  small  disturbance  theory,  similar  in  principle 
to  that  already  outlined  for  the  prediction  of  instability  of 
single  phase  laminar  flow.  The  theory  which  would  apply  in  a 
given  naturally  occurring  or  experimental  system  would  depend 
on  the  nature  of  the  observed  waves. 

1.  Pressure  fluctuation  theories 

The  generation  of  waves  by  random  fluctuations  of  pressure 
in  the  turbulent  air  stream  has  been  investigated  by  Eckart^^ 
and  Phillips (73) #  Eckart  assumed  a  random  distribution  of 
normal  pressure  acted  over  a  circular  region  of  the  ocean  which 
was  assumed  to  be  an  inviscid  fluid  and  incorporated  the  ideas 
of  statistical  correlation  analysis.  No  account  was  taken 
of  the  differences  in  pressure  between  the  wind-ward  and  lee¬ 
ward  slopes  of  the  waves,  nor  the  tangential  wind  stress.  In 
these  respects  his  approach  is  unrealistic,  and,  in  fact,  the 
heights  of  predicted  waves  are  less  than  those  observed. 
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Phillips's  theory  is  concerned  with  the  resonance 
effect  of  pressure  fluctuations  in  the  turbulent  air  stream 
on  an  initially  calm  water  surface.  Agreement  of  the  theory 
with  observations  is  good  and  Phillips  suggested  that  the 
pressure  resonance  mechanism  may  well  be  the  principal  means 
by  which  energy  is  transferred  to  waves. 

2.  Instability  of  mean  flow  theories 

The  Kelv in-Helmholtz  theory (^6)  ±s  based  on  the  simplest 
possible  model  in  which  uniform  flow  is  assumed  in  each  phase, 
the  velocity  distribution  being  discontinuous  at  the  inter¬ 
face.  Viscosity  is  neglected  so  that  the  motion  arising  at 
the  interface  from  a  small  arbitrary  disturbance  can  be  des¬ 
cribed  by  a  velocity  potential.  The  motion  is  controlled  by 
gravity  and  surface  tension.  The  critical  wind  velocity,  i.e. 
the  velocity  for  which  the  disturbance  will  neither  grow  nor 
decay,  is  predicted  to  be  about  650  cm/sec.  (15  miles/hr.), 
whereas  in  nature  waves  are  observed  at  wind  velocities  very 
much  below  this  value. 

Lock  and  Wuest  (^4)  have  extended  the  Kelvin- Helmholtz 

model  to  take  viscosity  into  account  i.e.  boundary  layers  are 
formed  on  both  sides  of  the  interface.  Outside  the  boundary 
layers  the  air  moves  with  uniform  velocity  and  the  water  is 
at  rest.  Lock  considered  a  velocity  profile  calculated  from  a 
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particular  point  upstream,  while  Wuest  arbitrarily  chose 
approximations  to  actual  boundary  layer  profiles.  Ursell 
remarks  that,  although  the  numerical  calculations  based  on 
Lock's  theory  are  far  from  complete,  the  theory  seems  to 
indicate  that  a  wind  velocity  of  100  cm/sec.  (2.5  miles/hr)  is 
capable  of  raising  waves,  and  that  it  appears  that  Wuest 's 
calculations  may  contain  some  errors. 

Feldman (^2)  considered  the  stability  of  a  horizontal 
liquid  film  in  contact  with  a  semi-infinite  air  stream.  He 
assumed  parallel  laminar  flow  with  a  linear  velocity  distri¬ 
bution  in  both  phases,  a  discontinuity  in  slope  existing  at 
the  interface.  The  critical  velocity  by  Feldman's  analysis 
is  an  order  of  magnitude  greater  than  the  experimental  velocity 
observed  by  Knuth^4)  for  annular  flow  of  air  and  water  in  a 
circular  pipe.  Feldman  suggested  the  discrepancy  could  be 
due  to  the  relative  stability  of  the  plane  Couette  flow  assumed 
in  his  theory  compared  with  flows  having  curved  velocity 
profiles . 

The  foregoing  instability  theories  have  assumed  laminar 
flow  in  both  phases.  Turbulence  has  been  taken  into  account 
by  Jeffreys  ,  Miles  and  Benjamin^,  but  only  to  the 

extent  that  it  determines  the  relevant  properties  of  the  mean 
flow.  Jeffreys  made  use  of  the  principle  that  a  wave  will 
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grow  in  amplitude  if  the  wind  supplies  energy  to  it  at  a  greater 
rate  than  that  of  viscous  dissipation  in  the  water.  He  neglected 
inertial  effects  and  tangential  stress  due  to  the  wind  and  con¬ 
sidered  only  the  pressure  difference  on  the  windward  and 
leeward  sides  of  the  waves  caused  by  separation  of  the  air  stream 
on  the  leeward  side.  To  express  the  magnitude  of  the  pressure 
difference  Jeffreys  introduced  a  sheltering  coefficient,  S,  the 
numerical  value  of  which  he  calculated  from  experimental  obser¬ 
vation  of  the  critical  wind  velocity. 

Jef f reys 1 s  criterion  for  the  generation  of  waves  is: 
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However,  experimental  measurements  of  S  using  fixed 
wave  models  have  shown  it  to  be  an  order  of  magnitude  smaller 
than  that  postulated  by  Jeffreys,  and  Ursell^^^  drew  the 
conclusion  that  the  sheltering  theory  in  its  simplest  form 
is  insufficient  to  explain  the  growth  of  waves.  However,  Miles 
has  removed  the  doubt  concerning  the  value  of  S  by  showing 
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theoretically  that  for  a  moving  wave  system  a  sufficiently 
large  value  of  S  may  be  attained  even  in  the  absence  of  flow 
separation.  Miles  concludes  that  his  model  is  an  important 
adjunct  to  Phillips's  theory.  Benjamin^)  calculated  the 
normal  and  tangential  stresses  on  a  wavy  surface  caused  by 
the  turbulent  flow  of  air  and  showed  that  they  are  distributed 
in  much  the  same  way  as  if  the  leeward  slopes  of  the  waves 
were  sheltered. 

Whatever  mechanism  is  responsible  for  the  initiation 
of  water  waves,  once  the  waves  are  formed  energy  is  transferred 
to  the  waves  by  several  mechanisms  and  any  simplified  theory 
is  unlikely  to  predict  exactly  the  characteristics  of  the 
waves  observed  in  nature. 

C.  THE  STRATIFIED  FLOW  OF  MISCIBLE  FLUIDS 

In  addition  to  the  case  of  wind  blowing  over  water 
other  naturally  occurring  stratified  systems  are  important  in 
meteorology  and  oceanography.  These  may  be  divided  into  two 
groups,  depending  on  whether  the  stratification  is  due  to  a 
continuous  variation  in  density  or  whether  it  is  caused  by 
fluid  layers  of  essentially  different  density.  An  example  of 
the  former  type  is  an  atmosphere  in  which  the  density  decreases 
continuously  with  height,  and  an  example  of  the  second  type  is 
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found  in  river  estuaries  when  the  river  efflux  tends  to  flow 
over  the  salt  water  of  the  ocean.  In  all  these  naturally 
occurring  phenomena  the  flow  of  the  different  layers  may  be 
co-current  or  counter- current.  In  the  case  of  a  river  estuary, 
for  example,  flow  will  tend  to  be  co- current  while  the  tide 
is  ebbing  and  counter- current  while  the  tide  is  rising. 

The  important  difference  between  these  stratified 
systems  and  the  air-water  system  is  that,  whereas  air  and  water 
are  immiscible,  the  layers  in  the  stratified  systems  cited 
above  are  miscible  with  each  other.  Thus,  while  at  low  differ¬ 
ential  velocities  the  layers  remain  distinct,  at  higher  velocities 
the  interface  becomes  unstable  and  waves  develop  in  much  the 
same  way  as  for  the  air-water  system  with  the  subsequent  diffusion 
of  one  layer  into  the  other. 

Motions  in  a  continuously  stratified  fluid  have  been 
investigated  theoretically  by  Taylor  ,  Yih^^'^^,  Drazin 
and  Howard  (22)  ^  and  Eliassen  et  al(^3)  an(^  both  theoretically 
and  experimentally  by  Long (^4) ^ 

Counter- current  flow  of  distinct  layers  has  been  studied 
by  Graebel^^^,  Mittendorf  (^9)  an(j  Macagno  and  Rouse  (^7)  ancj 
is  relatively  unstable  because  of  the  inflexion  point  which 
occurs  in  the  velocity  profile  at,  or  near,  the  interface.  The 
co- current  flow  of  two  or  more  fluid  layers  over  an  obstacle, 
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representing  a  mountain  barrier  to  the  movement  of  air  currents, 
has  been  studied  both  theoretically  and  experimentally  by 
Long »66).  T^e  motion  of  two  liquid  layers  towards  an 
orifice,  as  could  be  encountered  in  the  withdrawal  of  water 
from  the  ocean,  has  been  investigated  theoretically  by  Craya^-*-^) 
and  experimentally  by  Gariel^-*).  Esch(30,31)  considered 
the  stability  of  a  shear  layer  between  unbounded  parallel 
streams,  and  the  flow  of  a  liquid  with  a  free  surface  over  a 
slightly  heavier  liquid.  He  found  theoretically,  using  small 
disturbance  principles,  that  the  former  case  is  always  poten¬ 
tially  unstable  and  the  latter  may  be  stable,  depending  primarily 
upon  the  depth  and  velocity  of  the  upper  fluid.  Presumably 
the  presence  of  the  free  surface  has  a  stabilizing  effect. 

Investigations  of  essentially  co-current  flow  which  are 
of  particular  relevance  in  the  present  study  are  those  of 
Taylor  Keulegan  and  Ippen  and  Harlernan  ^ ^  .  All 

were  concerned  with  the  development  of  interfacial  instability 
and  the  subsequent  mixing  in  systems  in  which  one  fluid  flowed 
past  a  second  layer  of  miscible  fluid,  which  was  stationary 
except  for  the  movement  caused  by  the  drag  of  the  moving  phase. 
Lofquist^^)  has  recently  extended  Keulegan' s  work,  but  without 
further  clarification  of  the  criterion  for  interfacial 


instability. 
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Taylor  observed  the  flow  of  fresh  water  over  a  layer 
of  brine  in  a  closed  rectangular  conduit  and  described  how 
initially  the  interface  was  smooth  and  subsequently,  with 
increase  of  the  water  flow  rate,  wave  motions  and  mixing 
developed  at  the  interface.  Taylor  did  not  report  quantitatively 
the  conditions  under  which  the  instability  occurred. 

Keulegan  performed  experiments  similar  to  those  of 
Taylor  in  conduits  having  cross  sections  of  2  by  4  cms. ,  4 
by  3  cms.  and  11.3  by  23.5  cms.,  the  greater  dimension  being 
vertical  in  each  case.  The  fluids  used  were  water  and  sugar 
solution  and  the  depths  of  the  two  layers  were  approximately 
equal.  Experiments  were  conducted  in  each  of  the  three  conduits 
to  determine  the  effect  of  the  Reynolds  number  of  the  flowing 
water  layer  on  the  instability  of  the  interface,  while  tests 
were  conducted  in  the  medium  sized  conduit  to  determine  the 
effect  of  the  viscosity  of  the  lower  or  sugar  solution  layer 
on  the  instability.  He  recorded  the  maximum  velocity  of  the 
water  layer,  calculated  by  dividing  the  volumetric  flow  of 
water  by  the  area  of  cross-section  of  the  water  layer,  at  which 
the  interface  was  smooth  and  the  velocity  at  which  the  first 
signs  of  mixing  appeared.  The  mean  of  these  two  velocities 
was  defined  as  the  critical  velocity,  V'  . 
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Keulegan  derived  a  criterion  for  instability  of  the 


interface  as: 
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Experiments  in  the  medium  sized  conduit  confirmed  the 
criterion  to  the  extent  that  V'wc  was  approximately  proportional 
to  (to 1/3 .  Experiments  in  the  three  conduits  enabled  the 

b 

value  of  ©  to  be  determined  for  both  laminar  and  turbulent  flow 
of  the  upper  fluid.  For  turbulent  flow  ©  =  0.178  and  for 
laminar  flow  ©  --  0.127. 

Keulegan  pointed  out  that  his  criterion  has  the  same 
form  as  that  of  Jef f reys  and,  further,  the  value  of  © 

determined  from  Jeffreys ' s  observations  of  the  initiation  of 
waves  on  a  water  surface  is  0.215  -  not  very  different  in 
magnitude  from  Keulegan' s  turbulent  value  of  0.178. 

Ippen  and  Harleman  performed  similar  experiments  to 
those  of  Keulegan,  except  brine  under flowed  a  layer  of  fresh 
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water.  They  measured  the  critical  velocity  for  conditions  in 
which  the  brine  was  in  laminar  flow  at  Reynolds  numbers  less 
than  those  in  Keulegan's  experiments.  For  this  regime  they 
found,  as  they  predicted  theoretically,  that  ©  was  a  linear 
function  of  the  reciprocal  of  the  Reynolds  number  and  not  a 
constant  as  Keulegan  had  predicted.  However,  Keulegan  obtained 
few  experimental  data  in  the  laminar  regime  and  these  do 
agree  closely  with  the  results  of  Ippen  and  Harleman.  Most 
of  Keulegan's  data  were  obtained  for  the  transition  and  tur¬ 
bulent  regimes  and  there  is  little  doubt  that  in  the  turbulent 
regime  the  constancy  of  ©  was  verified. 

Since  both  Keulegan,  and  Ippen  and  Harleman,  observed 
waves  when  both  phases  were  in  laminar  motion,  it  would  appear 
that  the  waves  must  have  been  caused  by  a  fundamental  instability 
of  the  mean  flow  and  not  by  turbulent  pressure  fluctuations 
in  one  of  the  phases. 

D.  CO- CURRENT  TWO- PHASE  FLOW  THROUGH  CLOSED,  HORIZONTAL  CONDUITS 

1.  The  flow  of  gas  and  licruid  through  circular  pipes 

The  flow  of  two  phases  through  closed,  horizontal  conduits 
is  encountered  in  many  engineering  situations,  such  as  pipe¬ 
line  transportation,  tubular  reactors,  and  film  coolers.  In 
this  section  attention  is  focused  on  gas- liquid  flow  through 
circular,  horizontal  pipes. 
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Little  doubt  remains  concerning  the  nature  of  the  flow 
patterns  which  occur  in  co- current  gas- liquid  horizontal 
pipeline  flow.  Govier  and  Omer  (^8)  and  Hoogendoorn  (^4)  have 
covered  wide  ranges  of  individual  flow  rates  and  input  ratios. 
Govier  and  Omer  observed  the  flow  patterns  established  by  air 
and  water  flowing  in  a  1-inch  pipe,  while  Hoogendoorn  was 
concerned  with  air  and  water,  and  air  and  oil,  flowing  in 
pipes  with  diameters  ranging  from  1  to  5.5  inches.  At  low 
flow  rates  of  both  phases  stratified  flow  with  a  smooth  inter¬ 
face  takes  place.  If  the  flow  of  liquid  is  maintained  constant 
and  the  gas  flow  rate  increased  the  interface  becomes  unstable 
and  waves  develop.  Further  increase  in  the  gas  flow  rate 
causes  the  liquid  to  climb  up  the  pipe  wall  and  eventually 
"annular  flow"  is  established  when  a  liquid  annulus  adheres 
to  the  pipe  wall  and  completely  surrounds  the  gas  core. 

If  a  low  gas  flow  rate  is  initially  combined  with  a 
relatively  large  liquid  flow  rate  then  "bubble  flow"  -  in 
which  bubbles  of  gas  move  along  the  top  of  the  pipe  -  occurs. 
Subsequent  increase  in  the  gas  flow  rate  does  not  establish 
stratified  flow  -  but  rather  "plug  flow"  wherein  larger  gas 
bubbles  or  plugs  move  along  the  upper  part  of  the  pipe  and 
subsequently  "slug  flow"  in  which  splashes  or  slugs  of  liquid 
occasionally  pass  through  the  pipe  with  a  higher  velocity 
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than  the  bulk  of  the  liquid.  Slug  flow  may  be  considered  as 
an  extension  of  stratified  flow  in  which  the  waves  have  become 
very  large  compared  with  the  pipe  diameter  and  liquid  is  torn 
from  the  crests. 

Stratified  flow,  with  or  without  waves,  is,  therefore, 
but  one  of  the  flow  patterns  encountered  in  the  simultaneous 
flow  of  gas  and  liquid  through  a  circular  pipe.  Since  the  gas 
is  compressible  its  density  must  decrease  as  the  pressure 
decreases  along  the  pipeline.  The  consequent  acceleration  of 
the  gas  causes  an  increasing  drag  upon  the  liquid  surface  and 
the  liquid  also  accelerates.  This  manifests  itself,  for  constant 
gas  and  liquid  mass  rates  of  flow,  in  an  interfacial  gradient 
and  unequal  pressure  gradients  in  the  gas  and  liquid  phases. 

No  analytical  solution  exists  for  the  stratified  flow  of  two 
fluids  when  one  of  them  is  compressible. 

Govier  and  Omer (^8)  in  reviewing  the  literature  on  the 
flow  of  gas- liquid  mixtures  through  pipes  concluded  that  the 
nature  of  the  flow  pattern  greatly  influences  the  pressure 
gradient.  They  suggested,  therefore,  that  either  a  separate 
correlation  is  required  for  each  flow  pattern  or  in  the  case 
of  an  overall  correlation,  such  a  correlation  must  include 
all  the  variables  which  determine  the  flow  pattern.  They  also 
noted  that,  while  the  effects  of  the  individual  flow  rates  on 
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the  flow  pattern  are  reasonably  well  known,  the  effects  of 


the  densities  and  viscosities  of  the  fluids  and  the  pipe 


diameter  are  not. 

Lockhart  and  Martinelli 


proposed  an  overall  correlation, 


with  complete  disregard  for  the  flow  patterns,  in  terms  of  two 
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=  pressure  gradient  for  gas  flowing  alone 


pressure  gradient  for  liquid  flowing  alone 


and  the  four  flow  regimes:  turbulent- turbulent ,  turbulent- 
laminar,  laminar- turbulent ,  and  laminar- laminar.  A  phase  was 
assumed  to  be  in  turbulent  flow  if  the  superficial  Reynolds 
number  for  the  phase  exceeded  2000  and  to  be  in  laminar  flow 
if  the  Reynolds  number  was  less  than  1000.  Four  curves  were 
established,  one  for  each  regime,  with  which  their  experimental 
two- phase  pressure  gradients  agreed  to  within  about  -  55  percent. 
It  has  since  been  found  by  Gazley^"^  ,  Hoogendoorn  and 
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Govier  and  Omer (^8)  that  the  Lockhart-Martinelli  correlation 
considerably  and  consistently  over-estimates  the  two  phase 
pressure  gradient  for  stratified  and  wave  flow  in  all  four 
regimes.  A  basic  assumption  made  by  Lockhart  and  Martinelli 
was  that  the  pressure  gradients  in  both  phases  were  equal. 

That  this  is  definitely  not  so  for  stratified  flow  has  given 
some  investigators,  for  instance,  Gazley^^  a  ready  explanation 
for  the  discrepancy.  However  it  would  appear  that  the  basic 
assumption  of  equality  of  pressure  gradients  in  the  two  phases 
is  unlikely  to  hold  for  any  of  the  flow  patterns  and  that 
the  real  reason  for  the  discrepancy  is  that  the  correlation 
is  inadequate  to  cover  the  different  flow  patterns.  This  has 
been  strikingly  demonstrated  by  Hoogendoorn who  was  able, 
by  making  suitable  modifications  to  the  entry  section,  to 
obtain  either  stratified  flow  with  waves  or  slug  flow, 
other  conditions  being  the  same.  Pressure  gradients  for  slug 
flow  were  approximately  twice  those  for  the  corresponding 
wave  flow. 

Gazley(^)  investigated  the  co-current,  stratified  flow 
of  air  and  water  in  a  horizontal  2- inch  diameter,  2 2- foot 
long  pipe.  No  part  of  the  pipe  was  transparent  and  the 
details  of  the  flow  patterns  could  not  be  observed.  The  water 
phase  was  allowed  to  fall  freely  from  the  end  of  the  pipe  and 
because  of  this  the  inter facial  and  pressure  gradients  reported. 
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especially  at  low  flow  rates,  cannot  be  considered  to  be 
representative  of  those  to  be  expected  in  a  long  pipeline. 
Nevertheless,  certain  qualitative  effects  are  of  definite 
interest. 

Gazley  found  by  measuring  pressure  gradients  that  the 
presence  of  a  moving  layer  of  water  delayed  the  transition  to 
turbulence  in  the  gas  phase.  Whereas  for  the  flow  of  air 
alone  the  transition  from  laminar  to  turbulent  flow  occurred 
over  the  range  of  air  Reynolds  number  2000  to  4000,  for  the 
two  phase  system  the  transition  took  place  in  the  range  of 
superficial  air  Reynolds  number  3000  to  5000.  Conversely 
Gazley  found  that  the  presence  of  a  turbulent  air  layer  caused 
the  transition  to  turbulence  in  the  water  phase  to  occur  at 
slightly  lower  Reynolds  numbers  than  for  single  phase  flow  of 
water. 

The  point  at  which  the  interface  became  wavy  for  a 
constant  liquid  flow  rate  was  detected  by  the  presence  of  a 
change  in  slope  in  the  curve  representing  the  relationship 
between  the  two-phase  pressure  gradient  and  superficial  air 
velocity.  While  the  slope  was  1.8  with  both  phases  in 
turbulent  motion  and  no  interfacial  waves,  the  slope  increased 
to  2.75  when  waves  were  present.  It  was  found  that  the  onset 
of  a  wavy  interface  did  not  occur  at  a  definite  air  Reynolds 
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number  and  depended  on  the  water  flow  rate.  The  higher  the 
water  flow  rate,  the  lower  the  air  Reynolds  number  at  which 
the  waves  first  appeared.  However,  the  air  was  always  either 
in  transitional  or  turbulent  flow  when  waves  were  detected. 

Mologin^0^  determined  the  boundaries  of  the  wave  flow 
pattern  for  the  flow  of  air  and  water  in  pipes  ranging  in  diameter 
from  1  to  4  inches.  He  determined  the  period  and  amplitude 
of  the  waves  with  the  aid  of  a  motion  picture  camera  and 
presented  empirical  formulae  for  these  quantities.  No  pressure 
gradients  were  reported,  nor  were  the  transitions  from  laminar 
to  turbulent  flow  in  each  phase  indicated. 

2.  The  flow  of  gas  and  liquid  through  rectangular  conduits 

Theoretical  and  experimental  studies  of  the  stratified, 
co-current  flow  of  gas  and  liquid  in  a  rectangular  conduit 
have  been  undertaken  to  determine  the  basic  nature  of  the 
interaction  between  the  phases.  Stratified  flow  in  a  rectangular 
conduit  is  a  simpler  situation  than  stratified  flow  in  a  cir¬ 
cular  pipe.  In  a  rectangular  conduit  the  depths  of  both  phases 
are  constant  across  the  width  of  the  conduit  and,  provided  the 
ratio  of  width  to  depth  is  large,  the  velocity  distribution 

is  two-dimensional  in  the  central  region  of  the  conduit. 

(52)  (49) 

Keulegan v  '  and  Johnson  and  Rice '  investigated  the 

flow  of  air  over  a  "stationary"  body  of  water,  simulating  wind 
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blowing  over  the  ocean.  Keulegan  was  interested  in  the  inter¬ 
facial  gradients  developed  and  Johnson  and  Rice  measured 
wave  characteristics  at  various  points  along  the  conduit. 

The  significant  heights  and  periods  were  related  to  the  air 
velocity  and  fetch. 

Hanratty  and  Engen^0)  studied  the  interaction  between 
a  turbulent  air  stream  and  a  water  layer  flowing  co- currently 
in  a  12-inch  by  1-inch  rectangular  conduit.  Five  flow  patterns 
were  observed:  smooth  interface,  two-dimensional  waves, 
squall  surface  or  three-dimensional  waves,  roll  waves,  and 
dispersed  flow.  Initially  at  low  flow  rates  of  both  phases 
the  interface  was  smooth  and  no  disturbances  were  visible  to 
the  naked  eye.  If  the  liquid  flow  rate  was  maintained  constant 
and  the  gas  flow  steadily  increased  a  series  of  transitions 
took  place.  Two-dimensional  waves  first  appeared  which  were 
long-crested  and  had  a  wavelength  of  approximately  1  cm. , 
followed  by  three-dimensional  waves  in  which  the  lengths  of 
the  wave  crests  were  approximately  the  same  as  the  wavelengths. 
Three-dimensional  waves  gave  way  to  roll  waves,  which  had 
irregular  shaped  crests  running  across  the  width  of  the  conduit. 
In  this  flow  pattern  part  of  the  liquid  was  picked  up  and 
carried  over  the  liquid  surface.  Finally  dispersed  flow,  in 
which  droplets  of  liquid  were  removed  from  the  liquid  phase 
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and  carried  in  the  gas  stream,  was  observed.  As  would  be 
expected  these  flow  patterns  are  similar  to  those  observed  in 
pipe  flow  when  the  liquid  flow  rate  is  relatively  small. 

Hanratty  and  Engen  made  several  significant  observations. 
They  found  that  the  transition  from  a  smooth  interface  to  one 
with  two-dimensional  waves  did  not  necessarily  coincide  with 
the  transition  from  laminar  to  turbulent  flow  in  the  gas  phase. 
In  fact  they  noted  that  at  low  liquid  flow  rates  a  smooth 
interface  coexisted  with  a  turbulent  gas  phase  and,  at  high 
liquid  flow  rates  two-dimensional  waves  were  formed  in  the 
presence  of  a  laminar  gas  phase.  Thin  liquid  layers  were 
therefore  more  stable  than  thick  ones.  This  was  found  to  be 
so  for  all  the  transitions  in  surface  structure.  Further, 
they  found  the  presence  of  the  interfacial  waves  distorted 
the  gas  phase  velocity  profile  in  that  the  maximum  in  the 
profile  was  moved  away  from  the  interface.  This  effect,  which 
was  negligibly  small  for  two-dimensional  waves,  increased  with 
increasing  roughness  of  the  interface.  The  position  of  the 
maximum  in  the  velocity  profile,  together  with  the  depth  of 
the  water  layer  and  the  pressure  gradient  enabled  the  inter¬ 


facial  shear  stress  to  be  evaluated. 
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Ellis  and  Gay  ^ ^  performed  similar  experiments  to 
Hanratty  and  Engen  and  determined  the  interfacial  shear  stress 
from  a  knowledge  of  the  gas  phase  velocity  profile  and  the 
liquid  depth  only.  Ellis  and  Gay  compared  their  results  with 
those  of  Hanratty  and  Engen  and  found  good  agreement. 

Hanratty  and  Hershman  (^1)  derived  a  criterion  for  the 
initiation  of  roll  waves  on  a  liquid  layer  flowing  co- currently 
with  gas  in  a  closed  rectangular  conduit  by  extending  Jeffreys 's 
theory .  Excellent  agreement  between  theory  and  experiment 
was  obtained  in  that  the  theory  correctly  predicted  both  the 
trend  and  numerical  values  of  the  gas  phase  Reynolds  number 
at  which  roll  waves  first  appeared  for  air  and  water  flowing 
in  a  12-inch  by  1-inch  conduit.  If  the  effect  of  surface 
tension  is  neglected,  then  as  the  liquid  phase  Reynolds  number 
becomes  very  small,  their  criterion  becomes: 

V' 
gc 

actual  time  average  velocity  in  gas  phase 
gas  velocity  profile  factor  determined  from 
experimental  measurements 
depth  of  the  conduit 
density  ratio  (gas/liquid) 
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Lilleleht  and  Hanratty developed  an  experimental 
technique,  based  on  the  partial  absorption  of  a  light  beam 
by  a  dye  dissolved  in  the  water  phase,  to  determine  wave 
characteristics  during  the  co- current  flow  of  air  and  water 
in  a  rectangular  conduit.  At  a  constant  liquid  Reynolds 
number  the  root  mean  square  wave  height  was  found  to  increase 
rapidly  with  increasing  gas  phase  Reynolds  number  for  two- 
dimensional  waves  and  then  pass  through  a  maximum  for  three- 
dimensional  waves.  The  three-dimensional  waves  were  explored 
in  greater  detail  and  it  was  found  that  a  Gaussian  model 
approximately  described  the  wave  structure.  The  interfacial 
shear  stress  was  correlated  with  the  root  mean  square  wave 
height. 

(91) 

Van  Rossurn'  '  observed  the  flow  of  air  over  a  thin 
layer  of  water  in  a  6- inch  square  duct.  He  was  interested  in 
the  transition  from  a  smooth  interface  to  one  with  waves  and 
from  a  wavy  interface  to  dispersed  flow.  The  critical  air 
velocity  necessary  to  produce  waves  decreased  with  increasing 
film  thickness,  confirming  the  observations  of  Hanratty  and 
Engen^®)  and  Gazley(^),  and  tended  to  a  constant  value  for 
thick  water  layers.  This  constant  value  was  about  twice  the 
value  predicted  by  Jeffreys's  theory .  The  average  thickness 
of  the  water  layer  with  waves  present  at  the  interface  was 
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found  to  be  only  about  50  percent  of  the  thickness  predicted 
by  laminar  flow  theory  assuming  a  smooth  interface. 

The  fact  that  the  initiation  of  interfacial  waves 
during  the  flow  of  air  and  water  through  rectangular  channels 
does  not  coincide  with  the  transition  to  turbulent  flow  in 
the  gas  phase  leads  to  the  same  conclusion  drawn  from  the 
results  of  experiments  with  miscible  liquids,  i.e.  the  waves 
are  caused  by  a  fundamental  instability  of  the  mean  flow. 

3.  The  flow  of  two  immiscible  liquids  through  circular  pipes 

It  has  previously  been  mentioned  that  the  pressure 
gradient  and  power  requirement  necessary  for  the  flow  of  a 
viscous  liquid  at  a  given  flow  rate  in  a  horizontal,  circular 
pipe  may  be  substantially  reduced  by  the  presence  of  a  less 
viscous  liquid.  In  order  to  estimate  the  effectiveness  of 
different  flow  configurations  Russell  and  Charles  analysed 

theoretically  the  laminar  flow  of  two  incompressible  fluids 
flowing  stratified  between  infinitely  wide  parallel  plates 
and  concentrically  in  a  circular  pipe.  It  was  apparent  that 
concentric  flow,  in  which  the  less  viscous  fluid  occupies  the 
region  of  high  shear  rate  next  to  the  pipe  wall,  is  much  more 
effective  in  reducing  the  pressure  gradient  than  stratified 
flow  where  the  less  viscous  liquid  is  in  contact  with  only 
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part  of  the  conduit  boundary.  For  both  flow  patterns  the 
pressure  gradient  reduction  factor  -  defined  as  the  ratio 
of  the  pressure  gradient  for  the  flow  of  the  more  viscous 
fluid  at  a  given  flow  rate  to  the  pressure  gradient  for  the 
two  phase  system  with  the  more  viscous  fluid  flowing  at  the 
same  rate  -  was  found  to  be  a  function  of  only  the  viscosity 
ratio  and  the  input  ratio. 

For  concentric  flow  Russell  and  Charles  found  that  the 
maximum  pressure  gradient  reduction  factor  was  proportional 
to  the  viscosity  ratio  at  high  viscosity  ratios.  At  a  vis¬ 
cosity  ratio  of  1000  the  maximum  pressure  gradient  reduction 

(12) 

factor  was  predicted  to  be  500.  Charles,  Govier  and  Hodgson' 
reported  the  results  of  an  experimental  investigation  of  the 
co-current  pipeline  flow  of  two  liquids,  oil  and  water,  having 
the  same  densities.  The  liquids  were  introduced  concentrically 
into  the  pipe  and  concentric  flow,  with  the  oil  flowing  inside 
the  water,  was  found  to  be  stable  over  a  wide  range  of  flow 
rates.  Three  oils  having  viscosities  of  6.29,  16.8  and  65 
centipoise  were  used  in  the  experiments.  For  each  oil  the 
maximum  pressure  gradient  reduction  factor  observed  was  some¬ 
what  less  than  that  predicted  by  theory  because  concentric 
flow  was  not  stable  at  the  input  ratio  predicted  by  theory  to 
give  the  maximum  benef iciation,  the  oil  forming  slugs  moving 
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in  a  continuous  water  phase.  Nevertheless  the  experimentally 
determined  maximum  pressure  gradient  reduction  factors  were 
considerable  -  being  1.7  for  the  6.29  centipoise  oil,  6.0  for 
the  16.8  centipoise  oil  and  10.0  for  the  65  centipoise  oil. 

The  commercial  application  of  the  concentric  configuration  to 
reduce  pressure  gradients  during  the  pipeline  transportation 
has  been  outlined  in  patents  by  Isaacs  and  Speed  ,  Clark 
and  Shapiro  and  Chilton  and  Handley  .  However,  it  is 

unlikely  that  this  extremely  desirable  configuration  will 
find  widespread  use  because,  in  the  presence  of  even  a  small 
density  differential,  the  fluids  tend  to  stratify. 

The  results  obtained  by  Russell  and  Charles  for  stratified 
flow  between  parallel  plates  indicated  that  the  maximum  pressure 
gradient  reduction  factor  was  not  proportional  to  the  viscosity 
ratio  at  high  viscosity  ratios  as  for  concentric  flow  but 
appeared  to  approach  a  constant  value  at  high  viscosity  ratios. 

At  a  viscosity  ratio  of  1000  the  maximum  pressure  gradient 
reduction  factor  was  predicted  to  be  approximately  3.7,  compared 
with  500  for  concentric  flow. 

The  laminar,  stratified  flow  of  two  incompressible  fluids 
in  a  circular  pipe  has  since  been  treated  numerically  by 
Charles  and  Redberger^*^  and  by  Gemmell  and  Epstein  using 

slightly  different  numerical  methods.  The  analytical  solution 
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has  now  also  been  obtained  by  Bentwich ^  but  his  solution  has 
not  yet  been  evaluated  and  compared  with  the  numerical  results. 

The  numerical  results  of  Charles  and  Redberger,  and 
Gemmell  and  Epstein  agree  well  with  each  other  in  that  the 
pressure  gradient  reduction  factors  predicted  by  the  two 
methods  agree  within  about  7  percent,  those  of  Gemmell  and 
Epstein  being  consistently  greater  than  those  of  Charles  and 
Redberger.  It  is  doubtful  which  results  are  more  accurate; 
Charles  and  Redberger  used  considerably  smaller  distance 
increments  in  the  finite  difference  approximation  to  the  basic 
differential  equation,  while  Gemmell  and  Epstein  used  some¬ 
what  better  finite  difference  approximations  to  the  boundary 
conditions.  As  for  stratified  flow  between  parallel  plates, 
the  pressure  gradient  reduction  factor  is  a  function  of  the 
viscosity  ratio  and  input  ratio  only  and  the  maximum  pressure 
gradient  reduction  factor  tends  to  a  constant  value  at  high 
viscosity  ratios.  At  a  viscosity  ratio  of  1000  Charles  and 
Redberger  predicted  a  maximum  pressure  gradient  reduction 
factor  of  1.31  and  Gemmell  and  Epstein  a  value  of  1.41. 
Stratified  flow  in  a  circular  pipe  is,  therefore,  the  least 
attractive  configuration  from  the  point  of  view  of  pressure 
gradient  reduction.  This  is  a  direct  consequence  of  the 
proportion  of  the  conduit  boundary  which  is  wet  by  the  less 
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viscous  fluid.  The  addition  of  a  small  amount  of  the  less 
viscous  fluid  wets  the  complete  boundary  in  concentric  flow, 
half  the  boundary  for  stratified  flow  between  parallel  plates 
and  only  a  small  proportion  of  the  boundary  for  stratified 
flow  in  a  circular  pipe. 

In  the  above  theoretical  studies  the  fluids  were  assumed 
to  be  incompressible  and  have  Newtonian  characteristics.  The 
motion  was  assumed  to  be  laminar  with  a  smooth  interface. 
Integration  to  obtain  the  velocity  distribution  and  volumetric 
flow  rates  was  carried  out  using  the  appropriate  boundary 
conditions,  which  were 

(1)  the  velocity  of  the  more  viscous  liquid  in  contact 
with  the  boundary  was  zero 

(2)  the  velocity  of  the  less  viscous  liquid  in  contact 
with  the  boundary  was  zero 

(3)  the  velocities  of  both  phases  were  equal  at  the 
interface,  i.e.  there  was  no  slip  between  the  phases 

(4)  shear  stress  was  continuous  across  the  interface 

These  conditions  would  appear  to  be  quite  adequate 

although  any  effects  due  to  the  meniscus  at  the  pipe  wall  in 
stratified  flow  and  anomalous  interfacial  viscosity,  were 
neglected. 
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Several  experimental  investigations  have  produced 
data  which  may  be  compared  with  the  numerical  results  obtained 
for  stratified  flow  in  a  circular  pipe.  Russell,  Hodgson  and 
Govier^9)  reported  pressure  gradients  for  the  co-current 
flow  of  an  oil  having  a  viscosity  of  18  centipoise  and  water 
in  a  1-inch  diameter  pipe.  At  high  input  ratios  of  water  to 
oil  the  oil  formed  itself  into  bubbles  which  flowed  along  the 
top  of  the  pipe  in  a  continuous  water  phase.  Decrease  in  the 
water  to  oil  input  ratio  resulted  in  stratified  flow  with  a 
smooth  interface.  At  higher  flow  rates  the  interface  became 
wavy,  and  subsequently  a  mixed  type  of  flow  was  observed.  No 
attempt  was  made  to  describe  in  detail  the  characteristics  of 
the  wavy  interface.  The  pressure  gradient  reduction  factors 
computed  from  their  pressure  gradient  measurements  for  conditions 
under  which  the  flow  was  stratified  with  a  smooth  interface 
are  plotted  in  Figure  1,  along  with  the  predicted  curves  of 
Charles  and  Redberger,  and  Gemmell  and  Epstein.  The  experi¬ 
mental  data  scatter  approximately  t  30  percent  about  the 
predicted  curves  and  in  addition,  while  theory  predicts  that 
for  a  given  viscosity  ratio  the  pressure  gradient  reduction 
factor  is  a  function  of  the  input  ratio  or  input  oil  fraction 
only,  the  experimental  data  appear  to  be  a  function  of  the 
water  Reynolds  number  as  well.  Further,  at  input  fractions 
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FIGURE  1  -  COMPARISON  OF  THEORETICAL  AND  EXPERIMENTAL  PRESSURE 
GRADIENTS  FOR  OIL-WATER  LAMINAR  STRATIFIED  FLOW  IN  A  CIRCULAR 

PIPE. 
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where  theory  predicts  the  pressure  gradient  reduction  factor 
to  pass  through  a  maximum  the  experimental  data  show  distinct 
minima.  Govier  and  Omer  ^8)  measured  pressure  gradients  for 
the  stratified  flow  of  air  and  water.  Data  for  which  both 
phases  were  in  laminar  motion  are  plotted  in  Figure  2,  along 
with  the  curve  predicted  by  Charles  and  Redberger.  Under 
the  low  pressure  gradients  prevailing  when  both  phases  were 
in  laminar  flow  the  air  can  be  considered  to  be  incompressible. 

It  is  apparent  that  the  experimental  reduction  factors  are 
about  50  percent  of  the  predicted  values.  Charles 
reported  pressure  gradient  reduction  factors  for  the  flow  of 
viscous  crude  oil  with  water  in  both  1-  and  2-^2 -inch  diameter 
pipes.  While  a  significant  density  differential  existed  between 
the  oil  and  water  and  the  flow  was  assumed  to  be  stratified, 
it  was  not  possible  to  observe  the  flow  pattern.  The  measured 
pressure  gradient  reduction  factors  were  an  order  of  magnitude 
greater  than  those  predicted  by  theory. 

Assuming  that  the  experimental  flows  were  fully  developed, 
the  descrepancies  between  the  experimental  and  theoretical 
results  may  be  attributed  to  the  fact  that  the  actual  experi¬ 
mental  situations  were  not  exact  reproductions  of  the  theoretical 
model.  In  the  1-inch  diameter  pipes  used  by  Russell,  Hodgson 
and  Govier,  and  Govier  and  Omer  the  meniscus  effect  at  the  pipe 
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WATER  INPUT  VOLUME  FRACTION 


FIGURE  2  -  COMPARISON  OF  THEORETICAL  AND  EXPERIMENTAL  PRESSURE 
GRADIENTS  FOR  AIR-WATER  LAMINAR  STRATIFIED  FLOW  IN  A  CIRCULAR 

PIPE. 
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wall  may  well  have  had  a  significant  effect  on  the  pressure 
gradient,  while  the  extremely  high  pressure  gradient  reduction 
factors  obtained  by  Charles  with  the  crude  oil-water  system 
may  be  explained  by  the  fact  that  a  certain  amount  of  water 
adhered  to  the  upper  part  of  the  pipe  wall  and  aided  the  flow 
of  the  oil. 

E.  THE  EFFECT  OF  ADDITIVES  ON  THE  STABILITY  OF  A  FLUID-FLUID 
INTERFACE 

The  presence  of  a  surface  active  agent  or  a  solute 
soluble  in  both  fluids  may  have  a  profound  effect  on  the 
stability  of  a  fluid- fluid  interface. 

The  effect  of  the  presence  of  a  monolayer  of  surface 
active  agent  on  the  appearance  of  water  waves  generated  by 
wind  has  been  discussed  by  Davies  .  The  monolayer  effectively 
eliminates  the  small  ripples  which  are  superimposed  on  the 
longer  wavelength  gravity  waves.  This  causes  a  reduction  in 
the  drag  exerted  by  the  wind  and,  at  the  same  time,  greatly 
reduces  the  amount  of  energy  which  is  dissipated  in  the  wave 
motions.  The  combined  effect  is  that  more  energy  is  converted 
into  kinetic  energy  associated  with  the  bulk  movement  of  the 
water.  That  this  is  so  was  strikingly  demonstrated  by  Davies 
when  a  slick  of  hexadeconal  was  found  to  accelerate  across  a 
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lake  under  the  action  of  the  wind. 

Hanratty  and  Hershman found  that  the  addition  of 
less  than  0.01  percent  of  sodium  lauryl  sulphate  suppressed 
the  two-  and  three-dimensional  waves  normally  observed  in  the 
flow  of  air  and  water  in  a  rectangular  conduit,  but  did  not 
affect  the  transition  to  roll  waves. 

One  effect  of  the  presence  of  surface  active  agent  is 
to  reduce  the  surface  tension.  Such  an  effect  would  decrease 
the  stability  of  the  interface.  The  experimental  observations 
are  contrary  to  this  and  other  properties  of  surface  films, 
e.g.  anomalous  surface  viscosity  effects,  must  be  important. 

The  presence  of  a  soluble  solute  in,  and  hence  mass 
transfer  between,  two  liquid  phases  causes  local  variations 
in  solute  concentration  and,  therefore,  interfacial  tension. 
These  local  variations  cause  agitation  of  the  interface.  This 
effect,  known  as  the  Marangoni  effect,  has  been  studied  from 
the  hydrodynamical  point  of  view  by  Sternling  and  Scriven^5^ 
and  its  effect  on  mass  transfer  has  recently  been  reported  by 
Sawistowski  and  Goltz  . 
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F.  FORMULATION  OF  THE  PRESENT  INVESTIGATION 

The  following  conclusions  were  drawn  from  the  literature 

survey: 

(1)  The  nature  of  laminar,  transitional  and  turbulent 
flow  of  a  single  phase  is  reasonably  well  understood.  it 
appears  that  fully  developed  laminar  flows,  except  perhaps 
plane  Poiseuille  flow,  are  stable  for  infinitesimal  distur¬ 
bances.  Contemporary  theory  is,  therefore,  inadequate  to  pre¬ 
dict  the  transition  from  laminar  to  turbulent  flow.  Observed 
turbulence  probably  originates  either  in  the  boundary  layer 
present  in  the  entry  section  to  a  conduit  or  from  disturbances 
of  the  fully  developed  flow  which  have  amplitudes  greater 
than  a  threshold  value.  In  a  two-phase  system  the  effect  that 
the  moving  fluid- fluid  interface  has  on  the  transition  from 
laminar  to  turbulent  flow  in  the  phases  is  not  known. 

(2)  There  are  two  basic  approaches  to  the  prediction 
of  instability  at  a  fluid-fluid  interface.  One  is  based  on 
the  effect  of  turbulent  pressure  fluctuations  in  one  of  the 
phases  and  the  other  is  concerned  with  instability  of  the  mean 
flow,  whether  or  not  one  of  the  phases  is  turbulent.  Both 
methods  have  been  applied  with  some  success  to  the  air-water 
system.  However,  instability  analyses  are  exceedingly  complex 
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and  have  led  to  results  which  are  not  in  complete  agreement. 

This  disagreement  is  a  consequence  of  the  fact  that  it  has 
not  yet  been  possible  to  include  all  mechanisms  of  energy 
transfer  in  a  single  theory. 

(3)  The  fact  that  instability  of  the  interface  in  the 
stratified  flow  of  two  miscible  liquids  has  been  observed 

while  both  phases  are  in  laminar  flow  substantiates  the  hypothesis 
that  interfacial  instability  is  not  necessarily  caused  by  tur¬ 
bulent  fluctuations  and  may  be  caused  by  instability  of  the 
mean  flow. 

(4)  The  experimental  investigation  of  co- current  gas- 
liquid  flow  in  circular  pipes,  and  in  rectangular  conduits 
where  attention  has  been  focussed  on  the  stratified  flow  pattern, 
is  reasonably  well  advanced.  On  the  other  hand  the  experimental 
investigation  of  co- current  liquid- liquid  systems  is  not  far 
advanced  despite  the  possible  importance  of  such  systems  in 
pipeline  transportation,  tubular  reactors,  and  liquid- liquid 
extraction  processes.  Further,  the  few  available  data  for 
laminar- laminar  stratified  flow  of  two  immiscible  liquids  in 

a  circular  pipe  by  no  means  verify  theoretical  predictions. 

(5)  The  profound  effects  that  the  presence  of  dissolved 
solutes  and  surface  active  agents  may  have  on  the  stability 

of  a  fluid- fluid  interface  show  that  the  presence  of  such 
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substances  should  be  very  carefully  controlled  in  any  study 
of  the  hydrodynamics  of  two-phase  systems. 

It  was  decided,  therefore,  to  carry  out  an  experimental 
investigation  of  a  co-current  liquid- liquid  system,  and  to 
compare  theoretical  and  experimental  results  for  the  case  when 
both  fluids  were  in  laminar  flow.  A  rectangular  conduit  of 
large  aspect  ratio  was  selected  as  being  most  suitable  for  the 
experiments  because  it  would  result  in  the  depths  of  the  two 
phases  being  constant  across  the  width  of  the  conduit,  the 
velocity  distribution  in  the  central  part  of  the  conduit  being 
approximately  two-dimensional  and,  lastly,  meniscus  effects 
at  the  sides  of  the  conduit  having  a  relatively  small  effect 


on  the  flow. 


II  THEORY 


A.  SINGLE-PHASE  FLOW  THROUGH  A  RECTANGULAR  CONDUIT 

In  the  most  general  case  of  fluid  flow  the  variables 
velocity,  temperature,  pressure  and  physical  properties  of 
the  fluid  are  functions  of  the  space  co-ordinates  x,  y  and  z 
and  the  time  t.  The  basic  differential  equations  of  fluid 
flow  are  the  continuity  equation  which  is  based  on  the  law 
of  conservation  of  mass,  the  momentum  equation  which  is  based 
on  Newton's  second  law  of  motion,  and  the  energy  equation 
which  is  based  on  the  law  of  conservation  of  energy.  In 
addition,  two  other  relationships  are  usually  available  -  the 
equation  of  state  of  the  fluid  and,  in  the  case  of  Newtonian 
fluids,  the  relation  between  shear  stress  and  shear  rate  or 
velocity  gradient.  For  non- Newtonian  fluids  the  shear  stress- 
shear  rate  relationship  can  usually  be  approximated  by  a 
mathematical  model,  at  least  over  a  limited  range  of  shear 
rate.  The  most  general  forms  of  these  equations  have  been 
discussed  and  summarized  by  Knudsen  and  Katz  . 

When  the  fluid  is  incompressible  and  Newtonian  and  the 
flow  is  horizontal,  steady,  and  fully  developed  in  a  conduit 


of  constant  cross-sectional  area  then,  in  the  absence  of  heat 
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interchange  with  the  surroundings  and  body  forces  other  than 
gravity,  the  energy  equation  reduces  to  the  statement  that 
the  pressure  loss  along  the  conduit  is  completely  accounted 
for  by  viscous  energy  dissipation  effects  within  the  fluid. 
Since  the  amount  of  heat  generated  in  this  way  is  usually  ex¬ 
ceedingly  small  the  fluid  can  be  assumed  to  have  uniform 
temperature  and,  therefore,  constant  physical  properties. 


With  the  additional  assumption  that  the  flow  is  laminar 


then  the  continuity  and  momentum  equations  are  also  greatly 
simplified.  With  an  x,  y,  z  Cartesian  co-ordinate  system 
positioned  so  that  the  direction  of  flow  coincides  with  the 
z-axis  then,  if  the  velocity  in  the  z-direction  is  denoted 
by  w,  the  pressure  by  p  and  the  viscosity  of  the  fluid  by  jj, , 
the  continuity  equation  reduces  to 


(6) 


which,  when  substituted  in  the  momentum  equation,  together 
with  the  conditions  stated  above  gives 


x-direction: 


(7) 


y-direction:  =  0 

dy 


(3) 


b2w  gc  &P 

'  i  — 


bx2  &y2  fjL  bz 


z-direction: 


(9) 


I 
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where  gc  is  the  dimension  conversion  factor.  Equations  (7) 
and  (8)  simply  indicate  that  the  pressure  is  constant  in  the 


x-y  plane  for  any  value  of  z. 


The  velocity  distribution  and  flow  rate  for  laminar 
flow  in  a  rectangular  conduit  were  obtained  theoretically  by 
Cornish  .  He  integrated  equation  (9)  with  the  boundary 
condition  that  the  fluid  at  the  conduit  surface  had  zero 
velocity.  Cornish  considered  a  conduit  of  width  2a  and  depth 
2b  and  located  the  origin  of  his  co-ordinates  at  the  axis  of 
the  conduit.  With  the  x-axis  horizontal,  the  y-axis  vertical 


and  flow  in  the  z-direction  Cornish  obtained 

n-1 


w 


16  b2  gc  bp  \ 
~  3  ,  /  Tw 


V3  fJL  ^  n=  1 , 3 .  . 


(-D 


n* 


mrx 

cosh  2b  titty 
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for  the  velocity  distribution,  and 


4  a  b3  bp 
3/i  6z 


,  192  b 
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n=l , 3. . 


nJ 


tanh 


mra 

2b 


(11) 


for  the  volumetric  flow  rate. 

It  is  interesting  to  note  that  as  i.e.  the  conduit 

becomes  infinitely  wide,  both  equations  (10)  and  (11)  reduce 
to  the  corresponding  expressions  for  flow  between  infinitely 


wide  parallel  plates. 
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Cornish  performed  experiments  with  water  in  a  rec¬ 


tangular  conduit  with  an  aspect  ratio,  H  of  2.92  and  verified 
equation  (11) .  Equation  (11)  may  be  rearranged  to 


f  =  i 
2 


De 


1  - 


192  b 

t r5  a 


I 


n=  1 , 3 .  .  _ 


-i~  tanh 
n3 


mra 

2b 


(12) 


where  f  and  Re ' are  the  Fanning  friction  factor  and  Reynolds 
number  respectively  and  are  defined  by 


f  = 


::  _  -A  HE 


1  gc  d( 


2  dz  P  V2 


(13) 


and  Re’  = 


(14) 


in  which  V  is  the  average  fluid  velocity  computed  by  dividing 
the  volumetric  flow  rate  by  the  cross-sectional  area  of  the 
conduit,  P  is  the  fluid  density  and  De  is  the  equivalent 
diameter  of  the  conduit  defined  by  four  times  the  cross-sectional 
area  divided  by  the  wetted  perimeter. 


Hence 


_  4_a_b 
e  a  +  b 


(15) 


The  relationship  between  the  friction  factor  and  the 
Reynolds  number  for  laminar  flow  in  a  rectangular  conduit  is 
therefore  not  unique  but  depends  on  the  aspect  ratio.  As 
— >oo  equation  (12)  reduces  to 


f  = 


24 

Re ' 


(16) 


43 


which  is  the  relationship  for  infinitely  wide  parallel  plates. 


No  extensive  results  for  the  turbulent  velocity  distri- 

(53) 


bution  in  rectangular  conduits  have  been  reported 


However, 


Nikuradse 


(72) 


did  measure  point  velocities  in  non-circular 


ducts,  including  rectangular  conduits,  and  found  that  secondary 
motions  of  the  fluid  particles  into  the  corners  and  away  from 
the  walls  were  superimposed  upon  the  longitudinal  motion . 

Claibourne in  a  review  of  pressure  gradient  data 
obtained  by  many  investigators  for  turbulent  flow  in  smooth 
non- circular  conduits  concluded  that  the  results  could  be 
correlated  by  means  of  the  Blasius  equation 


f  =  0.0792  (Re')  ^ 


(17) 


The  conduits  had  annular,  triangular  and  rectangular  cross 
sections  and  the  scatter  of  the  data  was  t  15  percent  for 
Re  >  5000  and  -  35  percent  for  2100  < Re  < 5000 „ 


B.  LAMINAR  STRATIFIED  FLOW  OF  TWO  IMMISCIBLE  LIQUIDS  IN  A 
RECTANGULAR  CONDUIT 


1.  Expressions  for  velocity  distribution  and  flow  rates 

The  development  of  the  velocity  distribution  and  flow 
rate  expressions  for  the  co- current  laminar  stratified  flow 
of  two  immiscible,  Newtonian  incompressible  fluids  in  a 


' 
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rectangular  conduit  is  now  undertaken.  The  flow  is  assumed 
to  be  fully  developed  and  the  interface  between  the  fluid 
layers  is  assumed  to  be  smooth  and  horizontal.  The  flow 
model  is  depicted  in  Figure  3.  The  conduit  is  represented 
by  ABCD  and  the  interface  by  EF.  Liquid  2  flows  above  liquid 
1,  the  depths  being  b2  and  b]_  respectively.  The  width  of 
the  conduit  is  2a.  The  x,  y,  z  co-ordinate  system  has  the 
origin  at  0  with  the  x-axis  co- incident  with  the  interface 
and  flow  in  the  z- direction. 

The  basic  differential  equation  governing  the  flow  of 
both  fluids  is  equation  (9)  which  may  be  written  for  each 
phase  as 


62wj  d2wj  _  gc  dp 

dx2  dy2  fj^  dz 


(18) 


w  is  the  point  velocity  which  is  a  function  of  both  x  and  y. 
The  boundary  conditions  are: 


w-l 

=  0 

when 

x  = 

-  a,  -b]_  <  y  <  0 

(19) 

W1 

=  0 

when 

-a^x<a,  y  =  -bj_ 

(20) 

w2 

=  0 

when 

x  = 

±  a ,  0^y^b2 

(21)  * 

w2 

=  0 

when 

-a<x^a,  y  =  b2 

(22) 

w^ 

CM 

£ 

11 

when  -a<x<a, 

(23) 

Mi 

dwi 

y  =  0 

(24) 
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FIGURE  3  -  MODEL  OF  STRATIFIED  FLOW. 
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Let 


and 


(25) 

(26) 


If  the  variable  Wj  is  defined  by  the  relationship 


w .  =  w  •  + 


±1 


(x  -  az) 


(27) 


then  substitution  for  Wj  in  the  differential  equation  (18) 


gives 


d2W-}  d2W-j 
- 1  +  - 1  =  0 

V 


br2  >-2 


(28) 


The  boundary  conditions  (19)  through  (24)  now  become 


W1  = 

0  when  x  =  1  a,  -b^<y<0 

(29) 

II 

1 — 1 
& 

_  hi  (x2  -  a2)  when 

■a^x^a,  y  =  -b^ 

(30) 

W  2  = 

0  when  x  =  -a  0^y^b2 

(31) 

W2  = 

-  —  (x2  -  a2)  when  - 
2 

■a^x^a,  y  =  b2 

(32) 

w2  - 

W1  =  -|(ki  "  k2^  (x2"a2) 

>  when  -a^x^a. 

{33) 

dwi 

-  m  &W2  -  0 

•C 

II 

o 

(34) 

by 

by  J 

Boundary  conditions  (29)  and  (31)  give  Wj  =  0  when  x  =  la 
and  are  staisfied  by  terms  of  the  form 


cos 


(2i  +  1)  7T 

2a 


(35) 
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where  Yj  is  a  function  of  y  only  and  i  is  an  integer, 


Let 


n  = 


(2i  +  1)  7T 

2a 


(36) 


and  substitute 


W. 

3 


cos  nx 


in  equation 
i.  e. 


b2  Yj 

"V” 


2  „ 
n  Y. 


0 


(37) 


(38) 


from  which  it  follows  that 

Yj  =  Aj  (n)  sinh  ny  +  Bj  (n)  cosh  ny  (39) 

in  which  Aj (n)  and  Bj (n)  are  functions  of  n  and,  therefore, 
of  i  but  not  of  x  and  y.  The  solution  may  therefore  be 

written  in  the  form 

oo 

Wj  =  I[  Aj  (n)  sinh  ny  +  Bj  (n)  cosh  ny  j  cos  nx  (40) 

i=0 


which  may  be  considered  as  a  cosine  Fourier  series,  the  term 
coefficients  being  functions  of  y.  Two  of  the  six  boundary 
conditions  have  been  used.  Four  remain,  therefore,  for  the 
evaluation  of  Aj (n)  and  Bj (n) ,  j  =  1,2. 

Substitution  of  boundary  condition  (30)  into  equation 
(40)  gives 


oo 


-Zr  A-^  (n)  sinh  (-nb-jJ  +  B-^  (n)  cosh  (-nb] )  j  cos  nx 


(41) 
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Multiplication  of  both  sides  of  equation  (41)  by  cos  n'x  where 


n'  is  any  one  of  the  possible  values  of  n  and  integration 

with  respect  x  between  the  limits  -a  and  +a  gives 

+a 


_  kl ,  2 


+a 


(x  -a  ) cos  n'x  dx  = 


©o 

y 

(n)  sinh  (-nb-^)  + 


-a 


+  B1 (n) cosh (-nbj) J  cos  nx  cos  n'x  dx 


+a 


Since 


J cos  nx  cos  n'x  =  0  for  n  ^  n1 
-a 


equation  (42)  becomes 
+a 


(x2-a2) cos  nx  dx  = 


+a 


(n)  sinh  (-nb^)  +  (n)  cosh  nb-j_  J ~os“nx  dx 

-a  (44) 


(42) 

(43) 


from  which 


-  A-^(n)sinh  nbj_  +  B^  (n)  cosh  nb-^ 


2ki(-l) 
.3 


a  n 


(45) 


By  substituting  boundary  condition  (32)  into  equation  (40) 
and  following  a  procedure  similar  to  that  used  for  deriving 


equation  (45)  it  follows  that 


2k2  (“1) 

A2 (n) sinh  nb2  +  B2 (n) cosh  nb2  =  —  ~3~~ 


(46) 


Substitution  of  boundary  condition  (33)  in  equation  (40)  gives 

+a  +a 

^B2  (n)  -  Bj_(n)j  I  cos2nx  dx  =  -i  j  (k1~k2)  (x2-a2)cos  nx  dx  (47) 

-a  2  “'-a 
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from  which 


_  ,  ,  „  ,  ,  _  2(ki-ka)  (-1) 

a  n 


i+1 


(48) 


By  differentiating  equation  (40)  with  respect  to  y  and  using 
boundary  condition  (34)  it  follows  that 


A 


(n)  =  m  A2  (n) 


(49) 


Equations  (45) ,  (46) ,  (48)  and  (49)  may  be  solved 

simultaneously  for  the  coefficients  A^  (n)  ,  A2  (n)  ,  (n)  ,  Bp  (n) 

Thus  i 

^Hk2CQSh 


Ax  (n)  = 


A2  (n)  = 


2izlL 

a  n*5 


nbi-kicosh  nb2+  (ki~k2)  cosh  nb]  cosh  nb2j 

cosh  nb^  sinh  nb2  +  m  sinh  nb-j  cosh  nb2 

(50) 

k2COsh  nb^-k^cosh  nb2+  (k^-k2)  cosh  nb^  cosh  nb2j 


Bi(n)  = 


2(.Tl). 

Irr 


cosh  nbj_  sinh  nb2  +  m  sinh  nb^  cosh  nb2 

(51) 

k]sinh  nb2  +  mk2 s inh  nbi  +  m  (kq-k2)  sinh  nbi  cosh  nb2] 


cosh  nb-]  sinh  nb2  +  m  sinh  nb^  cosh  nb2 


B2  (n)  = 


(52) 


2i-i)„1r  v  '  1 

a  n3  Lmk2sinh  nbi+kisinh  nb2~ (ki~k2) sinh  nb2  cosh  nbij 


cosh  nb-]Sinh  nb2  +  m  sinh  nb^  cosh  nb2 


(53) 


The  desired  solution  for  the  velocity  distribution  is 


oo 


w- 


=  ^  [a-,  (n)  sinh  ny  +  B^  (n)  cosh  ny  1  cos  nx  +  — (x^-a2)  (54) 

i«0  J  2 
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with  A-j  (n)  and  (n)  given  by  equations  (50)  and  (52)  res¬ 


pectively  and 

oo 


V7r 


=  [a2  (n) sinh  ny  +  B2 (n) cosh  ny j  cos  nx  +  —  (x2-a2)  (55) 


i=0 


with  A2 (n)  and  B2 (n)  given  by  equations  (51)  and  (53). 

Immediate  checks  on  the  velocity  distribution  expressions 
were  obtained  in  two  ways: 

(i)  If  m  =  1  and,  therefore,  k-^  =  k2  then  fluids  1 
and  2  are  identical  as  far  as  this  problem  is  concerned.  If, 
in  addition,  b^  =  b2  =  b  then  the  x-axis  lies  equidistant 
between  the  top  and  bottom  of  the  conduit.  With  these 
substitutions  both  equations  (54)  and  (55)  reduce  to  the 
expression  derived  by  Cornish  for  the  velocity  distribution 
in  single- phase  flow,  equation  (10) . 

(ii)  If  y  =  0  and  b^  =  b2  =  b  then  both  equations  (54) 
and  (55)  reduce  to  the  equation 


2  9  ©o 

w  =  (x  -a  )  (ki+mk2)  + 


2  (1  +  m) 


Z 


i=0 


2  (-1)  (ki+mk2) 
a  n^  (1+m)  cosh  nb 


cos  nx 


(56) 


(24) 

which  has  been  obtained  by  Dumitrescu  and  Staniscu '  '  for 
the  interfacial  velocity  distribution  for  the  special  case 
when  the  depths  of  the  fluid  layers  are  equal. 
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Since  it  was  desired  to  evaluate  the  velocity  distribution 

by  automatic  digital  computation,  equations  (54)  and  (55)  were 

rearranged  into  more  convenient  forms  as  follows 
oo 

Wl  -  ad  l  -tiiL 
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and 


8a 
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oo 
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2 


(58) 


in  which 


A1  (n) 


(n) 


A2(n)  " 

k2  cosh  hfcq  -  cosh  nb2  +  (kl“k2^  cosk  nkl  cosh  nb2 
cosh  nb^  sinh  nb2  +  m  sinh  nb^  cosh  nb2 

(59) 

k]_  sinh  nb2  +  mk2  sinh  nb^  +  m  (kp-k2)  sinh  nb]_  cosh  nb2 
cosh  nbj_  sinh  nb2  +  m  sinh  nb^  cosh  nb2 

(60) 


B2(n) 


mk2  sinh  nb^  +  kj_  sinh  nb2  -  (k^-k2)  sinh  nb2  cosh  nbj 


cosh  nb^  sinh  nb2  +  m  sinh  nbj_  cosh  nb2 


(61) 


The  volumetric  flow  rates  and  Q2  are  obtained  by  evaluating 
the  integrals  o  +a 


(62) 
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bo  +  a 


and  Q2  =  I  I  w2  dx  dy  (63) 

0  -a 

Substitution  for  and  w2  from  equations  (54)  and  (55)  gives 
finally 


oo 


Q 


i=  I  ^ 


i-0  n" 


Ajl  (n)  (1  -  cosh  nb^)  +  B^  (n)  sinh  nb-^ 


2  ,  3  , 

3  kl  a  bl 


(64) 


and 


oo 


o  =  y  2  (-D 1 

2  1=0  "2 


A2  (n)  (cosh  nb2  -  1)  +  B2  (n)  sinh 


2  k  a3 

3  k2  a  b2 


nb 


(65) 


with  A^  (n)  ,  A2  (n)  ,  B-^  (n)  and  B2  (n)  given  by  equations  (50) 
through  (53) . 

Alternatively  equations  (64)  and  (65)  may  be  written 


in  the  form 


oo 


A 


Qi  = 


_  128a~*  Y  X 


it5  to  (2i+D5 


mA^  (n)  (1-cosh  nb-jJ +B  j  (n)  sinh  nb^ 


“  i  kl  *  bi 


(66) 


and 


oo 


Qo  = 


123a 


2  -n-5 


tt5  i-0  (2i+l)  5 


A2 (n)  (cosh  nb2~l) +B2  (n) sinh  nb2 


-  T  k2  a3  b2 


(67) 


with  Aj^  (n)  ,  A2  (n)  ,  B^  (n)  and  B^  (n)  given  by  equations  (59) 


through  (61) . 
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The  following  check  was  obtained  for  the  flow  rate 
expressions.  If  b-^  =  I02  =  b  and  m  =  1  then  fluids  1  and  2  are 
identical  and  Q-^  =  Q2*  The  sum  of  and  C>2  should  be  identical 
with  the  flow  rate  expression  derived  by  Cornish,  equation  (11) 
for  single  phase  flow.  This  is,  in  fact,  so. 

Expressions  have  been  derived  in  open  form  for  the 
velocity  distribution  and  volumetric  flow  rates  for  the  strati¬ 
fied,  laminar  flow  of  two  fluids  in  a  closed  rectangular  conduit. 
The  velocity  distribution  and  flow  rates  can,  therefore,  be 
calculated  from  a  knowledge  of  the  physical  dimensions  of  the 
conduit,  the  fluid  properties,  the  pressure  gradient  and  the 
position  of  the  interface. 

It  is  now  convenient  to  define  several  terms  which  will 
be  useful  in  the  presentation  of  the  theoretical  results. 

Reference  was  made  in  the  introduction  to  the  pressure  gradient 
reduction  factor  for  the  simultaneous  flow  of  two  fluids  through 
closed  conduits.  The  pressure  gradient  reduction  factor  is 
defined  as  the  ratio  of  the  pressure  gradient  for  the  flow  of 
the  more  viscous  phase  alone  to  the  pressure  gradient  for  the 
flow  of  the  two-phase  system,  the  flow  rate  of  the  more  viscous 
phase  being  the  same  in  each  case.  In  presenting  the  results 
of  the  present  investigation  it  will  be  assumed  that  the  viscosity 
of  the  upper  liquid  2  is  equal  to  or  greater  than  the  viscosity 


59 


of  the  lower  liquid  1  i.e.  m>l. 


Hence 


pressure  gradient  reduction  factor  = 


>z  /  2 


dp] 

dzj 


for  same  Q2  (68) 


tp 


Now,  since  for  a  given  interface  position  flow  rate  and  pressure 
gradient  are  directly  proportional  to  each  other  the  pressure 


gradient  reduction  factor  may  also  be  written  in  the  form 

(Q2) tP 


pressure  gradient  reduction  factor  = 


for  same 


(69) 


Q- 


equation  (67) _ 

“  equation  (11)  written 
for  liquid  2 

for  same  (70) 

By  dividing  equation  (67)  by  equation  (11)  it  is  possible 
to  show  that  the  pressure  gradient  reduction  factor  is  a 
function  of  the  interface  position,  the  viscosity  ratio  and 
the  conduit  aspect  ratio  only.  Hence  for  different  aspect 
ratios  it  is  possible  to  represent  the  pressure  gradient  reduction 
factor  as  a  function  of  the  interface  position  and  viscosity 
ratio. 

Hie  position  of  the  interface  will  be  defined  by 

b2 

relative  interface  position  =  - 

bl  +  b2 


(71) 
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The  relative  interface  position  varies  from  0  when  liquid  1 
flows  alone  to  1.0  when  liquid  2  flows  alone. 

In  general,  a  characteristic  of  two-phase  flow  is  that 
the  average  velocities  of  the  two  phases  in  the  conduit  are 
not  equal.  This  causes  the  insitu  ratio  of  the  phases  to  be 
different  from  the  input  ratio  or  flow  rate  ratio. 


If 


and 


Ql 

input  ratio  =  — 

q2 

.  .  .  bl 

insitu  ratio  =  - — 

b2 


(72) 

(73) 


then  the  holdup  ratio  is  defined  by 

,  _  _  ,  .  input  ratio 

holdup  ratio  =  ; - ; —  (74) 

insitu  ratio 

2 .  Computation 

Representative  velocity  profiles  were  evaluated  using 
equations  (57)  and  (58)  to  illustrate  the  effects  of  interface 
position  and  aspect  ratio  on  the  velocity  distribution.  Pressure 
gradient  reduction  factors  and  holdup  ratios  were  evaluated 
from  equations  (11) ,  (66)  and  (67)  for  the  entire  range  of 

interface  position  for  viscosity  ratios  from  1  to  1000  and 
aspect  ratios  from  8:1  to  1:3.  Equation  (11)  was  evaluated 
with  the  aid  of  a  desk  calculator  and  the  other  equations  with 
the  aid  of  the  University  of  Alberta  IBM  1620  computer  having 
punched  card  input  and  output.  Programs  were  written  in 
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Fortran  II. 

All  of  the  equations  evaluated  consist  of  two  parts  - 
one  in  open  form  and  the  other  in  closed  form.  In  each  case 
the  sum  of  the  open  form  expression  was  obtained  until  the 
value  of  the  next  term  was  less  than  0.01  percent  of  the  sum 
already  accummulated.  That  this  was  a  reliable  procedure 
leading  to  negligible  error  was  confirmed  in  two  ways: 

(i)  Interfacial  velocities  evaluated  from  equations  (57) 
and  (58)  always  agreed  to  within  approximately  1  percent. 

Usually  between  four  and  ten  terms  were  required  in  the  infinite 
series  to  give  the  stated  tolerance.  The  series  in  the  flow 
rate  expressions  converged  faster  than  the  series  in  the 
velocity  expressions  showing  that  the  error  introduced  by 
neglecting  terms  in  the  flow  rate  expressions  having  a  magni¬ 
tude  less  than  0.01  percent  of  the  sum  would  cause  even  less 
error  than  in  the  velocity  expressions. 

(ii)  The  sum  of  the  flow  rates  calculated  from  equations 
(66)  and  (67)  for  any  interface  position  and  a  viscosity 
ratio  of  unity  agreed  to  within  approximately  0.5  percent  with 
the  flow  rate  for  single-phase  flow  evaluated  from  equation  (11) . 

It  is  worthwhile  to  consider  the  convergence  of  the 
series  in  more  detail.  Consider,  for  example,  equation  (67) 


I 

t 

* 
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which  is  the  flow  rate  expression  for  liquid  2.  The  series 
converges  because  as  n  increases  the  two  terms  A'^Cn) (cosh 
nb2  -  1)  and  (n) sinh  nb2  become  more  nearly  equal  in  magni¬ 
tude  and  opposite  in  sign.  Further ,  for  a  given  value  of  n, 
the  values  of  cosh  nb2  and  sinh  nb2  become  very  large  for 
large  values  of  b2/a,  i.e.  for  small  aspect  ratios.  Consequently 
it  would  be  expected  that  machine  round-off  error  would  be 
more  troublesome  for  computations  involving  small  aspect 
ratios.  This  was,  in  fact,  found  to  be  so.  Whereas  ordinary 
precision  floating  point  arithmetic  in  which  8  figures  were 
carried  was  sufficient  for  an  aspect  ratio  of  8:1  extended 
precision  floating  point  arithmetic  in  which  28  figures  were 
carried  was  necessary  for  an  aspect  ratio  of  1:3.  Round-off 
error  manifested  itself  in  extremely  "wild"  results  being 
obtained  which  disappeared  when  sufficient  figures  were  carried. 

3.  Results 

(i)  Velocity  distributions 

The  evaluation  of  complete  velocity  distributions  for 
wide  ranges  of  interface  position,  viscosity  ratio  and  aspect 
ratio  was  not  warranted  in  the  present  study.  Rather,  repre¬ 
sentative  profiles  were  calculated  to  illustrate 

(1)  the  effect  on  the  profile  of  different 
interface  positions  for  given  viscosity  and  aspect  ratios  and 
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(2)  the  effect  on  the  profile  of  different 
aspect  ratios  for  a  given  interface  position  and  viscosity- 
ratio. 

The  effect  of  interface  position  is  illustrated  in 
Figure  4  in  which  velocity  profiles  on  the  y-axis  are  shown 
for  a  viscosity  ratio  of  5.326  and  an  aspect  ratio  of  7.95:1. 
The  velocities  are  represented  by  the  dimensionless  ratio 
w/Vpp  where  w  is  the  local  velocity  calculated  for  a  given 
pressure  gradient  and  Vpp  is  the  average  velocity  for  flow  of 
the  more  viscous  phase  under  the  same  pressure  gradient  between 
parallel  plates  (b^  +  b 2 )  apart.  In  this  way  the  results  are 
independent  of  the  pressure  gradient.  The  relative  interface 
position  is  shown  for  each  profile.  A  relative  position  of 
0  corresponds  to  the  flow  of  the  less  viscous  phase  alone  and 
a  relative  position  of  1  corresponds  to  the  flow  of  the  more 
viscous  phase  alone. 

It  is  apparent  that  the  addition  of  the  less  viscous 
phase  extends  the  profile  within  the  more  viscous  phase.  In 
doing  so  the  maximum  in  the  profile  always  occurs  at  a  point 
below  the  mid- point  between  the  top  and  bottom  of  the  conduit. 
For  relative  interface  positions  less  than  0.7  approximately 
the  maximum  occurs  in  the  less  viscous  phase;  for  relative 
positions  greater  than  approximately  0.7  the  maximum  occurs 
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FIGURE  4  -  VERTICAL  VELOCITY  PROFILES 
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in  the  more  viscous  phase.  At  a  relative  position  of  approxi¬ 
mately  0.7  the  maximum  velocity  occurs  at  the  interface.  It 
also  appears  from  the  figure,  although  it  cannot  be  considered 
as  proven,  that  for  a  given  pressure  gradient  and  this  parti¬ 
cular  aspect  ratio  the  interface  velocity  is  greatest  when 
the  maximum  in  the  complete  profile  occurs  at  the  interface. 

The  effect  of  aspect  ratio  is  illustrated  in  Figure  5 
in  which  interfacial  velocity  profiles,  i.e.  at  y  =  0,  are 
shown  for  a  relative  interface  position  of  0.9  and  a  viscosity 
ratio  of  5.326.  The  velocities  are  represented  by  the  same 
dimensionless  ratio  used  in  Figure  4  and  are,  therefore, 
independent  of  pressure  gradient.  Aspect  ratios  range  from 
oo:l,  i.e.  infinitely  wide  parallel  plates,  to  1:3.  The 
interface  velocity  for  an  aspect  ratio  ofoo:l  was  computed 
from  the  equation 

w  =  _  Is  fol  +  b2)bl  b2  (7 

1  2  bz  (b2  tlx  +  bx'  fl2) 

which  may  be  obtained  by  extending  the  analysis  of  Russell  and 
Charles  (78) . 

It  is  apparent  that  for  an  aspect  ratio  of  7.95:1  the 
velocity  profile  for  x/a  less  than  0.5  is  coincident  for  all 
practical  purposes  with  the  profile  for  flow  between  parallel 
plates  and  is,  therefore,  a  function  of  y  only.  For  such  an 


' 


66 


a 


FIGURE  5  -  INTERFACIAL  VELOCITY  PROFILES 
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aspect  ratio  the  effect  of  the  sides  of  the  conduit  are  only 
apparent  for  x/a  greater  than  0.5.  For  an  aspect  ratio  of  3:1 
the  effect  of  the  sides  is  transmitted  across  the  width  of  the 
conduit  so  that  no  part  of  the  profile  is  two-dimensional  and, 
in  addition,  the  value  of  the  velocity  at  x/a  =  0  is  approxi¬ 
mately  5  percent  less  than  the  parallel  plate  value.  The 
profiles  for  aspect  ratios  of  1:1  and  1:3  show  that,  as 
would  be  expected,  the  effect  of  the  sides  of  the  conduit 
increases  as  the  aspect  ratio  decreases. 

The  numerical  data  used  to  prepare  Figures  4  and  5  are 
tabulated  in  Table  A1  of  Appendix  A. 

(ii)  Pressure  gradient  reduction  factors 
Pressure  gradient  reduction  factors  were  evaluated  over 
wide  ranges  of  interface  position  and  viscosity  ratio  for 
aspect  ratios  of  7.896:1,  3:1,  1:1  and  1:3.  Viscosity  ratios 
ranged  from  1  to  1000.  Pressure  gradient  reduction  factors 
are  presented  in  Figure  6  for  an  aspect  ratio  of  7.896:1. 

The  curves  for  each  viscosity  ratio  have  the  same  general 
shape.  For  a  relative  interface  position  of  1.0  the  interface 
is  at  the  bottom  of  the  conduit,  the  more  viscous  phase  flows 
alone  and  the  pressure  gradient  reduction  factor  is  unity. 

With  increasing  amounts  of  the  less  viscous  phase  the  pressure 
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gradient  reduction  factor  increases  to  a  maximum  and  then 
decreases  to  0  at  a  relative  interface  position  of  0.  Under 
conditions  when  the  pressure  gradient  reduction  factor  exceeds 
unity  the  two- phase  pressure  gradient  is  less  than  the  single¬ 
phase  pressure  gradient,  despite  the  fact  that  the  more 
viscous  phase  flows  at  the  same  rate  in  both  cases.  The 
factor  by  which  the  pressure  gradient  can  be  reduced  increases 
with  viscosity  ratio  but  becomes  less  sensitive  to  variation 
in  viscosity  ratio  at  high  viscosity  ratios. 

Similar  plots  showing  the  pressure  gradient  reduction 
factor  as  a  function  of  interface  position  and  viscosity  ratio 
are  shown  in  Figures  7,  8  and  9  for  aspect  ratios  of  3:1,  1:1 
and  1:3  respectively.  It  is  apparent  that  as  the  aspect  ratio 
decreases  the  addition  of  a  less  viscous  phase  to  decrease 
the  pressure  gradient  becomes  increasingly  unattractive.  Thus 
for  a  viscosity  ratio  of  1000,  while  the  maximum  reduction  factor 
is  3.38  for  an  aspect  ratio  of  7.896:1,  it  is  only  about  1.13 
for  an  aspect  ratio  of  1:3. 

Superimposed  on  the  curves  for  a  square  conduit  in  Figure 
8  are  curves  for  a  circular  pipe  obtained  by  Charles  and  Red- 

(13) 

berger v  '  .  It  is  immediately  apparent  that  for  a  given 
viscosity  ratio  the  maximum  reduction  factor  is  less  for  a 
circular  pipe  than  it  is  for  a  square  conduit.  In  addition. 


C- 


-  70  - 


0  0.2  0.4  0.6  0.8  1.0 

RELATIVE  INTERFACE  POSITION 


PRESSURE  GRADIENT  REDUCTION  FACTOR 


71 


FIGURE  8  -  PRESSURE  GRADIENT  REDUCTION  FACTORS  FOR  AN  ASPECT 

RATIO  OF  1:1 
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FIGURE  9  -  PRESSURE  GRADIENT  REDUCTION  FACTORS  FOR  AN  ASPECT 

RATIO  OF  1:3 
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the  depth  of  the  less  viscous  phase  required  to  produce  the 
maximum  value  of  the  reduction  factor  is  considerably  greater 
for  the  pipe  than  for  the  square  conduit.  As  a  consequence 
the  proportion  of  the  less  viscous  phase  in  the  flowing  stream 
for  a  maximum  pressure  gradient  reduction  is  considerably 
greater  for  flow  in  a  pipe. 

Figure  10  compares  the  maximum  reduction  factors  for 
the  different  aspect  ratios  and  a  circular  pipe.  Included  in 
this  figure  is  the  curve  evaluated  from  the  expressions 
developed  by  Russell  and  Charles  ^8)  for  stratified  flow 
between  infinitely  wide  parallel  plates,  the  curve  for  a 
circular  pipe  evaluated  by  Charles  and  Redberger ^3)  and  the 
curves  for  aspect  ratios  of  7.396:1,  3:1,  1:1  and  1:3  for 
which  the  maximum  pressure  gradient  reduction  factors  were 
read  from  the  curves  shown  in  Figures  6,  7,  8  and  9.  The 
steady  decrease  in  the  maximum  reduction  ractor  with  decreasing 
aspect  ratio  is  apparent.  For  a  given  aspect  ratio  the 
maximum  reduction  factor  tends  to  a  constant  value  at  high 
viscosity  ratios. 

It  is  worthwhile  to  consider  the  cases  of  a  square 
conduit  and  a  circular  pipe  further.  At  a  viscosity  ratio  of 
1000  not  only  is  the  maximum  reduction  factor  1.595  for  a 
square  conduit  and  1.31  for  a  circular  pipe  but  the  greater 
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reduction  is  achieved  with  a  considerably  less  proportion  of 
the  less  viscous  phase  -  only  1.5  percent  of  the  total  flow 
in  a  square  conduit  compared  with  90  percent  in  a  circular 
pipe.  This  is  a  consequence  of  the  relatively  higher  inter¬ 
face  position  in  pipe  flow  allowing  the  less  viscous  phase 
a  relatively  easy  passage  through  the  pipe. 

The  curves  shown  in  Figure  10  for  the  different  aspect 
ratios  are  cross  plotted  in  Figure  11  showing  in  more  detail 
how  the  maximum  reduction  factor  varies  with  aspect  ratio. 

Also  included  are  the  maximum  reduction  factors  for  an  infinite 
aspect  ratio. 

(iii)  Holdup  ratios 

The  holdup  ratio  has  been  defined  as  the  input  ratio 
(less  viscous  phase  to  more  viscous  phase)  divided  by  the 
insitu  ratio.  As  a  consequence,  a  holdup  ratio  greater  than 
unity  indicates  that  the  actual  average  velocity  of  the  less 
viscous  phase  in  the  conduit  is  greater  than  the  actual 
average  velocity  of  the  more  viscous  phase. 

Curves  of  holdup  ratio  plotted  against  insitu  ratio 
with  viscosity  ratio  as  parameter  are  presented  for  an  aspect 
ratio  of  7.896:1  in  Figure  12.  As  would  be  expected  the 
holdup  ratio  is  unity  when  both  the  insitu  ratio  and  viscosity 
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ratio  are  unity.  It  is  apparent  that  the  higher  the  viscosity 
ratio,  the  greater  is  the  tendency  for  the  average  velocity 
of  the  less  viscous  phase  to  exceed  that  of  the  more  viscous 
phase.  At  the  higher  viscosity  ratios  it  is  only  at  very  low 
insitu  ratios  that  the  average  velocity  of  the  more  viscous 
phase  exceeds  that  of  the  less  viscous  phase.  The  manner  in 
which  the  holdup  ratio  varies  with  insitu  ratio  is  readily 
appreciated  when  reference  is  made  to  the  velocity  profiles 
presented  in  Figure  4. 

Similar  plots  of  holdup  ratio  against  insitu  ratio  for 
aspect  ratios  of  3:1,  1:1  and  1:3  are  presented  in  Figures  13, 
14  and  15.  While  all  the  families  of  curves  are  basically 
similar  a  definite  trend  is  evident  as  the  aspect  ratio 
decreases:  for  a  given  viscosity  ratio,  at  relatively  low 

insitu  ratios  the  holdup  ratio  increases,  while  at  relatively 
high  insitu  ratios  the  holdup  ratio  decreases.  The  following 
table  illustrates  this  trend  for  a  viscosity  ratio  of  100. 
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TABLE  I 

Variation  of  Holdup  Ratio  with  Aspect  Ratio 
for  a  Viscosity  Ratio  of  100 


Aspect  Ratio 
7.896:1 


Insitu  Ratio 

0.30 

10.0 


Holdup  Ratio 

3.20 

300 


1:3 


0.30 

10.0 


29.0 

115 


The  trend  reflects  the  increasing  effect  of  the  sides 
of  the  conduit  as  the  aspect  ratio  decreases.  At  low  insitu 
ratios  the  sides  tend  to  hold  back  the  more  viscous  phase, 
while  at  high  insitu  ratios  they  tend  to  hold  back  the  less 
viscous  phase. 

All  the  numerical  pressure  gradient  reduction  factor 
and  holdup  ratio  data  used  to  prepare  the  foregoing  figures 
are  tabulated  in  Table  A2  of  Appendix  A. 

4.  Analogues  of  stratified  flow 

Other  physical  processes  are  governed  by  the  same  form 
of  differential  equation  as  governs  steady  fully  developed 
unidirectional,  laminar  flow  of  a  fluid.  Thus,  for  example, 
the  differential  equations  for  the  stress  distribution 
resulting  from  the  torsion  of  a  long  solid  bar  and  the  tempera- 
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ture  distribution  created  by  uniform  heat  generation  within  a 
long  bar  have  the  same  form  as  equation  (9) .  The  latter 
example  is  of  more  interest  in  chemical  engineering  and  will 
be  considered  in  some  detail. 

If  heat  is  generated  uniformly  in  a  long  solid  bar  of 
constant  cross-sectional  area  by,  for  instance,  ionizing 
radiation  at  the  rate  of  q  BTU/unit  volume/unit  time  then  the 
differential  equation  governing  the  steady  state  temperature 
distribution  is 


d2T  +  62T  =  JL 

dx2  by2  ^ 


(76) 


where  T  is  temperature  and  k  is  the  thermal  conductivity. 


Consequently  the  thermal  conductivity  and  rate  of  heat  genera¬ 
tion  are  equivalent  to  viscosity  and  pressure  gradient  respect¬ 


ively. 


Consider,  now,  a  long  composite  bar  of  rectangular 
cross  section  consisting  of  two  materials  of  different  thermal 
conductivity  joined  together  such  that  no  temperature  differ¬ 
ence  exists  across  the  interface  and  no  heat  is  generated  or 
removed  in  the  interface  itself.  If  the  outside  of  the  bar  is 
maintained  at  some  reference  temperature  then  the  solution  of 
the  differential  equation  (76)  is  identical  in  form  to  the 
solution  of  the  laminar,  stratified  flow  problem. 


■ 
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The  velocity  profiles  already  presented  are,  therefore, 
representative  of  the  corresponding  temperature  profiles  of 
the  analogue.  Further,  the  pressure  gradient  reduction  factor 
is  representative  of  the  ratio  of  the  amount  of  heat  stored 
in  the  material  having  the  higher  thermal  conductivity  in  the 
composite  bar  to  the  amount  which  would  be  stored  in  a  bar 
having  the  same  outside  dimensions  as  the  composite  bar  but 
composed  only  of  the  material  having  the  higher  thermal 
conductivity.  Finally,  the  holdup  ratio  is  representative  of 
a  ratio  which  would  indicate  which  part  of  the  bar  contained 
the  more  heat,  provided,  of  course,  the  specific  heats  of  the 


two  materials  were  the  same. 


Ill  EQUIPMENT 


Equipment  was  designed  and  constructed  so  that  two 
immiscible  liquids  could  be  pumped  simultaneously  through  a 
rectangular  conduit  of  large  aspect  ratio.  In  this  section 
the  selection  of  the  fluids  and  the  design  of  the  rectangular 
conduit  and  auxiliary  equipment  will  be  discussed.  The 
arrangement  of  the  experimental  equipment  is  illustrated 
diagrammatically  in  Figure  16. 

A.  FLUIDS 

The  liquids  were  selected  with  the  following  three 
factors  in  mind: 

(1)  the  liquids  should  have  a  sufficient  viscosity 
differential  for  the  laminar- turbulent  regime  to  cover  wide 
ranges  of  individual  flow  .rates, 

(2)  the  liquids  should  separate  readily  under  the 
action  of  gravity  and 

(3)  both  liquids  should  be  colourless  so  that  the 
interfacial  structure  would  not  be  obscured. 

Water  and  a  refined  mineral  oil  with  a  viscosity  of  approxi¬ 
mately  five  centipoise  were  chosen.  The  oil  was  a  5:1  blend 
of  Mentor  29  and  Bayol  35,  both  supplied  by  Imperial  Oil  Ltd. 
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FIGURE  16  -  EXPERIMENTAL  EQUIPMENT 
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The  physical  properties  of  the  liquids  at  the  average  tempera¬ 
ture  of  the  experiments  are  given  in  Table  II  and  an  account 
of  their  determination  and  variation  with  temperature  is 
given  in  Appendix  B. 


Table  II 


Physical  Properties 

of  Liquids 

Oil 

Water 

Specific  gravity 

0.818 

0.997 

Viscosity  (centipoise) 

4.82 

0.905 

Density  ratio  =  0.821 
Viscosity  ratio  =  5.326 
Interfacial  tension  =  42  dynes/cm. 

B.  RECTANGULAR  CONDUIT 

The  rectangular  conduit  had  an  overall  length,  including 
the  entry,  calming,  test  and  discharge  sections,  of  37  feet 
and  a  nominal  1-inch  by  8- inch  cross  section.  Figure  17  is  a 
photograph  of  the  conduit. 

The  top  and  bottom  plates  of  the  conduit  were  1/2- inch 
thick  transparent  Plexiglas  plastic  and  the  sides  were  3/4-inch 
thick.  Except  for  a  photographic  section  the  top,  bottom  and 
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FIGURE  17-  RECTANGULAR  CONDUIT 
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sides  of  the  channel  were  continuous  and  formed  by  joining 
6-foot  lengths  with  a  cement  consisting  of  plastic  chips 
dissolved  in  ethylene  dichloride.  The  faces  to  be  joined 
were  milled  to  a  45°  slope  to  give  the  joints  greater  strength. 
The  manner  in  which  the  top,  bottom  and  sides  of  the  conduit 
were  assembled  together  is  shown  in  Figure  18.  The  1/4- inch 
diameter  bolts  were  spaced  every  two  inches  along  the  entire 
length  of  the  conduit  and  electrical  insulation  tape  was  a 
satisfactory  gasket  material. 

The  conduit,  except  for  the  photographic  section, 
rested  freely  upon  two  2-inch  x  1-inch  channel  iron  runners 
which  were  bolted  to  supports  attached  rigidly  to  the  wall. 

The  conduit  was  thus  free  to  expand  and  contract  under  the 
influence  of  changes  in  ambient  temperature  since  all  connections 
to  the  conduit  contained  flexible  linkages.  The  conduit  was 
levelled  to  within  0.02  inches  over  its  entire  length  and  was 
completely  free  from  external  vibrations. 

The  sheet  metal  entry  section  was  designed  to  bring 
the  two  phases  into  contact  across  a  horizontal  interface  with 
as  little  turbulence  as  possible  and  was  similar  to  the  one 
used  by  Gazley'  '  to  introduce  air  and  water  to  a  circular 
pipe.  The  section  is  illustrated  in  Figure  19  and  was  found 
to  be  completely  satisfactory  over  the  entire  range  of  input 


Did .  bolt 


FIGURE  18  -  CROSS  SECTION  THROUGH  RECTANGULAR  CONDUIT 


Baffle  Flange  8  thick 
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FIGURE  19  -  ENTRY  SECTION 


92 


water-oil  ratios.  The  two  liquids  entered  the  section  at  the 
top  and  bottom  of  a  horizontal  cylinder  which  was  9  inches 
in  diameter  and  9  inches  long.  Interior  baffle  plates  prevented 
the  phases  from  impinging  directly  upon  each  other.  The 
circular  cross  section  gradually  changed  to  the  rectangular 
cross-section  over  a  distance  of  two  feet.  A  removable, 
circular  end  plate  held  in  position  by  twenty  four  1/4- inch 
bolts  enabled  the  section  to  be  cleaned  periodically.  Dial 
thermometers  were  inserted  through  openings  in  this  end  plate 
and,  therefore,  indicated  the  temperatures  of  both  phases  at 
the  inlet  to  the  rectangular  conduit.  The  entry  section  was 
joined  to  the  rectangular  conduit  by  a  flange  containing 
fourteen  1/4- inch  bolts  and  was  provided  with  a  vent  and  drain. 

The  calming  section  was  12.5  feet  or  84  equivalent 
diameters  long.  This  length  was  considered  to  be  sufficient 
for  the  velocity  profile  to  become  fully  developed  after  a 
study  of  the  experimental  results  obtained  by  Eckert  and 
Irvine who  measured  local  pressure  gradients  along  both 
a  rectangular  and  triangular  conduit.  The  rectangular  conduit 
had  an  aspect  ratio  of  3:1  and  pressure  gradients  were 
measured  with  both  abrupt  and  rounded  entrance  sections.  With 
the  abrupt  entrance  pressure  gradients  measured  for  the  flow 
of  air  alone  in  the  rectangular  conduit  at  Reynolds  numbers 


!ko  moztiod 
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ranging  from  300  to  11,600  indicated  that  the  calming  length 
amounted  to  no  more  than  40  equivalent  diameters.  With  the 
rounded  entrance  pressure  gradients  measured  for  Reynolds 
numbers  ranging  from  6,000  to  10,000  indicated  that  the  calming 
length  was  no  greater  than  70  equivalent  diameters.  Since  the 
entry  section  used  in  the  present  study  could  be  considered 
intermediate  between  the  abrupt  and  rounded  entrances  used 
by  Eckert  and  Irvine,  84  equivalent  diameters  were  judged  to 
be  entirely  sufficient  for  the  establishment  of  fully  developed 
flow. 

The  test  section  was  20  feet  or  135  equivalent  diameters 
long  and  the  discharge  section  was  2  feet  or  13.5  equivalent 
diameters  long.  Pressure  taps,  formed  by  drilling  1/16- inch 
diameter  holes  through  the  bottom  plate  along  the  centre  line 
of  the  conduit,  were  situated  every  five  feet  along  the  test 
section. 

The  18- inch  photographic  section  was  carefully  constructed 
so  that  its  internal  dimensions  matched  those  of  the  parts  of 
the  conduit  to  which  it  was  joined  to  within  0.001  inches. 

It  was  necessary  to  have  a  removable  photographic  section  so 
that  it  could  be  readily  cleaned.  To  ensure  dimensional 
stability  the  top  and  bottom  plates  of  the  photographic 
section  were  made  from  1-inch  thick  plastic.  Flanged  connections 


. 
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each  with  twelve  1/4- inch  diameter  bolts  and  locating  pins, 
joined  the  removable  section  to  the  upstream  and  downstream 
portions  of  the  conduit.  A  photograph  of  the  removable  section 
in  position  is  shown  in  the  upper  part  of  Figure  20. 

C.  AUXILLIARY  EQUIPMENT 

(79  12) 

Previous  experimental  studies  '  '  of  the  pipeline 

flow  of  oil-water  systems  involved  a  separation  procedure 
requiring  intermittent  operation.  In  the  present  study  more 
convenient  continuous  operation  was  achieved  by  using  a  combined 
storage  and  separator  tank.  The  two  phases  were  withdrawn 
individually  from  the  separator  by  positive  displacement  pumps 
and  rotameters  indicated  the  flow  of  the  two  phases  to  the 
test  section.  The  unwanted  portions  of  the  total  pump  dis¬ 
charges  were  either  returned  to  the  separator  tank  or  cir¬ 
culated  through  bypass  tanks.  Upon  leaving  the  rotameters 
the  liquids  passed  through  a  concentric  tube  heat  exchanger, 
which  promoted  thermal  equilibrium  between  them,  and  subsequently 
were  admitted  to  the  rectangular  conduit.  After  flowing  through 
the  conduit  the  oil  and  water  discharged  tangentially  into 
the  separator  where  they  settled  under  the  action  of  gravity 
to  form  two  distinct  layers.  Copper  pipe  was  used  throughout 
so  that  corrosion  would  be  kept  to  a  minimum.  Figure  16  shows 
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FIGURE  20  -  PHOTOGRAPHIC  SECTION  AND  CAMERAS 
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the  relative  positions  of  the  different  items  of  the  equipment. 
The  separator  consisted  of  a  52-inch  diameter,  75-inch 
high  cylindrical  tank  with  a  flat  bottom  and  removable  cover. 

The  inside  of  the  tank  was  baffled  in  order  to  aid  the 
separation  process.  A  sectional  elevation  and  a  plan  view 
of  the  separator  are  given  in  Figure  21.  The  baffles  were 
cylindrical,  the  outer  one  being  42  inches  in  diameter  and  52 
inches  high  and  the  inner  ones  being  27  inches  in  diameter  and 
26  inches  high.  The  two- phase  mixture  entered  the  separator 
tangentially  at  the  mid-point  of  the  vertical  side  through  a 
3- inch  diameter  line.  The  tangential  entry  minimized  the 
turbulence  and  consequent  emulsion  formation  which  would  have 
resulted  had  the  entry  been  normal  to  the  periphery  of  the  tank. 
The  separation  of  the  oil  and  water  took  place  in  the  annular 
space  between  the  outer  baffle  and  the  tank  wall.  The  oil 
then  overflowed,  and  the  water  under flowed,  the  outer  baffle 
into  the  relatively  undisturbed  region  inside  the  outer 
baffle.  It  was  found  that  the  oil  tended  to  carry  minute 
droplets  of  water  over  the  outer  baffle  and  these  were  removed 
by  a  bed  of  coarse  sodium  chloride  supported  on  a  perforated 
tray  between  the  outer  and  inner  baffles. 


■ 


FIGURE  21  -  OIL-WATER  SEPARATOR 
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Inside  the  outer  baffle  the  phases  first  flowed  towards 
their  common  interface  and  then  subsequently  away  from  the 
interface  to  the  respective  outlets  -  the  outlet  for  the  oil 
being  at  the  centre  of  the  upper  inner  baffle,  about  five 
inches  below  the  oil  surface  and  the  outlet  for  the  water 
being  at  the  centre  point  of  the  bottom  of  the  tank. 

Bypass  flows  entered  the  separator  tank  tangentially, 

12  inches  above  and  12  inches  below  the  two- phase  entry, 
through  2-inch  diameter  lines.  These  bypass  lines  were  positioned 
so  that  both  liquids  could  be  completely  circulated  to  ensure 
uniformity  of  temperature  without  passing  them  through  the 
rectangular  conduit.  Fresh  make-up  water  was  admitted  through 
three  equally  spaced  1/2- inch  diameter  openings  on  the  periphery 
of  the  tank.  Dial  thermometers  inserted  through  the  side  of 
the  tank  indicated  the  temperatures  of  the  two  phases  and  a 
sight  glass  indicated  the  level  of  the  oil-water  interface. 

The  separator  tank  was  manufactured  completely  from 
14  gauge  galvanised  iron  sheet,  contained  about  300  U.S. 
gallons  of  each  phase,  and  was  positioned  on  a  wooden  plat¬ 
form  about  18  inches  high. 

A  sheet  metal  adapter  was  used  to  convert  the  1-inch 
x  8-inch  conduit  cross-section  into  the  3-inch  diameter  circular 
cross-section  of  the  pipe  which  ran  directly  to  the  separator. 
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This  pipe  contained  a  3- inch  A.C.F.  ball  valve.  This  valve, 
when  fully  open,  provided  a  smooth  passage  for  the  oil  and 
water  leaving  the  conduit.  A  valve  of  the  gate  or  globe 
type  would  have  created  turbulence  tending  to  emulsify  the 
oil  and  water  together. 

Two  bypass  tanks,  one  for  each  phase,  were  installed 
to  allow  the  unwanted  portions  of  the  total  flow  from  each 
pump  to  be  circulated  without  returning  the  flows  to  the 
separator.  High  bypass  flows  returned  directly  to  the  separator 
tended  to  hinder  the  separation  process.  The  tanks  were 
identical  with  capacities  of  approximately  50  U. S.  gallons 
and  had  cylindrical  bodies,  35  inches  high  and  22  inches 
diameter,  and  ends  coned  outwards.  Tanks  with  substantial 
capacities,  rather  than  just  pipe  connections,  were  installed 
so  that  the  rise  in  temperature  of  the  liquids  on  continued 
recirculation  through  the  pumps  would  be  kept  to  a  minimum. 

Both  pumps  were  of  the  positive  displacement  gear  type. 
The  oil  pump,  model  K30  manufactured  by  G.D.  Roper  Corporation, 
had  a  capacity  of  36  U.S.  gallons  per  minute  at  1750  R.P.M. 
and  the  water  pump,  model  FL  manufactured  by  Jabsco  Pump  Co. , 
had  a  capacity  of  30  U.S.  gallons  per  minute  at  1750  R.P.M. 

Both  pumps  were  driven  by  3- phase,  two-horse  power  electric 
motors  running  at  1750  R.P.M. 
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Rotameter  type  flow  meters  had  been  found  to  be 
entirely  satisfactory  in  previous  investigations and 
similar  instruments  were  used  in  the  present  study.  Oil  flow 
rates  in  the  range  0.3  to  36  U.S.  gallons  per  minute  were 
covered  by  three  rotameters  in  which  the  floats  were  of  the 
sharp  edged  type  to  provide  immunity  to  slight  viscosity  varia¬ 
tions.  Water  flow  rates  from  0.3  to  30  U.S.  gallons  per  minute 
were  also  covered  by  three  rotameters.  Since  the  temperature 
coefficient  for  the  viscosity  of  water  is  relatively  small, 
sharp  edged  floats  were  not  used  in  the  water  rotameters  and 
consequently  use  could  be  made  of  smaller  sized  rotameters. 
Details  of  the  rotameters  and  their  calibration  are  given  in 


Appendix  C. 


IV  EXPERIMENTAL  TECHNIQUES  AND  PROCEDURE 

The  independent  variables  in  the  experimental  investi¬ 
gation  were  the  flow  rates  of  the  two  phases  which  were  indi¬ 
cated  by  the  rotameters  described  in  the  previous  section  and 
Appendix  C.  For  a  given  combination  of  flow  rates  the  input 
water-oil  volume  ratio  was,  of  course,  also  fixed.  Under  each 
set  of  input  conditions  it  was  desired  to  determine  the  following 
dependent  variables 

(1)  the  flow  condition,  i.e.  laminar,  transitional  or 
turbulent  in  each  phase, 

(2)  the  insitu  water-oil  volume  ratio, 

(3)  the  pressure  gradient, 

(4)  local  velocities,  and 

(5)  wave  velocity,  wavelength  and  wave  amplitude  for 
conditions  under  which  interfacial  waves  were 
present. 

A.  TECHNIQUES 

1.  Flow  conditions 

Flow  conditions  were  determined  by  injecting  soluble 


dyes  into  the  flowing  fluids  by  means  of  extremely  small 
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diameter  hypodermic  needles  inserted  through  the  top  and 
bottom  of  the  conduit  at  the  upstream  end  of  the  photographic 
section.  Under  laminar  conditions  the  dye  filaments  retained 
their  identity,  while  in  turbulent  flow  the  eddies  caused 
the  filaments  to  break  up  and  the  dye  diffuse  throughout  the 
fluid.  In  general,  the  dye  traces  were  observed  visually 
although  a  few  representative  traces  were  recorded  photogra¬ 
phically. 

2.  Insitu  ratios 

Insitu  ratios  were  determined  from  measurements  of 
the  depth  of  the  oil  phase.  Relatively  large  oil  depths  were 
measured  with  the  aid  of  a  scale,  which  was  read  to  0.005 
inches,  attached  to  the  outside  of  the  conduit.  Small  depths 
were  recorded  by  means  of  a  pointer  gauge  which  consisted 
of  a  hypodermic  needle  attached  to  the  shaft  of  a  micrometer. 
The  micrometer  was  adjusted  until  the  tip  of  the  needle  just 
touched  the  interface  and  was  read  to  0.001  inches. 

For  conditions  when  no  waves  were  present  at  the  inter¬ 
face,  or  waves  of  very  small  amplitude,  it  was  possible  to 
obtain  the  oil  depth  under  flowing  conditions.  However,  when 
pronounced  waves  were  present,  the  oil  depth  was  recorded 
immediately  after  the  flows  of  both  phases  were  stopped  simul- 
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taneously  and  the  A.C.F.  valve  between  the  test  section  and 
the  separator  closed.  In  this  way  the  oil  depth  was  measured 
after  the  waves  had  disappeared  but  before  any  imperfections 
in  the  levelling  of  the  conduit  caused  the  level  of  the  inter 
face  to  change. 

3.  Pressure  gradients 

Pressure  drops  in  earlier  studies (^9 ' -^)  Qf  the  hori¬ 
zontal  pipeline  flow  of  oil-water  mixtures  were  measured  by 
means  of  a  manometer  system  in  which  the  pressure  in  the  pipe 
line  was  transmitted  to  the  manometer  by  air  adjusted  to 
the  correct  pressure  by  pressure  reduction  valves.  Such  a 
system  prevented  oil-water  mixture  being  present  in  the 
manometer  lines,  which  would  have  caused  erroneous  results. 
The  method  was,  however,  cumbersome  in  operation  since  it 
required  an  operator  to  be  present  at  each  pressure  tap. 

Since  the  present  study  was  restricted  to  an  investi¬ 
gation  of  stratified  flow  it  was  not  necessary  to  adopt  the 
method  used  previously.  Pressure  taps  located  in  the  bottom 
of  the  conduit  allowed  water  to  be  used  as  the  fluid  trans¬ 
mitting  pressure  to  the  manometers  with  no  danger  of  the  oil 
entering  the  manometer  lines. 
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Pressure  drops  over  the  20-foot  test  section  ranged 
from  about  0.1  inches  of  water  at  low  flow  rates  to  about  4 
inches  of  water  at  high  flow  rates.  Pressure  drops  up  to  about 
one  inch  of  water  were  measured  by  the  use  of  a  specially 
constructed  micromanometer,  while  greater  pressure  drops  were 
measured  by  means  of  an  inclined  manometer. 

Single  phase  pressure  drops  were  recorded  by  the  use 
of  the  micromanometer  and  inclined  manometer  in  the  case  of 
the  flow  of  water,  and  by  the  use  of  an  inclined  manometer 
only  in  the  case  of  the  flow  of  oil. 

The  micromanometer  is  shown  diagrammatically  in  Figure 
22  and  was  based  on  an  instrument  marketed  by  Flow  Corporation, 
Massachusetts,  which  was  designed  to  measure  pressure  drops 
in  gas  flow.  The  design  was  modified  so  that  it  could  be  used 
as  a  two- liquid  manometer  -  the  liquids  being  carbon  tetra¬ 
chloride  and  water.  The  instrument  consisted  basically  of  a 
U-tube,  one  side  of  which  was  a  4- inch  diameter  reservoir  and 
the  other  a  4mm- diameter  glass  tube  inclined  at  30°  to  the 
horizontal.  A  flexible,  lram  inside  diameter,  polyethylene 
tube  joined  the  reservoir  to  the  inclined  tube,  which  itself 
was  held  in  a  brass  slide  raised  and  lowered  by  adjustment  of 


the  micrometer. 
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FIGURE  22  -  MICROMANOMETER 
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Pressure  drops  were  obtained  in  the  following  manner. 
Initially,  a  "zero"  reading  was  obtained  by  adjusting  the 
micrometer  so  that  the  meniscus  in  the  inclined  tube  just 
touched  the  cross-wire  with  tap  A  open  and  taps  B  and  C 
closed.  Subsequently,  with  tap  A  closed  and  taps  B  and  C 
open  the  meniscus  was  returned  to  the  cross-wire  by  again 
adjusting  the  micrometer.  The  pressure  drop  was  therefore 
represented  by  a  column  of  fluid  having  a  density  equal  to  the 
density  differential  between  the  carbon  tetrachloride  and  water, 
and  a  height  equal  to  the  difference  of  the  micrometer  readings. 

Since  the  meniscus  was  always  returned  to  the  same 
position  in  the  inclined  tube  no  errors  were  introduced  by 
local  variations  in  tube  diameter.  The  only  possible  error 
arose  from  the  very  small  change  in  volume  of  the  flexible 
tube  joining  the  reservoir  with  the  inclined  tube  during  the 
vertical  movement  of  the  tube.  However,  the  large  cross- 
sectional  area  of  the  reservoir  caused  any  variation  in  the 
volume  of  the  flexible  tube  to  be  entirely  negligible.  The 
micromanometer  was  found  to  be  sensitive  to  a  change  in 
pressure  corresponding  to  0.001  inches  of  water. 

The  successful  operation  of  the  instrument  depended 
primarily  on  the  inside  of  the  inclined  tube  being  perfectly 
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clean,  otherwise  the  meniscus  would  not  move  freely  and 
erroneous  readings  were  obtained. 

Pressure  drops  were  also  obtained  with  the  micromanometer 
modified  to  the  form  of  a  water- air  inverted  manometer.  Although 
somewhat  less  sensitive,  this  arrangement  was  found  to  be  less 
troublesome  in  operation  than  the  two- liquid  arrangement. 

The  inclined  manometer  used  for  measuring  pressure  drops 
in  excess  of  one  inch  of  water  was  of  the  inverted  form  and 
had  a  1  in  4  slope.  Air  separated  the  liquid  in  the  two  legs 
of  the  manometer.  In  experiments  with  water  flowing  alone 
and  with  the  two-phase  system  water  was  admitted  to  the  mano¬ 
meter  from  the  test  section,  while  in  experiments  with  oil 
flowing  alone  oil  was  admitted  to  the  manometer. 

The  pressure  transmission  lines  consisted  of  semi-trans¬ 
parent,  flexible  polyethylene  tubing.  Such  tubing  readily 
allowed  any  air  bubbles  to  be  observed  and  removed.  Poly¬ 
ethylene  tubing  was  found  to  be  preferable  to  vinyl  tubing, 
supplied  by  Fisher  Scientific  Co. ,  the  position  of  which  was 
found  to  affect  the  pressure  differential  indicated  by  the 
oil- air  manometer.  The  cause  of  this  interaction  between  the 


vinyl  tubing  and  the  oil  was  not  determined. 
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4.  Local  Velocities 

Plans  were  originally  made  to  determine  the  complete 
velocity  profile  on  the  vertical  centre  line  of  the  conduit„ 

For  laminar  conditions  in  both  phases  the  experimental  velocity 
profiles  could  then  be  compared  with  profiles  similar  to 
those  shown  in  Figure  3.  With  one,  or  both,  of  the  phases  in 
turbulent  flow  the  experimental  measurements  would  yield 
information  concerning  the  nature  of  the  interaction  between 
the  phases.  The  use  of  pitot  tubes  was  rejected  on  the  basis 
that  the  velocity  heads  under  laminar  conditions  would  be  of 
the  order  of  0.01  inches  of  water  and  therefore  too  small  to 
be  measured  accurately  without  extremely  elaborate  pressure 
measuring  equipment.  A  flow  visualization  technique  appeared, 
therefore,  to  be  most  suitable. 

Flow  visualization  techniques  have  been  discussed 
by  Hyzer  and  Hinze  a  Measurement  of  the  velocity 

of  particles  carried  in  a  flowing  fluid  to  determine  velocity 
profiles  has  been  used  by  Van  Driest  t  wh0  was  concerned 
with  turbulent  diffusion  in  open  channels,  Gaffyn  and  Under¬ 
wood^'*),  who  were  interested  in  velocity  profiles  for  liquid 
flow  in  pipes,  and  Nedderman who  measured  velocities  in 
the  wall  region  for  turbulent  liquid  flow  in  pipes.  All 
these  investigators  used  particles  which  did  not  drift  signi- 
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ficantly  under  the  action  of  gravity,  i.e„  the  particle  and 
fluid  densities  were  either  equal  or  the  particle  was  so  small 
that  its  drift  velocity  was  negligible. 

In  the  present  investigation  it  was  proposed  to  use  a 
modified  flow  visualization  method  in  which  a  small  air  bubble, 
injected  through  the  bottom  of  the  conduit,  would  rise  through 
the  flowing  liquids  under  the  action  of  gravity.  At  the  same 
time,  however,  the  bubble  would  be  carried  horizontally  by  the 
moving  fluid.  Stroboscopic  illumination  of  the  air  bubble 
would  allow  the  complete  trajectory  of  the  bubble,  viewed 
from  the  side  of  the  conduit,  to  be  recorded  on  a  single  sheet 
of  photographic  film.  Subsequent  analysis  of  the  trajectory 
would  yield  information  concerning  the  local  horizontal  velocity 
at  each  vertical  elevation. 

A  theoretical  prediction  of  the  trajectory  of  a  small 
air  bubble  having  a  diameter  of  0.5  mm  rising  through  a  para¬ 
bolic  velocity  field  indicated  that  inertia  effects  of  the 
bubble  and  its  associated  fluid  would  be  negligible  for  liquids 
having  viscosities  in  excess  of  about  5  centipoise.  On  the 
other  hand,  the  inertia  effect  for  the  bubble  rising  through 
water  would  not  be  insignificant.  The  theoretical  analysis 
and  the  predicted  bubble  velocities  are  presented  in  Appendix  D» 


' 


110 


In  addition  to  the  doubt  raised  about  the  suitability 
of  the  method  by  the  theoretical  analysis  experimental  diffi¬ 
culties  were  experienced  in  obtaining  suitable  photographs. 

While  reasonably  good  photographs  were  obtained  with  an  air 
bubble  rising  through  a  clear  single  phase,  it  was  found  that 
in  the  two-phase  system  the  two  liquids  contaminated  each 
other  slightly  and  neither  was  perfectly  clear.  Consequently 
the  definition  of  the  photographic  image  was  extremely  poor. 
Further,  in  the  two  phase  system,  the  air  bubble  was  found  to  be 
held  at  the  interface  by  interfacial  tension,  thus  preventing 
the  photographing  of  the  trajectory  within  the  upper  phase. 

It  was  concluded  that  the  method  would  be  very  suitable 
for  single  phase  flow  of  a  liquid  with  a  viscosity  greater 
than  about  5  centipoise  but  as  far  as  the  present  project  was 
concerned  the  method  was  felt  to  be  unsuitable. 

Velocity  measurements  were,  therefore,  restricted  to 
measurements  of  the  interface  velocity.  Very  fine  aluminium 
powder,  supplied  by  Fisher  Scientific  Co. ,  was  introduced 
through  the  top  of  the  conduit  and  allowed  to  settle  to  the 
interface.  Measurement  of  the  time  required  for  the  powder 
to  travel  a  known  distance  enabled  the  interface  velocity  to 
be  calculated.  These  rather  limited  velocity  measurements  did, 
however,  present  a  very  clear  picture  of  the  nature  of  the 
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interaction  between  the  phases  and  they  could,  of  course,  be 
compared  with  the  theoretical  values  when  both  phases  were 
in  laminar  flow. 

5.  Wave  Characteristics 

Wave  velocities  were  obtained  simply  by  timing  the 
movement  of  a  particular  wave  over  a  distance  of  three  feet. 
Measurements  of  wavelength  and  amplitude  were  not  so  easily 
perf ormed. 

Several  experimental  investigations  have  been  reported 
in  which  film  thicknesses  or  wave  profiles  were  measured. 

Amongst  these,  Bunt ^  measured  film  thickness  using  fluorescence, 
Dukler  and  Bergelin^3)  use<3  a  capacitometer  in  conjunction 
with  an  oscillograph  to  determine  wave  forms  on  liquids  flowing 
down  inclined  plates  and  Lilleleht  and  Hanratty used  a 
light  absorption  technique  to  determine  wave  profiles  at  the 
interface  for  the  co-current  flow  of  air  and  water  in  a 
rectangular  conduit.  The  use  of  stereophotogrammetry  in 
contouring  ocean  waves  has  been  discussed  by  Schumacher 
and  Brobikov  ^  . 

It  was  decided  to  investigate  the  use  of  stereophoto¬ 
grammetry  to  obtain  wavelengths  and  amplitudes  on  the  relatively 
small  laboratory  scale.  The  method  has  the  advantage  that 
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no  probes  of  any  kind  interfere  with  the  flow  of  the  fluids 
and  the  disadvantage  that  only  a  relatively  small  part  of  the 
interface  can  be  analysed.  The  only  requirement  of  the  method 
is  that  the  interface  must  be  visible  to  the  cameras. 

The  two  cameras  were  identical  small  field  cameras 
manufactured  by  Instruments  Ltd. ,  Ottawa.  The  focal  lengths 
of  the  lenses  were  112  mm  and  in  order  that  the  cameras  could 
be  focused  on  the  interface  special  brass  adapters  were  manu¬ 
factured  to  extend  the  lens  to  photographic  plate  distance. 

The  holders  took  either  3  ^/4-inch  x  4^/4- inch  sheet  film  or 
glass  plates.  Fiducial  marks  within  the  cameras  appeared 
on  the  exposed  film. 

The  cameras  were  arranged  as  shown  in  Figure  20  so 
that  photographs  could  be  taken  from  beneath  the  conduit. 

They  were  rigidly  attached  to  a  1/2- inch  thick  brass  plate 
by  means  of  90°  brackets  and  thumb  screws,  so  that  the  photo¬ 
graphic  planes  were  parallel  to  the  plate.  The  horizontal 
distance  between  the  cameras  could  be  varied  in  1/2- inch 
intervals.  The  brass  plate  was  supported  by  three  vertical 
rods,  the  upper  ends  of  which  were  threaded.  Nuts  on  the 
threaded  parts  beneath  the  plate  allowed  the  vertical  position 
of  the  plate,  and  hence  the  cameras,  to  be  adjusted.  Initially 
the  plate  was  approximately  levelled  using  the  circular  levels 
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attached  to  the  upper  surface  of  the  plate.  The  final  levelling 
was  performed  with  the  aid  of  a  level  sensitive  to  0.001  inches 
over  a  distance  of  one  foot. 

The  positions  of  the  cameras  relative  to  each  other  and 
to  the  interface  are  important  in  that  they  affect  the  para¬ 
llax  readings  which  are  used  in  the  calculation  of  the  eleva¬ 
tion  of  different  points  on  the  interface.  For  a  given 
elevation  difference,  parallax  difference  is  greater  when  the 
cameras  are  further  apart  and  closer  to  the  interface.  However, 
under  such  conditions  the  overlapping  field  of  view  is 
relatively  small.  Consequently  it  is  necessary  to  compromise 
accuracy  and  field  of  view.  In  order  to  obtain  an  overlapping 
field  of  view  amounting  to  about  six  inches  in  the  direction 
of  flow,  the  cameras  were  positioned  approximately  27  inches 
below  the  conduit  and  with  about  6.5  inches  between  the  optical 
axes.  With  this  arrangement  it  was  found  that  only  the 
wavelength  and  amplitude  of  waves  having  relatively  large 
amplitudes  could  be  obtained  accurately.  Had  cameras  been 
available  with  wider  fields  of  view  then  a  higher  degree  of 
accuracy  could  have  been  obtained  by  positioning  the  cameras 
further  apart  and  closer  to  the  interface.  Nevertheless, 
valuable  information  concerning  the  larger  amplitude  waves  was 


obtained. 
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Photographs  were  obtained  on  Kodak  Tri-X  glass  plates 
which  are  extremely  insensitive  to  ambient  temperature  and 
humidity  changes ^ .  Ordinary  flexible  sheet  film  was  found 
to  be  completely  inadequate  owing  to  its  dimensional  instability 
and  to  the  fact  that  it  would  not  lie  perfectly  flat  in  the 
holders.  The  interface  was  made  visible  by  injecting  aluminium 
powder  into  the  oil  phase  such  that  it  settled  to  the  inter¬ 
face  before  reaching  the  photographic  section.  Simultaneous 
photographs  were  obtained  by  using  an  electronic  flash  with 
1/1200  seconds  duration  and  200  watt- secs  power  and  by  having 
the  camera  shutters  open  and  the  laboratory  in  otherwise  total 
darkness.  The  camera  aperture  setting  was  f/16. 

The  stereo  pairs  of  negatives  were  placed  on  a  light 
table  and  viewed  under  a  Wild  ST3  mirror  stereoscope. 

Parallax  differences  were  obtained  by  the  use  of  a  parallax 
bar. 

Complete  details  of  the  theory  of  the  method  and  the 
calibration  procedure  are  given  in  Appendix  E. 
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B.  PROCEDURE 

The  conduit  was  periodically  cleaned  by  closing  the 
valves  at  both  ends  of  the  conduit,  removing  the  photographic 
section  and  the  end  plate  of  the  entry  section,  and  pushing 
absorbant  material  through  the  conduit  on  the  end  of  a  12- foot 
rod.  No  detergent  was  used  in  the  cleaning  operation  because 
of  the  profound  effects  which  the  presence  of  even  a  very  small 
amount  of  surface  active  agent  may  have  on  the  stability  of 
a  fluid-fluid  interface. 

After  replacing  the  removable  section  and  the  end 
plate  on  the  entry  section  adjustment  of  the  amount  of  water 
in  the  separator  enabled  the  level  of  the  oil  surface  in  the 
separator  to  be  positioned  approximately  five  inches  above  the 
oil  outlet  and  the  level  of  the  oil-water  interface  to  be 
positioned  level  with  the  mid- point  of  the  conduit.  This 
was  done  in  order  to  provide  a  sufficient  head  above  the  oil 
outlet  so  that  air  would  not  be  withdrawn  with  the  oil  and 
also  so  that  approximately  a  constant  back-pressure  would  be 
maintained  on  the  conduit  during  the  experiments. 

In  the  case  of  single-phase  experiments  the  particular 
liquid  was  pumped  slowly  into  the  conduit  with  the  A.C.F. 
valve  between  the  conduit  and  the  separator  closed  and  the 
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vent  on  the  entry  section  open.  When  the  conduit  was  full 
the  vent  was  closed  and  the  A.C.F.  valve  opened.  The  flow 
rate  was  then  adjusted  to  the  desired  value  by  manipulating 
the  rotameter  and  bypass  valves.  Such  a  procedure  prevented 
the  other  phase  from  entering  the  conduit  and  adhering  to  the 
conduit  walls. 

A  similar  procedure  was  followed  in  the  case  of  two- 
phase  experiments.  Water  was  pumped  into  the  conduit  first, 
and  when  the  water  completely  covered  the  bottom  of  the  conduit 
oil  was  pumped  in  until  it  made  contact  with  the  top  of  the 
conduit.  The  vent  was  then  closed,  the  A.C.F.  valve  opened 
and  the  flow  rates  adjusted  to  the  desired  valves.  A  final 
adjustment  was  made  to  the  level  of  the  oil-water  interface 
in  the  separator  so  that  it  coincided  as  nearly  as  possible 
with  the  level  of  the  interface  in  the  conduit. 

At  start-up  the  oil  and  water  had  approximately  the 
same  temperature  as  the  laboratory  which  was  usually  very 
close  to  76°F.  With  continuous  circulation  of  the  liquids 
through  the  pumps  their  temperature  rose  slightly,  the 
temperature  of  the  oil  phase  being  about  0.5°F  higher  than 


that  of  the  water  phase. 
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After  visual  observation  and  the  measurement  of  the 
pressure  gradient  indicated  that  conditions  had  become  steady, 
the  pressure  gradient  was  recorded  over  the  20-foot  test  section, 
the  oil  depth  was  measured  at  the  mid-point  of  the  test  section 
(no  evidence  of  an  interfacial  gradient  was  found  outside  the 
initial  part  of  the  calming  section)  and  the  interface  velocity 
in  the  central  region  of  the  conduit  was  measured  by  timing 
the  movement  of  the  aluminium  dust  over  an  8- foot  length 
in  the  downstream  half  of  the  test  section.  When  waves  were 
present  at  the  interface  their  velocity  was  obtained  by  timing 
their  movement  over  a  3- foot  section  and,  in  certain  cases, 
stereo-pairs  of  photographs  were  made  of  the  interfacial 


structure. 


V  RESULTS  AND  DISCUSSION 


A.  SINGLE- PHASE  FLOW 

The  rectangular  conduit  was  constructed  using  a  wooden, 
internal  spacer  block  8  inches  wide  and  1  inch  deep.  Sub¬ 
sequent  physical  measurements  of  the  width  of  the  conduit 
gave  an  average  value  of  8.01  inches.  Direct  physical  measure¬ 
ments  of  the  exact  depth  of  the  conduit  could,  however,  not 
be  made.  Although  the  sides  of  the  conduit  were  exactly 
1  inch  deep,  the  gasket  material  contributed  an  unknown 
amount  to  the  conduit  depth.  Further,  it  was  found  that  the 
top  and  bottom  of  the  conduit  warped  slightly,  possibly  due 
to  stresses  caused  by  the  tightening  of  the  bolts  or  ambient 
temperature  changes.  In  addition,  the  depth  varied  slightly 
with  the  internal  pressure.  On  applying  the  normal  operating 
pressure  of  approximately  three  feet  of  water  the  top  and 
bottom  moved  apart  by  about  0.01  inches.  In  view  of  the 
fact  that  the  depth  of  the  conduit  was  a  very  important 
variable,  the  pressure  gradient  being  proportional  to  the 
cube  of  the  depth  in  single- phase  flow,  it  was  decided  to 
determine  the  average  depth  over  the  twenty- foot  test  section 
under  operating  conditions  by  measuring  pressure  drops  for  the 
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laminar  flow  of  oil  alone  since  theoretical  expressions  for 
single-phase  laminar  flow  in  non-circular  ducts  have  been 
verified  by  Cornish  and  Sparrow  „ 

Pressure  drops  were  recorded  using  the  inclined  mano¬ 
meter  for  oil  Reynolds  numbers  in  the  range  1250  to  2300, 
for  which  dye  injection  showed  the  oil  to  be  in  laminar  flow, 
and  were  reproducible  to  within  0.5  percent.  Two  sets  of  data 
are  given  in  Table  IF  of  Appendix  F.  Flow  rate  was  carefully 
plotted  against  pressure  gradient  on  logarithmic  paper  and 
a  straight  line  of  unit  slope  drawn  through  the  data. 

Equation  (11),  verified  by  Cornish  ,  represented  the  line 
and  thus  enabled  the  depth  of  the  conduit  to  be  calculated 
as  1.007  inches.  The  aspect  ratio  and  equivalent  diameter 
of  the  conduit  were,  therefore,  7.95  and  1.786  inches  respect¬ 
ive  ly . 

Substitution  of  the  physical  dimensions  of  the  conduit 
into  equation  (12)  gave  the  friction  factor- Reynolds  number 
relationship  as 

f  =  20  -  (7' 

Re  1 

This  relationship  is  shown  in  Figure  23  along  with  the  data 
for  the  oil  in  laminar  flow.  Both  sets  of  data  are  included 
but  only  at  a  Reynolds  number  of  3240  is  there  a  sufficient 
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difference  between  the  friction  factors  for  two  points  to  be 
plotted.  Also  included  in  the  figure  for  comparison  purposes 
are  the  theoretical  relationships  for  laminar  flow  between 
parallel  plates  and  in  circular  pipes. 

Pressure  drops  were  also  recorded  for  the  turbulent 
flow  of  water  at  Reynolds  numbers  up  to  20,000.  The  data, 
tabulated  in  Table  2F,  Appendix  F,  and  plotted  in  Figure  23, 
show  excellent  agreement  with  the  Blasius  equation,  equation 
(17)  . 

Data  were  obtained  in  the  transition  region  for  the 
flow  of  both  fluids.  It  is  apparent  from  Figure  23  that  the 
lower  critical  Reynolds  number  was  2300  for  both  fluids  and 
that  the  transition  to  turbulent  flow  was  complete  at  a 
Reynolds  number  of  3500.  These  values  compare  with  2000 
and  3800  found  by  Cornish for  the  flow  of  water  in  a 
rectangular  conduit  having  an  aspect  ratio  of  2.92,  and  with 
2700  and  3300  obtained  by  Walker  et  al for  t^e  f iow  Qf 
water  in  a  conduit  having  an  aspect  ratio  of  20. 
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B.  TWO- PHASE  FLOW 

1.  Flow  regimes  and  interfacial  structure 

The  flow  regimes,  i.e.  the  ranges  of  input  conditions 
for  which  the  phases  were  in  laminar  or  turbulent  flow,  were 
determined  by  injecting  dye  into  each  phase  through  very  fine 
hypodermic  needles.  The  interfacial  structure  was  observed 
visually  by  directing  a  flood  lamp  obliquely  at  the  interface. 
In  this  way  the  slightest  interfacial  disturbances  were  readily 
visible. 

In  order  to  delineate  the  regimes  and  indicate  the 
different  forms  of  interfacial  structure  Figure  24  was 
prepared  in  which  the  superficial  oil  Reynolds  number  is  the 
ordinate  and  the  superficial  water  Reynolds  number  is  the 
abscissa.  The  Reynolds  numbers  are  superficial  in  the  sense 
that  they  are  based  on  the  average  velocities  the  phases  would 
have  if  they  flowed  in  the  absence  of  the  second  phase,  and 
on  the  equivalent  diameter  of  the  conduit.  Since  the  physical 
properties  of  the  fluids  were  constant  the  Reynolds  numbers 
were  simply  proportional  to  the  flow  rates  of  the  phases.  The 
Reynolds  number  approach  is,  however,  useful  because  it  allows 
a  direct  comparison  to  be  made  with  the  single-phase  results. 
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FIGURE  24  -  FLOW  REGIMES  AND  INTERFACIAL  STRUCTURE 
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In  the  two- phase  experiments  oil  Reynolds  numbers  ranged 
from  46  to  3,330  and  water  Reynolds  numbers  from  617  to  22,000. 
The  single- phase  experiments  indicated  that  the  lower  critical 
Reynolds  number  was  2300.  Hence,  for  purposes  of  comparison, 
a  vertical  broken  line  is  drawn  on  Figure  24  at  a  water  Reynolds 
number  of  2300,  and  a  similar  horizontal  line  for  the  oil 
phase. 

Both  phases  were  in  laminar  flow  at  relatively  low 
Reynolds  numbers.  Under  these  conditions  no  waves  were 
perceptible  at  the  interface,  although  slight  inter facial 
disturbances  were  caused  by  the  presence  of  the  meniscus 
at  the  sides  of  the  conduit.  The  meniscus  was  never  perfectly 
smooth  and  the  slight  roughnesses  caused  disturbances  to 
spread  across  the  interface  downstream  of  the  roughness.  These 
disturbances  were,  however,  damped  out  within  two  or  three 
feet  of  their  source. 

Conditions  within  the  two  phases  in  the  laminar- laminar 
regime  are  illustrated  in  Figure  25  (i)  for  ReQ  =  649  and 
Rew  =  1225.  Dye  filaments  injected  simultaneously  into  the 
upper  oil  phase,  which  shows  as  the  darker  band  in  the  figure, 
and  the  water  phase  both  retained  their  identity  and  showed 


the  unidirectional  nature  of  the  motion. 
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FIGURE  25  -  DYE  TRACES 
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Consider  now  a  steady  increase  in  the  water  flow  rate 
at  a  constant  oil  flow  rate  with  ReQ  =  649 *  At  a  water  Reynolds 
number  of  2150  the  water  phase  showed  the  first  signs  of 
turbulence,  while  the  oil  phase  remained  in  laminar  motion* 
Turbulent  patches  in  the  water  phase  travelled  along  the  con¬ 
duit  and  were  separated  by  sections  in  which  the  flow  was 
laminar.  Slight  increases  in  the  water  flow  rate  caused  the 
turbulent  patches  to  be  generated  more  and  more  closely  to¬ 
gether  until  the  transition  to  turbulence  was  complete* 

These  observations  of  the  initiation  of  turbulence  are  quali¬ 
tatively  very  similar  to  those  of  Rotta^7^  for  the  flow  of 
a  single  phase  in  a  circular  pipe.  No  attempt  was  made  in 
the  present  investigation  to  evaluate  intermittency  factors* 

It  was  also  observed  that  the  depth  of  the  water  phase  momen¬ 
tarily  increased  slightly  by  0.02-0*03  inches  as  a  turbulent 
patch  passed  an  observation  point. 

Coincident  with  the  turbulent  patches  the  interface 
was  observed  to  be  disturbed  by  two-dimensional  waves,  i*e. 
waves  having  crests  several  times  longer  than  the  wavelength* 

The  wave  crests  were  normal  to  the  direction  of  flow*  The 
waves  had  a  relatively  short  wavelength  of  about  0*5  inches 
and  a  very  small  amplitude.  There  was  absolutely  no  doubt 
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that  the  first  interfacial  waves  were  associated  with  the 
turbulent  patches  in  the  water  phase.  Between  the  turbulent 
patches  the  interface  was  free  of  waves  and  therefore  the 
groups  of  waves  travelled  at  the  same  velocity  as  the  tur¬ 
bulent  patches. 

The  groups  of  waves  were  first  observed  about  6-12 
inches  inside  the  calming  section  and  it  is,  therefore,  reasonable 
to  conclude  that  the  turbulent  patches  in  the  water  phase 
originated  there  too.  The  turbulent  patches,  as  evidenced  by 
the  interfacial  waves,  filled  the  entire  width  of  the  conduit 
soon  after  they  were  formed  and  increased  in  length  as  they 
travelled  along  the  conduit  -  an  observation  also  in  accord 
with  the  observations  of  Rotta  on  initiation  of 

turbulence  in  circular  pipes. 

The  different  conditions  which  existed  in  the  two 

phases  when  the  first  signs  of  turbulence  were  observed  are 

illustrated  in  Figure  25,  (iia)  and  (iib) .  Both  photographs 

were  taken  with  Re  =  649  and  Re  =  2150.  In  (iia)  both 

o  w 

phases  were  in  laminar  flow,  the  dye  traces  being  similar  to 
those  shown  in  (i) ,  and  no  waves  were  present  at  the  inter¬ 
face.  In  (iib)  a  patch  of  turbulence  existed  in  the  water 
phase  and  associated  with  this  turbulence  were  waves  at  the 
interface.  It  is  also  apparent  that  the  oil  was  little 
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affected  by  the  presence  of  the  waves. 

An  increase  of  about  30  percent  in  the  water  flow 
rate  caused  the  turbulence  in  the  water  phase  to  become 
continuous  with  an  uninterrupted  wave  system  at  the  inter¬ 
face. 

The  lower  critical  water  Reynolds  number  was  located 

at  a  number  of  oil  Reynolds  numbers  and  the  data  plotted  in 

Figure  24.  At  all  oil  Reynolds  numbers  the  transition  was 

similar  to  that  already  described  for  an  oil  Reynolds  number 

of  649.  It  is  apparent  that  at  low  oil  flow  rates,  i.e.  small 

oil  depths,  the  lower  critical  water  Reynolds  number  was  greater 

than  that  for  the  flow  of  water  by  itself.  A  similar  effect 

was  noted  in  Gazley's^^  data  for  the  co-current  flow  of  air 

and  water  in  a  circular  pipe  where  the  gas  phase  was 

stabilized  by  the  presence  of  a  moving  water  layer.  Whether 

the  increased  stability  was  simply  a  result  of  the  presence 

of  a  moving  boundary,  or  whether  it  was  due  to  the  flexibility 

of  the  boundary,  which  has  been  shown  to  be  important  by 
.  (5) 

Benjamin'  ,  is  not  clear.  At  ReQ  =  46  the  lower  critical 
water  Reynolds  number  was  increased  to  3000.  As  the  oil 
Reynolds  number  was  increased  the  critical  water  Reynolds 
number  decreased  until  at  ReQ  =  300  the  critical  Reynolds 
number  was  again  2300. 
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On  a  further  increase  in  the  oil  Reynolds  number  the 
trend  of  a  decreasing  critical  water  Reynolds  number  was  tem¬ 
porarily  arrested.  For  oil  Reynolds  numbers  between  400  and 
1000  the  critical  Reynolds  number  remained  at  approximately 
2200.  At  an  oil  Reynolds  number  of  900  the  relative  interface 
position  was  approximately  0.7.  Reference  to  the  velocity 
profiles  presented  in  Figure  4,  which  applies  for  the  particular 
aspect  and  viscosity  ratios  of  the  experimental  investigation, 
indicates  that  for  such  an  interface  position  the  maximum 
velocity  occurred  at  or  near  the  interface.  It  is  suggested 
that  this  is  the  cause  of  the  relative  stability  of  the  two- 
phase  system.  At  lower  oil  Reynolds  numbers  the  maximum 
velocity  occurred  within  the  water  phase. 

Further  increase  in  the  oil  flow  rate  caused  a  sudden 
decrease  in  stability  with  the  phases  being  mutually  unstable 
in  the  region  of  ReQ  =  2300  and  Rew  =  1000.  In  this  region 
the  small  water  depth  made  dye  injection  studies  of  the  condi¬ 
tions  in  the  water  phase  difficult.  Further,  since  the  plastic 
conduit  was  preferentially  wet  by  the  oil,  the  oil  tended  to 
displace  the  water  from  the  bottom  of  the  conduit  close  to 
the  sides.  The  depths  of  the  phases  were,  therefore,  not 
constant.  In  view  of  this  a  tentative  continuation  of  the 
locus  of  the  critical  water  Reynolds  number  is  shown  dotted 
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in  this  region. 

The  locus  of  the  lower  critical  oil  Reynolds  number 
was  determined  in  a  similar  way  to  that  for  water.  Again, 
turbulent  patches  were  observed  first  and  their  effect  on  the 
interfacial  structure  was  to  tend  to  destroy  the  two-dimensional 
nature  of  the  waves.  The  presence  of  the  turbulent  water  phase 
caused  the  critical  oil  Reynolds  number  to  be  less  than  that 
for  single  phase  flow,  the  difference  being  more  pronounced 
at  high  water  flow  rates.  At  a  water  Reynolds  number  of 
22,000  the  critical  oil  Reynolds  number  was  1600, 

Consider  again  conditions  for  which  ReQ  =  649,  The 

appearance  of  the  two-dimensional  waves  remained  much  the 

same  on  increasing  the  water  flow  rate  up  to  the  capacity  of 

the  equipment,  although  the  wavelength  increased  slightly  to 

about  0.75  inches  and  the  amplitude  increased  slightly  too. 

Conditions  within  the  phases  at  Re„  =  649  and  Rew  =  9060 

are  illustrated  in  Figure  25  (iii) .  Under  these  conditions, 

although  the  water  was  highly  turbulent  and  a  continuous  wave 

system  present  at  the  interface,  the  oil  phase  was  still 

in  basically  unidirectional  laminar  motion.  The  general 

appearance  of  the  two-dimensional  waves  is  indicated  in  Figure 

26  for  Re^  =  1000  and  Re  = 
o  w 


6000. 


>  - 
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FIGURE  26  -  SMALL  TWO- DIMENSIONAL  WAVES 
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At  oil  Reynolds  numbers  in  the  neighbourhood  of  1500 
and  with  the  water  phase  in  turbulent  flow  a  definite  oscillatory 
motion  was  observed  in  the  oil  phase.  This  is  illustrated 
in  Figure  25  (iv)  in  which  the  dye  filament  shows  that  the 
interfacial  waves  caused  a  more  profound  effect  on  the  motion 
of  the  oil  phase. 

At  oil  Reynolds  numbers  below  about  500  profound 

changes  were  noted  in  the  nature  of  the  two-dimensional  waves 

as  the  water  flow  rate  was  increased.  Large  two-dimensional 

waves  appeared  with  a  wavelength  of  approximately  3  inches. 

The  original  small  amplitude,  short  wavelength  waves  were 

superimposed  upon  the  larger  waves.  These  pronounced  waves 

are  shown  in  Figure  27  for  Re  =  115.5  and  Re  =  13700.  The 

3  o  w 

wave  crests  ran  across  the  complete  width  of  the  conduit 
and  the  waves  were  very  regular  in  appearance. 

Further  increase  in  the  water  flow  rate,  with  a  conse¬ 
quent  thinning  of  the  oil  layer,  caused  roll  waves  to  appear. 
Again  the  small  amplitude  waves  were  superimposed  upon  them. 

The  roll  waves  moved  relatively  slowly,  had  sharp  irregular 
crests,  and  a  highly  variable  wavelength  up  to  approximately 
3  inches.  The  roll  waves  are  shown  in  Figure  28  for  ReQ  =  46.5 
and  Rew  =  18130.  While  the  term  "roll  waves"  appears  to  most 
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FIGURE  27  -  LARGE  TWO-DIMENSIONAL  WAVES 
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FIGURE  28  -  ROLL  WAVES 
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aptly  describe  the  waves  observed  they  did  differ  from  the 
roll  waves  observed  by  Hanratty  and  Engen^4*^  for  the 
air-water  system.  In  the  present  investigation  none  of  the 
more  viscous  phase  was  actually  picked  up  and  carried 
forward  by  the  less  viscous  phase  as  in  the  air-water  system. 

It  has  already  been  remarked  that  the  transition  to 
turbulent  flow  in  the  oil  phase  tended  to  destroy  the  two- 
dimensional  nature  of  the  inter facial  waves.  As  the  Reynolds 
numbers  of  both  phases  were  increased  in  the  turbulent- turbulent 
regime  the  waves  became  more  and  more  three-dimensional  in 
character.  Thus,  three-dimensional  waves,  with  the  crest 
length  approximately  equal  to  the  wavelength,  were  observed 
when  the  flow  rates  of  both  phases  were  close  to  the  capacities 
of  the  respective  pumps.  The  three-dimensional  waves  are 
illustrated  in  Figure  29  for  ReQ  =  3030  and  Re^7  =  18,130. 

In  order  to  clarify  the  transitional  phenomena  further, 
pressure  gradient  traverses  are  shown  in  Figures  30  and  31. 

In  Figure  30  pressure  gradient  is  shown  as  a  function  of  water 
Reynolds  number  for  ReQ  =  282.  For  water  Reynolds  numbers 
less  than  about  2200  both  phases  were  in  laminar  flow  and  the 
pressure  gradient  showed  a  gradual  increase  with  increasing 
water  Reynolds  number.  At  water  Reynolds  numbers  greater 


than  2200  the  water  was  in  turbulent  flow  with  two-dimensional 
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FIGURE  29  -  THREE-DIMENSIONAL  WAVES 
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FIGURE  30  -  PRESSURE  GRADIENTS  AT  A  CONSTANT  OIL  REYNOLDS  NUMBER 
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FIGURE  31  -  PRESSURE  GRADIENTS  AT  A  CONSTANT  WATER  REYNOLDS  NUMBER 
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waves  at  the  interface  and  the  pressure  gradient  increased 
more  and  more  rapidly  as  the  conditions  were  approached  for 
which  large  two-dimensional  waves  were  present.  Also 
indicated  in  the  figure  is  the  pressure  gradient  for  the  flow 
of  the  oil  alone.  It  is  apparent  that  for  relatively  low 
water  flow  rates  in  the  region  where  the  water  was  in  laminar 
flow  the  pressure  gradient  was  reduced  by  the  addition  of  the 
water. 

The  traverse  shown  in  Figure  31  is  for  a  constant  water 
Reynolds  number  of  18130.  At  low  oil  Reynolds  numbers  roll 
waves  were  present  and  the  pressure  gradient  was  relatively 
high.  As  the  oil  Reynolds  number  was  increased  and  the  inter¬ 
face  became  less  violently  disturbed  the  pressure  gradient 
actually  decreased  slightly.  Further  increase  in  the  oil 
flow  rate  caused  the  pressure  gradient  to  increase  relatively 
sharply  at  the  transition  from  laminar  to  turbulent  flow  in 
the  oil  phase.  The  pressure  gradient  for  the  flow  of  the 
water  alone  is  also  indicated  for  purposes  of  comparison  and 
it  is  apparent  that  the  two-phase  pressure  gradient  always 
exceeded  the  single- phase  pressure  gradient. 

The  general  nature  of  the  observed  interfacial  structure 
agrees  qualitatively  with  the  observations  of  Scott  Russell 
and  Benjamin ^  for  the  flow  of  air  over  liquids  in  that  the 
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first  waves  observed  appeared  to  be  caused  by  the  pressure 
fluctuations  in  the  turbulent  phase  and  the  large  two-dimensional 
waves  appeared  to  be  caused  by  an  instability  of  the  mean  flow* 
However,  for  the  liquid- liquid  system,  the  waves  generated  by 
the  pressure  fluctuations  in  the  less  viscous  phase  had  a  definite 
two-dimensional  character  not  observed  for  the  air-water  system. 

In  addition  to  the  general  agreement  with  the  observa¬ 
tions  of  Scott  Russell  and  Benjamin,  the  transitional  phenomena 
in  the  present  investigation  are  in  accord  with  the  observa¬ 
tions  of  Hanratty  and  Engen  for  the  flow  of  air  and  water 

in  a  rectangular  conduit.  The  observations  of  Hanratty  and 
Engen  were  made  predominantly  in  the  turbulent- turbulent  regime 
and  waves  similar  to  the  two-dimensional  and  three-dimensional 
waves  which  they  observed  were  observed  in  the  present  study. 
Presumably  if  the  water  flow  rate  could  have  been  increased 
further  in  the  present  investigation  roll  waves,  similar  to 
those  observed  by  Hanratty  and  Engen,  would  have  developed. 

It  is  also  interesting  to  note  a  further  point  of  agreement 
between  the  two  systems;  at  the  transition  from  two-dimensional 
to  three-dimensional  waves  thinner  layers  of  the  more  viscous 
phases  were  more  stable  than  thicker  ones. 

The  different  regimes  and  interfacial  structures  will 


now  be  considered  in  more  detail. 
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2.  Laminar- laminar  regime 

In  the  laminar- laminar  regime  no  waves  were  present 
at  the  interface  and  the  experimental  situation  did,  therefore, 
closely  approximate  the  idealized  model  assumed  in  the  fore¬ 
going  theory  for  stratified  flow.  In  the  laminar- laminar 
regime  experimental  measurements  of  the  interfacial  velocity, 
the  insitu  ratio,  and  the  pressure  gradient  reduction  factor 
were  made  for  oil  Reynolds  numbers  from  111.,  5  to  1413,  i .  e. 
over  more  than  a  tenfold  range,  and  water  Reynolds  numbers 
from  approximately  700  up  to  the  critical  Reynolds  number. 

The  theoretical  expressions  were  evaluated  for  the  experimental 
aspect  ratio,  7.95,  and  viscosity  ratio,  5.326,  which  was 
the  average  value  during  the  experiments.  The  theoretical 
quantities  are  tabulated  in  Table  A3,  Appendix  A  and  the 
experimental  data  in  Table  F3,  Appendix  F. 

The  comparison  of  the  theoretical  and  experimental 
quantities  are  shown  in  Figures  32,  33  and  34;  the  interface 
velocities  are  compared  in  Figure  32,  the  insitu  ratios  in 
Figure  33,  and  the  pressure  gradient  reduction  factors  in 
Figure  34.  Both  the  theoretical  and  experimental  interface 
velocities  are  plotted  as  the  dimensionless  ratio  of  w^,  the 
interface  velocity,  to  Vrc,  the  average  velocity  for  the  flow 
of  the  more  viscous  oil  phase  alone  in  the  rectangular  conduit 


■ 


142 


0^0 
r\i  —  _• 


o 


G) 

O 

00 

o 


O 


CD 

t 

O 

iD 

t 

O 

O 


r<~) 

O 


<\J 

o 


o 


o 


u 

< 

cr 


UJ 

D 

_J 

O 

> 

I— 

D 

a 


_j 

O 


cn 

n 


o 


CO 


FIGURE  32  -  COMPARISON  OF  EXPERIMENTAL  AND  THEORETICAL  LAMINAR-LAMINAR  INTERFACE  VELOCITIES 
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FIGURE  33  -  COMPARISON  OF  EXPERIMENTAL  AND  THEORETICAL  LAMINAR- LAMINAR  INSITU  RATIOS 
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FIGURE  34  -  COMPARISON  OF  EXPERIMENTAL  AND  THEORETICAL  LAMINAR- 
LAMINAR  PRESSURE  GRADIENT  REDUCTION  FACTORS 
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under  the  theoretical  two-phase  pressure  gradient,.  In  this 
way  the  pressure  gradient  is  eliminated  as  a  variable  and, 
for  given  aspect  and  viscosity  ratios,  the  velocity  ratio  is 
predicted  by  theory  to  be  a  function  of  the  input  oil  fraction 
only,  as  is  the  case  with  the  insitu  ratio  and  the  pressure 
gradient  reduction  factor. 

The  experimental  data  confirm  the  unique  relationship 
between  the  three  ratios  and  the  input  oil  fraction,  or  input 
ratio,  and,  in  addition,  verify  the  trends  of  the  theoretical 
curves.  It  will  be  recalled  that  this  was  not  so  for  the 
experimental  data  of  Russell  et  al  for  stratified  flow  in  a 
circular  pipe  shown  in  Figure  1.  However,  small  differences 
between  the  theoretical  and  experimental  quantities  do  exist. 
The  experimental  velocity  ratios,  while  agreeing  very  well 
with  the  theoretical  values  at  both  relatively  high  and  low 
input  fractions,  are  about  4  percent  low  for  input  fractions 
in  the  range  0.4  to  0.7.  The  insitu  ratios  agree  well  at 
high  input  ratios  but  the  experimental  values  are  about  5 
percent  high  at  low  input  ratios.  The  pressure  gradient 
reduction  factors  are  consistently  approximately  10  percent 
lower  than  the  theoretical  values. 

While  experimental  errors  and  slight  variation  in  the 
viscosity  ratio  no  doubt  account  for  the  slight  scatter  of 
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the  experimental  data,  it  is  unlikely  that  they  account  for 
the  deviations  noted  above.  The  following  departures  from 
the  theoretical  model  could  possibly  have  contributed  to  these 
deviations : 

(1)  the  presence  of  a  large  interfacial  viscosity, 

(2)  not  fully  developed  velocity  distribution  and 

(3)  the  presence  of  a  meniscus 

Very  high  interfacial  viscosities  have  been  observed 
in  some  fluid-fluid  systems  (^0)  .  The  effect  of  a  high  inter¬ 
facial  viscosity  would  tend  to  reduce  the  interface  velocity 
at  the  centre  line  of  the  conduit,  especially  in  conduits 
of  small  aspect  ratio  where  the  velocity  distribution  is  not 
two-dimensional  in  the  central  region.  This  effect  would, 
however,  not  be  limited  to  certain  interface  positions  and 
would  not  be  expected  to  affect  the  interface  velocity  in  the 
central  region  of  a  conduit  with  the  relatively  large  aspect 
ratio  of  7.95:1.  It  is,  therefore,  unlikely  that  high  inter¬ 
facial  viscosity  effects  account  for  the  deviations. 

Despite  the  provision  of  a  long  calming  section,  the 
length  of  which  was  determined  from  the  experimental  data  of 
Eckert  and  Irvine  for  single  phase  flow  in  a  rectangular 

conduit,  it  is  possible  that  the  velocity  distribution  was 
not  in  fact  fully  developed.  In  general,  if  this  were  so. 
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experimental  velocity  measurements  in  high  velocity  regions 
would  tend  to  be  lower  than,  and  experimental  pressure  gradients 
higher  than  the  theoretical  values.  These  two  effects  are, 
in  fact,  apparent  in  the  data.  However,  the  deviation  of 
the  insitu  ratio  data  from  the  theory  is  contradictory.  At 
low  water  to  oil  input  ratios  the  maximum  velocity  occurred 
in  the  oil  phase.  Under  these  conditions,  if  the  velocity 
distribution  was  not  fully  developed  the  tendency  would  be 
for  the  depth  of  the  water  phase  to  be  less  than  the  theoretical 
value  and  thus  the  experimental  insitu  ratios  to  be  less 
than  the  theoretical  values  rather  than  the  converse. 

The  presence  of  a  meniscus  at  the  sides  of  the  conduit 
was  a  definite  departure  from  the  theoretical  model.  Although 
the  plastic  walls  of  the  conduit  were  preferentially  oil  wet 
the  direction  of  the  curvature  depended  on  the  previous  history 
of  the  flow  rates.  If,  for  example,  the  water  flow  rate  was 
increased  to  a  new  value  the  interface  curved  down  to  the 
sides  of  the  conduit.  Data  were  obtained  with  the  meniscus 
curved  both  upwards  and  downwards  and  while,  no  doubt,  the 
curvature  contributed  to  the  scatter  of  the  data,  it  could 
not  account  for  the  overall  deviations.  However,  the  dis¬ 
turbances  caused  by  the  roughnesses  in  the  meniscus  resulted 
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in  energy  losses  which  certainly  caused  the  pressure  gradient 
to  exceed  the  theoretical  value,  although  they  would  probably 
have  little  or  no  effect  on  the  interfacial  velocity  and 
insitu  ratio.  The  lower  pressure  gradient  reduction  factors 
are,  therefore,  at  least  partially  accounted  for  by  these 
disturbances . 

It  would  appear,  therefore,  that  no  single  explanation 
can  be  advanced  to  explain  all  the  small  deviations  observed. 
Despite  these  deviations,  the  experimental  data  verify  the 
theory  reasonably  well  and  the  numerical  data  evaluated  from 
the  theoretical  expressions  can  be  used  with  confidence. 
Further,  the  numerical  results  obtained  by  Charles  and  Red- 
berger '  '  and  Gemmell  and  Epstein1  ;  for  laminar,  stratified 

flow  in  a  circular  pipe  should  also  represent  the  practical 
situation  well,  provided  the  velocity  distribution  is  fully 
developed  and  the  pipe  diameter  is  sufficiently  large  for 
meniscus  effects  at  the  pipe  wall  to  be  insignificant. 

Lockhart  and  Martinelli correlated  pressure  gradient 
data  for  the  flow  of  gas- liquid  mixtures  through  horizontal 
circular  pipes  by  the  four  flow  regimes  and  the  parameters 
defined  by  equations  (3)  and  (4)  in  the  introduction.  As 
was  pointed  out  in  the  introduction  subsequent  experimentally 
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measured  pressure  gradients  for  stratified  flow  in  the  laminar- 
laminar  and  laminar- turbulent  regimes  have  been  less  than 
those  calculated  from  the  relevant  Lockhart-Martinelli  curves. 
Nevertheless,  it  is  interesting  to  explore  what  use  could 
be  made  of  the  basic  Lockhart-Martinelli  parameters  in  correl¬ 
ating  pressure  gradient  data  for  stratified  flow. 

When  the  two  phases  in  stratified  flow  are  considered 
as  the  more  viscous  and  less  viscous  phases,  then  the  Lockhart- 
Martinelli  parameter  0^  ,  defined  by  equation  (3)  is  simply 
the  reciprocal  of  the  pressure  gradient  reduction  factor 

defined  in  the  present  study.  Further,  the  Lockhart-Martinelli 
2 

X  ,  defined  by  equation  (4) ,  is  equivalent  to  the  reciprocal 


of  the  ratio 
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where 


is  the  pressure  gradient 
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for  the  flow  of  the  water  alone  and 


gradient  for  the  flow  of  the  oil  alone, 


is  the  pressure 
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For  laminar  flow  of  both  phases,  it  is  readily  shown 
by  the  use  of  equations  (12) ,  (13)  and  (14)  that 
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Obviously,  therefore,  the  pressure  gradient  reduction  factor 
and  the  ratio  of  the  single  phase  pressure  gradients  could 
not  correlate  data  for  different  aspect  ratios.  Further,  the 


complexity  of  the  curves  presented  in  Figures  6,  7,  8,  and  9 
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in  which  the  pressure  gradient  reduction  factor  is  plotted 
against  the  interface  position  for  different  aspect  ratios 
indicates  that,  even  for  a  given  aspect  ratio,  a  Lockhart- 
Martinelli  type  correlation  would  be  inadequate  over  the 
complete  range  of  input  conditions.  However,  it  is  evident 
from  these  figures  that  the  pressure  gradient  reduction 
factor  tends  to  become  independent  of  the  viscosity  ratio  at 
low  relative  interface  positions.  Under  such  conditions  an 
analysis  of  the  numerical  data  indicates  that  the  input  ratio 
is  approximately  proportional  to  the  viscosity  ratio  and, 
therefore,  a  Lockhart-Martinelli  type  correlation  would 
represent  the  data  well.  Low  relative  interface  positions 
correspond  to  large  input  fractions  of  the  less  viscous  phase 
Equation  (78)  would  be  unchanged  for  laminar- laminar 
flow  in  a  circular  pipe.  The  numerical  results  obtained  by 
Charles  and  Redberger ( -^)  and  Gemmell  and  Epstein ^7)  for 
this  flow  situation  also  show  that  at  high  input  fractions 
of  the  less  viscous  phase  the  input  fraction  is  again  approxi 
mately  proportional  to  the  viscosity  ratio.  Further,  the 
numerical  results  were  independent  of  the  pipe  diameter. 

There  does,  therefore,  appear  to  be  definite  theoretical 
justification  for  the  Lockhart-Martinelli  correlation  in  the 
stratified  laminar- laminar  regime  provided  the  input  fraction 
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of  the  less  viscous  phase  is  relatively  high. 

3.  Laminar- turbulent  regime 

When  the  oil  was  in  laminar,  and  the  water  in  turbulent 
flow,  waves  were  always  present  at  the  interface,  the  three 
wave  types  observed  in  this  regime  being  small  two- dimensional 
waves,  large  two-dimensional  waves,  and  roll  waves.  Turbulence 
in  the  water  phase  established  a  more  uniform  velocity  distri¬ 
bution  and  increased  the  rate  of  dissipation  of  energy  in 
the  water  phase  compared  with  the  corresponding  laminar  flow. 
The  effect  which  these  two  factors  had  on  the  experimental 
insitu  ratios,  velocity  ratios,  and  pressure  gradient  reduction 
factors  will  be  illustrated  by  plotting  representative  data 
on  the  same  basic  co-ordinates  as  were  used  to  represent  the 
laminar- laminar  data. 

(i)  Insitu  ratios 

It  has  been  remarked  that,  on  the  transition  to  tur¬ 
bulence  in  the  water  phase,  the  level  of  the  interface  in  a 
turbulent  patch  was  higher  than  the  level  for  the  corresponding 
laminar- laminar  flow,  i.e.  the  insitu  water  to  oil  ratio  was 
greater  for  laminar- turbulent  flow  than  for  the  corresponding 
laminar- laminar  flow.  This  effect,  which  was  a  consequence 
of  the  increased  rate  of  energy  dissipation  in  the  water  phase. 


152 


was  apparent  throughout  the  laminar- turbulent  regime .  To 
illustrate  this,  representative  data  are  plotted  in  Figure  35. 
The  data  for  ReQ  =  649  are  typical  of  the  results  obtained 
for  conditions  when  small  two-dimensional  waves  were  present 
at  the  interface.  The  insitu  ratios  scattered  slightly  about 
a  curve  parallel  to,  and  about  20-25  percent  above,  the 
theoretical  laminar- laminar  curve  and  showed  no  definite  trend 
with  oil  Reynolds  number. 

At  oil  Reynolds  numbers  for  which  the  large  two-dimen¬ 
sional  and  roll  waves  were  observed,  the  insitu  ratios  showed 
a  further  increase  at  the  transition  to  large  two-dimensional 
waves,  reflecting  the  increased  shear  force  at  the  interface 
caused  by  the  presence  of  the  larger  waves.  For  oil  Reynolds 
numbers  in  the  range  46.5  to  183  a  definite  effect  of  oil 
Reynolds  number  existed.  At  a  particular  oil  Reynolds  number 
the  insitu  ratio  increased  rapidly  with  increase  in  the  input 
ratio  throughout  the  region  for  which  large  two-dimensional 
waves  were  present.  At  oil  Reynolds  numbers  of  46.5  and  73.8 
insitu  ratios  obtained  with  roll  waves  at  the  interface 
increased  less  rapidly  with  increasing  input  ratio. 


FIGURE  35  -  LAMINAR-TURBULENT  INSITU  RATIOS 
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(ii)  Interface  velocities 

Some  laminar- turbulent  velocity  ratios  are  compared 
with  the  theoretical  laminar- laminar  curve  in  Figure  36.  The 
pressure  gradient  used  to  calculate  the  denominator  of  the 
velocity  ratio  was  again  computed  from  the  laminar- laminar 
theory.  While  the  velocity  ratio  is  now  a  fictitious  quantity 
as  far  as  the  experimental  data  are  concerned  it  does  allow 
a  valid  and  interesting  comparison  to  be  made  with  the  data 
previously  presented  for  laminar- laminar  flow.  Two  definite 
effects  are  apparent  in  the  data.  At  low  input  volume  fractions 
the  interface  velocities  were  higher  than  the  laminar- laminar 
values  and  at  relatively  high  input  fractions  the  converse  was 
true.  These  deviations  are  readily  explained  with  reference 
to  the  insitu  ratio  data  just  presented.  At  oil  Reynolds 
numbers  of  46.5  and  111.5  the  oil  was  present  as  a  relatively 
thin  layer  and  the  maximum  velocity  occurred  in  the  water 
phase.  The  increased  insitu  ratio  under  laminar- turbulent 
conditions  necessitated  that  the  actual  average  velocity  of 
the  oil  phase  be  greater  than  the  laminar- laminar  value. 
Consequently  a  higher  interface  velocity  was  observed.  At 
oil  Reynolds  numbers  of  649  and  930  the  maximum  velocity 
occurred  in  the  oil  phase  when  the  transition  to  turbulence 
took  place  in  the  water  phase.  The  increase  in  insitu  ratio 
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which  accompanied  the  transition  lowered  the  actual  average 
water  velocity  and  thus  caused  a  lower  interface  velocity. 

With  increasing  water  fraction  at  a  constant  oil  Reynolds 
numbers  of  649  and  930  the  oil  layer  became  progressively 
thinner  and  eventually  interface  velocities  greater  than  the 
theoretical  laminar- laminar  values  were  again  observed. 

It  would  appear  from  the  data  that,  at  the  transition 
to  turbulence  in  the  water  phase,  either  an  increase  or  a 
decrease  in  the  interface  velocity  relative  to  the  laminar- 
laminar  value  is  possible,  depending  upon  which  phase  contained 
the  maximum  velocity.  The  deviations  of  the  interface  velocities 
from  the  laminar- laminar  curve  were  a  direct  consequence  of 
the  increased  interfacial  shear  force  and  the  more  uniform 
nature  of  the  velocity  distribution  under  turbulent  conditions 
in  the  water  phase.  It  is  important  to  note  that  for  laminar- 
laminar  flow  the  maximum  in  the  curve  of  velocity  ratio  versus 
oil  input  fraction  indicates  the  oil  input  fraction  for  which 
the  interface  velocity  is  a  maximum  for  a  given  pressure 
gradient  but  this  is  not  so  for  the  experimental  data  in  the 
laminar- turbulent  regime  because,  for  a  given  interface 
position,  the  flow  rates  are  no  longer  directly  proportional 
to  the  pressure  gradient. 
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(iii)  Pressure  gradients 

Representative  laminar- turbulent  pressure  gradient 
data  are  plotted  in  Figure  37  in  terms  of  the  pressure  gradient 
reduction  factor  and  the  oil  input  volume  fraction.  The 
laminar- laminar  theoretical  curve  and  the  position  of  the 
laminar- laminar  experimental  data  are  included  for  comparison. 

The  onset  of  turbulence  in  the  water  phase  caused  a  sharp  in¬ 
crease  in  the  two-phase  pressure  gradient  with  a  consequent 
decrease  in  the  pressure  gradient  reduction  factor.  In  the 
laminar- turbulent  regime,  for  given  aspect  and  viscosity  ratios, 
the  pressure  gradient  reduction  factor  is  no  longer  a  function 
of  the  input  fraction  only,  but  of  the  oil  Reynolds  number 
as  well. 

As  the  oil  input  fraction  decreased  the  data  deviated 
further  and  further  from  the  theoretical  laminar- laminar  curve. 

At  an  oil  Reynolds  number  of  74.8  and  at  an  oil  Reynolds  number 
of  183  with  input  fractions  less  than  about  0.10  large  two- 
dimensional  waves  or  roll  waves  were  present  at  the  interface. 
This  is  reflected  in  the  increased  steepness  of  the  curves 
relative  to  the  other  curves  for  which  the  small  two-dimensional 
waves  were  present  at  the  interface. 

The  increase  in  pressure  gradient  accompanying  the 
transition  to  turbulence  in  the  water  phase  limits  the  fractional 
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FIGURE  37  -  LAMINAR- TURBULENT  PRESSURE  GRADIENT  REDUCTION  FACTORS 
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reduction  in  pressure  gradient,  which  can  be  achieved  for 
the  flow  of  the  more  viscous  phase  at  a  given  rate  by  adding 
a  less  viscous  phase  in  turbulent  flow,  to  a  value  less  than 
the  laminar- laminar  value.  This  observation  invalidates  the 
tentative  conclusion  advanced  by  Charles  and  Redberger^^ 
that  the  high  pressure  gradient  reduction  factors  obtained 
by  Charles  were  caused  by  much  of  the  more  viscous  oil 

phase  being  transported  in  fast  moving  wave  crests  adjacent 
to  the  high  velocity  water  layer. 

All  the  pressure  gradient  data  obtained  in  the 
laminar- turbulent  regime  are  represented  in  Figure  38  by  a 
Lockhart-Martinelli  type  correlation  with  the  pressure  gradient 

bz 

data  fall  very  close  to  a  single  curve,  the  maximum  deviation 
of  any  point  being  approximately  15  percent.  The  slope  of 
the  curve  is  relatively  steep  when  large  two-dimensional  and 
roll  waves  were  present  at  the  interface,  reflecting  the 
higher  energy  losses  under  these  conditions.  The  Lockhart- 
Martinelli  curve  for  the  corresponding  regime  of  gas- liquid 
flow  in  circular  pipes  is  also  plotted  for  comparison.  On 
this  type  of  a  plot  the  experimental  data  for  stratified 
gas-liquid  flow  in  circular  pipes  reported  by  Gazley(36), 
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reduction  factor  plotted  against  the  ratio 
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FIGURE  38  -  CORRELATION  OF  LAMINAR  -TURBULENT  PRESSURE  GRADIENT  DATA 
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Hoogendoorn  and  Govier  and  Omer^^  would  fall  between 

the  Lockhart-Martinelli  curve  and  the  curve  representing  the 
present  data. 

It  is  readily  shown,  by  using  equations  (12) ,  (13) , 

(14)  and  (17)  that  for  the  laminar- turbulent  regime 


£b] 

b  z 

'w 


function  of  aspect 
ratio  and  equivalent 
diameter 


-0.25 


input 

ratio 


0.75 


(79) 

is  actually  a  function  of  the 
o 

aspect  ratio,  the  equivalent  diameter,  the  viscosity  ratio, 

the  input  ratio,  and  the  water  Reynolds  number.  In  the  present 

study  the  aspect  and  viscosity  ratios,  and  the  equivalent 

diameter  were  constant.  Under  these  conditions,  therefore, 

it  appears  that  the  pressure  gradient  reduction  factor  and 

0  75 

the  product  of  the  input  ratio  and  (Rew)  *  should  be 

reasonably  good  correlating  parameters.  Whether  or  not  the 

includes  the  correct  functions  of  aspect 
o 

ratio,  equivalent  diameter  and  viscosity  ratio  can  only  be 
determined  by  further  experiments.  However,  other  variables 
which  may  be  important,  notably  the  density  ratio  and  the 
interfacial  tension,  would  not  be  included  in  such  a  correla¬ 
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tion. 
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(iv)  Wave  properties 

Wave  velocities  were  recorded  for  the  three  different 
wave  types  -  the  small  two-dimensional  waves,  the  large  two- 
dimensional  waves  and  the  roll  waves.  The  coincidence  between 
the  onset  of  turbulence  in  the  water  phase  and  the  initiation 
of  small  two-dimensional  waves  suggested  that  these  waves 
were  generated  and  sustained  by  the  turbulent  pressure  fluctua¬ 
tions  in  the  water  phase.  It  would  be  expected  that  the 
velocity  of  these  waves  would  be  a  function  of  the  linear 
velocity  of  the  pressure  fluctuations  and,  therefore,  of  the 
actual  average  velocity  of  the  water  phase.  The  wave  velocities 
were,  therefore,  plotted  against  the  average  water  velocity, 
as  in  Figure  39  which  includes  data  for  all  wave  types. 
Superimposed  on  the  data  is  a  line  of  unit  slope.  The  data 
for  the  small  two-dimensional  waves  fall  very  close  to  this 
line,  indicating  that  the  wave  velocity  was  essentially  the 
same  as  the  actual  water  velocity.  This  provides  further 
evidence  that  the  small  two-dimensional  waves  were  directly 
related  to  the  pressure  fluctuations  in  the  water  phase.  The 
theoretical  approach  of  Phillips should,  therefore,  be 
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actual  average  water  VELOCITY  (FT/SEC) 


FIGURE  39  -  LAMINAR-TURBULENT  WAVE  VELOCITIES 
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useful  in  characterizing  these  waves.  That  the  waves  observed 
were  two-dimensional  in  character,  and  not  three-dimensional 
as  might  be  expected,  is  probably  explained  by  the  theoretical 
work  of  Jeffreys  and  Squire  who  showed  that  two- 

dimensional  disturbances  are  more  unstable  than  three-dimensional 
ones.  Two-dimensional  waves  are,  therefore,  more  likely  to 
be  observed. 

Only  the  velocity  data  obtained  for  the  large  two-dimen¬ 
sional  waves  and  roll  waves  deviate  significantly  from  the 
line  of  unit  slope.  The  velocities  of  the  large  two-dimensional 
and  roll  waves  were  very  much  less  than  the  average  water 
velocity.  There  is  definite  evidence  in  the  data  that  the 
velocity  of  these  waves  increased  with  the  average  water 
velocity  at  a  constant  oil  Reynolds  number.  It  will  be 
recalled  that  the  small  two-dimensional  waves  remained 
superimposed  on  the  large  two-dimensional  waves  and  roll 
waves,  although  their  two-dimensional  character  was  then  not 
so  obvious.  In  one  instance,  at  ReQ  =  46.5  and  an  average 
water  velocity  of  0.854  ft. /sec.  it  was  possible  to  measure 
the  velocity  of  both  two-dimensional  wave  types.  In  Figure 
39  the  upper  point  at  0.596  ft. /sec.  represents  the  velocity 
of  the  small  waves,  and  the  lower  one  at  0.249  ft. /sec.  the 
velocity  of  the  large  waves. 
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The  large  two-dimensional  waves  appeared  to  be  generated 
by  a  fundamental  instability  of  the  mean  flow.  if  such  is  the 
case  then  theoretical  approaches  of  the  types  offered  by 
Lock  and  Miles  for  the  air-water  system  should  be  useful 

in  predicting  the  initiation  of  these  waves.  For  liquid- liquid 
systems  Keulegans ' s  criterion,  equation  (2),  is  available 
for  predicting  interfacial  instability.  Large  two-dimensional 
waves  were  first  observed  when  the  actual  average  water 
velocity  was  approximately  1  ft. /sec.  and  this  value  compares 
with  0.4  ft. /sec.  calculated  from  Keulegan's  criterion.  The 
large  discrepancy  is,  presumably,  a  consequence  of  the  larger 
density  ratio  of  the  present  study  which  would  tend  to  stabilize 
the  flow,  for  in  other  respects  the  experimental  arrangements 
were  similar. 

The  wavelength  and  amplitude  of  the  large  two-dimensional 
and  roll  waves  were  determined  by  the  use  of  stereophotogrammetry 
the  theory  and  calibration  procedure  for  which  are  described 
in  Appendix  E.  Two  representative  stereo  pairs  of  photographs 
are  reproduced  as  positives  in  Figures  40  and  41.  When  viewing 
the  negatives  of  these,  and  of  the  other,  stereo  pairs  under 
the  stereoscope  a  magnification  of  eight  times  was  used. 
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The  lighter  bands  in  the  photographs  are  caused  by 
the  preferential  reflection  of  light  from  the  aluminium 
powder  at  the  interface  when  it  coincided  with  the  wave 
slopes  facing  the  electronic  flash.  The  flash  was  situated 
to  the  left  of  the  photographs. 

The  regular,  equally  spaced  nature  of  the  large 

two-dimensional  waves  is  evident  in  Figure  40  for  which 

ReQ  =  74.8  and  Re^  =  13700,  and,  in  comparison,  the  irregular 

nature  of  the  roll  waves  is  apparent  in  Figure  41  for  which 

Re_  =  74.8  and  Re  =  18130. 
o  w 

The  wavelengths  and  amplitudes  for  the  different  wave 
types  are  tabulated  in  Table  III.  At  a  given  oil  Reynolds 
number  the  wavelength  tended  to  decrease  from  about  three 
inches  to  two  inches  as  the  water  Reynolds  number  increased, 
although  the  wavelength  was  highly  variable  for  the  roll  waves. 
The  decreasing  wavelength  is  compatible  with  the  increasing 
wave  velocity  previously  noted  in  Figure  39.  The  amplitude 
of  the  waves  also  decreased  slightly  as  the  water  Reynolds 
number  increased,  i.e.  as  the  thickness  of  the  oil  layer 


decreased. 
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Table  III 


Data 

for  Large  Two- 

Dimensional 

and  Roll  Waves 

Superficial 

Superficial 

Wave 

Wavelength 

Wave 

Oil  Reynolds 
Number 

Water  Reynolds 
Number 

Type 

(inches) 

Amplitude 

(mm) 

183 

18130 

large  2-d 

1.8 

183 

15800 

large  2-d 

2.0 

0.21 

183 

13700 

large  2-d 

2.4 

0.36 

111.5 

18130 

roll 

1 . 2- 1 . 6 

— 

111.5 

15800 

large  2-d 

1.8 

- 

111.5 

13700 

large  2-d 

2.0 

1.13 

74.8 

18130 

roll 

0 . 4- 2 . 6 

0.70 

74.8 

15800 

large  2-d 

1.96 

0.60 

74.8 

13700 

large  2-d 

2.2 

0.53 

74.8 

11800 

large  2-d 

3.0 

1.11 

46.5 

18130 

roll 

1. 4-3.0 

*_> 

46.5 

15800 

roll 

2.6 

0.58 

46.5 

13700 

large  2-d 

1.8 

0.61 

46.5 

11800 

large  2-d 

2.0 

0.85 

46.5 

10000 

large  2-d 

2.7 

average  amplitude  for  large  two-dimensional  waves  = 

average  amplitude  for  roll  waves  =  0.64  mm 

0 . 70  mm 

■ 
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4.  Turbulent- turbulent  regime 

With  the  water  in  turbulent  flow  the  transition  to 
turbulence  in  the  oil  phase  caused  a  further  sharp  increase 
in  the  pressure  gradient  and  tended  to  destroy  the  two- 
dimensional  character  of  the  interfacial  waves.  At  the 
highest  flow  rates  attainable  with  the  equipment  the  wave 
pattern  was  essentially  three-dimensional,  the  amplitude  of 
the  waves  being  considerably  greater  than  that  of  the  small 
two-dimensional  waves.  The  transition  to  a  three-dimensional 
wave  structure  was  gradual  and  appeared  to  be  caused  by  the 
interaction  of  the  pressure  fluctuations  in  the  two  phases. 

No  quantitative  measurements  of  the  three-dimensional  wave 
properties  were  obtained. 

Pressure  gradient  data  for  oil  Reynolds  numbers  of 
2330  and  3330  are  plotted,  along  with  the  theoretical  laminar- 
laminar  curve,  in  Figure  42  which  is  directly  comparable 
with  Figure  37.  At  the  lower  of  the  oil  Reynolds  numbers  the 
oil  was  in  laminar  flow  when  flowing  alone  and  consequently 
the  pressure  gradient  reduction  factors  were  relatively  low. 

On  the  other  hand,  at  the  higher  Reynolds  number  the  oil  was 
in  turbulent  flow  when  flowing  alone  and  the  pressure  gradient 
reduction  factors  were  more  favorable.  The  curves  are  relatively 
steep  compared  with  those  in  the  laminar- turbulent  regime,  and 


. 


PRESSURE  GRADIENT  REDUCTION  FACTOR 


FIGURE  42  -  TURBULENT- TURBULENT  PRESSURE  GRADIENT  REDUCTION  FACTORS 
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the  pressure  gradient  reduction  factors  were  always  much  less 
than  the  laminar- laminar  values. 


VI  CONCLUSIONS 


1.  The  theoretical  expressions  obtained  for  the  velocity 
distributions  and  volumetric  flow  rates  for  the  laminar, 
stratified  flow  of  two  immiscible  liquids  in  a  rectangular 
conduit  were  reasonably  well  verified  for  arbitrary  viscosity 
and  aspect  ratios  of  5.326  and  7.95:1  respectively.  The 
computed  theoretical  results  which  cover  viscosity  ratios  from 
1  to  1000  and  aspect  ratios  from  8:1  to  1:3  can  therefore  be 
used  with  confidence.  Theory  predicted  that  substantial 
reductions  in  pressure  gradient  can  be  achieved  for  the  flow 
of  a  viscous  liquid  at  a  given  rate  in  a  rectangular  conduit 
by  the  addition  of  a  second  less  viscous  liquid.  For  a 
viscosity  ratio  of  1000  the  maximum  pressure  gradient  reduction 
factor  was  3.38  for  an  aspect  ratio  of  7.896:1  and  1.13  for 

an  aspect  ratio  1:3. 

2.  The  critical  Reynolds  number  for  single-phase  flow  in  the 
rectangular  conduit  was  2300.  With  the  two-phase  system  the 
transition  to  turbulence  in  the  water  phase  occurred  at  a 
higher  Reynolds  number  in  the  presence  of  a  laminar  oil  layer 
provided  the  input  water  to  oil  ratio  was  relatively  high, 
while  the  transition  in  the  oil  phase  took  place  at  a  lower 
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Reynolds  number  in  the  presence  of  a  turbulent  water  layer. 

3.  Small  two-dimensional  waves  appeared  at  the  oil-water 
interface  upon  the  transition  to  turbulence  in  the  water  phase 
and  the  velocity  of  these  waves  was  essentially  the  same  as 
the  actual  average  water  velocity.  It  appeared  therefore 
that  these  small  two-dimensional  waves  were  generated  and 
sustained  by  the  pressure  fluctuations  in  the  turbulent  less 
viscous  phase. 

4.  At  relatively  low  oil  flow  rates  and  relatively  high  water 
flow  rates  large  two-dimensional  and  roll  waves  were  observed. 
The  velocity  of  these  waves  was  considerably  less  than  the 
actual  average  water  velocity  and  they  appeared  to  be  the 
result  of  a  basic  instability  of  the  mean  flow. 

5.  In  the  laminar- turbulent  regime  the  experimentally  deter¬ 
mined  interface  velocities,  insitu  ratios,  and  pressure  gradient 
reduction  factors  deviated  considerably  from  the  laminar- 
laminar  values.  In  particular,  the  pressure  gradient  reduction 
factors  were  always  much  lower  than  the  laminar- laminar  values. 
Lockhart-Martinelli  type  parameters  were  found  to  correlate 

the  laminar- turbulent  pressure  gradient  data  very  well. 
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6.  The  stereophotogrammetric  technique  yielded  good  estimates 
of  the  amplitudes  of  the  large  two-dimensional  and  roll  waves. 
The  accuracy  of  the  method  could  be  greatly  improved  by  the 
use  of  cameras  having  wider  fields  of  view.  The  cameras  could 
then  be  positioned  closer  to  the  interface  and  farther  apart 
so  that  parallaxes  would  be  increased. 

7.  It  would  be  interesting  to  extend  the  present  investiga¬ 
tion  to  cover  higher  water  flow  rates  so  that  the  transition 

to  dispersed  flow  could  be  located,  to  obtain  data  for  different 
viscosity  and  aspect  ratios,  and  to  predict  theoretically 
the  onset  of  the  large  two-dimensional  waves  from  small  dis¬ 


turbance  principles. 


.  ■ 


VII  NOMENCLATURE 


l 

2a 


coefficients  in  equation  (104) 


width  of  rectangular  conduit 


actual  average  water  velocity  = 


water  flow  rate 


(depth  of  water) (width  of  conduit) 


A 

Ax(n) 
A2  (n) 
AJ  (n) 

A2(n) 

2b 

bl 

b2 

B 

B1 

B2  (n) 
B’  (n) 
B^(n) 


ci 


c 

d 

D 


quantity  defined  by  equation  (85) 
coefficient  defined  by  equation  (50) 
coefficient  defined  by  equation  (51) 
coefficient  defined  by  equation  (59) 
coefficient  defined  by  equation  (59) 
depth  of  rectangular  conduit 

depth  of  lower  less  viscous  phase  in  stratified  flow 

depth  of  upper  more  viscous  phase  in  stratified  flow 

quantity  defined  by  equation  (86) 

coefficient  defined  by  equation  (52) 

coefficient  defined  by  equation  (53) 

coefficient  defined  by  equation  (60) 

coefficient  defined  by  equation  (61) 

coefficients  in  equation  (103) 

quantity  defined  by  equation  (87) 

diameter  of  particle 

diameter  of  circular  pipe 
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De 

equivalent  diameter  of  non-circular  conduit 

defined  as  four  times  the  cross-sectional  area 

divided  by  the  wetted  perimeter 

E 

difference  between  fluid  and  particle  velocities 

f 

Fanning  friction  factor  defined  by  equation  (13) 

g 

gravitational  acceleration 

dimension  conversion  factor 

h 

elevation  of  point  above  datum  plane 

Ah 

a 

actual  elevation  difference 

Ahb 

basic  elevation  adjustment 

Ahc 

calculated  crude  elevation  difference 

Ah 

r 

refraction  elevation  adjustment 

H 

height  of  datum  plane  above  camera  lenses 

holdup  ratio 

=  input  ratio  divided  by  insitu  ratio 

i  = 

integer 

input  ratio  = 

Ql 

Q2 

insitu  ratio 

=  ^1 
b2 

k 

g q 

77"  6z  excePb  i-n  equation  (76)  where  k  =  thermal 

conductivity 

m  = 

viscosity  ratio  ^2.  and  actual  mass  of  particle 

Ml 

m*  = 

virtual  mass  of  particle 

m' 

induced  mass  of  particle 
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M 


number  of  sets  of  data  in  equation  (109) 


n 


N 

P 

P 


(2i+l)TT 

except  in  equations  (10) ,  (11)  and 

2a 

(12)  where  n  =  1,3,5.... 

number  of  sets  of  data  in  equation  (108) 
pressure 

parallax  defined  by  equations  (96)  and  (97) 


pressure  gradient  reduction  factor  = 


q 

Q 

Q1 

Q2 

R 

Re 

Re 


Re 


Re. 


w 


$p] 

/ 

I$e1 

dz 

2 

l*J 

tp 


heat  generated/unit  volume/unit  time 
flow  rate 

flow  rate  of  less  viscous  phase 
flow  rate  of  more  viscous  phase 
density  ratio  (air/water) 

Reynolds  number  for  flow  in  a  circular  pipe  = 
Reynolds  number  for  flow  in  a  non- circular 
conduit  = 

M- 

superficial  oil  Reynolds  number  in  two-phase 

.  PoVoDe 
Uo 


flow  = 


superficial  water  Reynolds  number  in  two- phase 

_  PwVe 

M-w 


flow 


relative  interface  position  = 


b' 


bl  +  b2 


sheltering  coefficient  and  transformed  variable 
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time 

temperature 

average  velocity  for  single-phase  flow  in 
a  circular  pipe  or  rectangular  conduit  defined 
by  flow  rate  divided  by  cross-sectional  area 
of  conduit 

vertical  velocity  of  particle 

actual  average  fluid  velocity  in  two-phase  flow 
defined  by  flow  rate  divided  by  cross-sectional 
area  occupied  by  fluid 

superficial  oil  velocity  in  two-phase  flow  defined 
by  oil  flow  rate  divided  by  cross-sectional 
area  of  conduit 

superficial  water  velocity  in  two-phase  flow 
defined  by  water  flow  rate  divided  by  cross- 
sectional  area  of  conduit 

average  velocity  for  the  flow  of  the  more 
viscous  phase  between  parallel  plates  under  the 
theoretical  two-phase  pressure  gradient 
average  velocity  for  flow  of  more  viscous 
phase  in  a  rectangular  conduit  under  the 
theoretical  two-phase  pressure  gradient 
local  fluid  velocity  in  z-direction 


-  5  »V  ,  1  I 


. 


local  velocity  of  less  viscous  fluid 
local  velocity  of  more  viscous  fluid 
local  interface  velocity 
maximum  velocity  in  parabolic  profile 
local  fluid  velocity  in  z-direction 
particle  velocity  in  z-direction 
variable  defined  by  equation  (27) 

Cartesian  co-ordinate 

Lockhart-Martinelli  parameter  defined  by 
equation  (4) 

Cartesian  co-ordinate 

value  of  y  for  which  E  is  a  maximum 

variable  defined  by  equation  (35) 

Cartesian  co-ordinate 

distances  on  stereo- negatives  defined  in  Figure  49 
intermittency  factor 
constant  defined  by  equation  (2) 
dynamic  viscosity 

dynamic  viscosity  of  fluid  particle 

kinematic  viscosity 

density 

gas  phase  velocity  profile  factor  in  equation  (5) 
Lockhart-Martinelli  parameter  defined  by 


equation  (3) 


. 

' 


181 


Subscripts 

1  =  fluid  1  or  object  1 

2  =  fluid  2  or  object  2 

a  =  air 

f  =  fluid 

g  =  gas 

o  =  oil 

p  =  particle 

w  -  water 

j  =  1  or  2 
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APPENDIX  A 

NUMERICAL  DATA  EVALUATED  FROM  LAMINAR  STRATIFIED  FLOW  THEORY 

Table  A1  contains  the  velocity  data  used  to  prepare 
Figures  4  and  5.  The  data  for  the  vertical  velocity  profiles, 
i.e.  at  x  =  0  appear  in  Table  A1 (a)  and  the  data  for  the 
interfacial  velocity  profiles,  i.e.  at  y  =  0,  appear  in 
Table  A1 (b)  . 

Table  A2  contains  the  numerical  data  used  to  prepare 
Figures  6-15.  The  table  is  subdivided  according  to  aspect 
ratio  and  viscosity  ratio. 

Table  A3  gives  the  numerical  data  for  an  aspect  ratio 
of  7.95:1  and  a  viscosity  ratio  of  5.326.  These  data  were 
used  to  plot  the  theoretical  laminar- laminar  curves  in  Figures 
32-37  and  42. 
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Table  Al (b)  Velocity  ratios,  -  for  interfacial  velocity  profiles  at 
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TABLE  A2 

Theoretical  pressure  gradient  reduction  factors  and 
holdup  ratios  for  viscosity  ratios  from  1  to  1000  and  aspect 
ratios  from  7.896:1  to  1:3. 


Table  A2 (a) 
Aspect  Ratio  =  7.896:1 


(i)  Viscosity 

Ratio  = 

Relative 

Interface 

Insitu 

Position 

Ratio 

1.000 

0 

0.900 

0.1111 

0.300 

0.2500 

0.700 

0.4286 

0.600 

0.6667 

0.500 

1.000 

0.400 

1.500 

0.300 

2.333 

0.200 

4.000 

0.100 

9.000 

(ii)  Viscosity  Ratio 

1.000 

0 

0.950 

0.05263 

0.900 

0.1111 

0.350 

0.1765 

0.800 

0.2500 

0.750 

0.3333 

0.700 

0.4286 

0.650 

0.5335 

0.600 

0.6667 

Pressure 

Gradient 

Reduction 

Factor 

Holdup 

Ratio 

1.000 

0.970 

0.2649 

0.893 

0.4686 

0.782 

0.6480 

0.646 

0.8177 

0.499 

1.000 

0.353 

1.223 

0.217 

1.543 

0.1047 

2.134 

0.0287 

3.776 

1.000 

- 

1.413 

0.4364 

1.532 

0.5995 

1.607 

0.7100 

1.554 

0.8078 

1.453 

0.9084 

1. 326 

1.021 

1. 135 

1.153 

1.038 

1.  311 

Input 

Ratio 


0 

0.02943 

0.1172 

0.2777 

0.5451 

1.000 

1.334 

3.601 

8.536 

33.93 


=  5.259 
0 

0.02297 

0.06661 

0.1253 

0.2019 

0.3028 

0.4376 

0.6207 

0.3740 
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Relative 

Interface 

Insitu 

Input 

Pressure 

Gradient 

Reduction 

Holdup 

Position 

Ratio 

Ratio 

Factor 

Ratio 

0.550 

0.8132 

1.232 

0.891 

1.  505 

0.500 

1.000 

1.748 

0.750 

1.748 

0.450 

1.222 

2.512 

0.615 

2.055 

0.400 

1.500 

3.632 

0.491 

2.455 

0.350 

1.357 

5.567 

0.379 

2.998 

0.300 

2.333 

8.698 

0.230 

3.723 

0.250 

3.000 

14.  33 

0.1951 

4.773 

0.200 

4.000 

2  5.54 

0.1252 

6.384 

0.150 

5.667 

51.38 

0.0702 

9.155 

0.100 

9.000 

132.3 

0.0312 

14.70 

0.050 

19.000 

590.6 

0.00790 

31.09 

(iii) 

Viscosity  Ratio 

=  10 

1.000 

0 

0 

1.000 

0.975 

0.02564 

0.01112 

1.430 

0.4338 

0.950 

0.05263 

0.02939 

1.735 

0.5679 

0.925 

0.08103 

0.05269 

1.370 

0.6498 

0.900 

0.1111 

0.07920 

1.925 

0.7123 

0.375 

0.1429 

0.1097 

1.933 

0.7632 

0.850 

0.1765 

0.1452 

1.920 

0.8227 

0.300 

0.2500 

0.2345 

1.310 

0.9381 

0.700 

0.4236 

0.5327 

1.435 

1.243 

0.600 

0.6667 

1.145 

1.132 

1.713 

0.500 

1.000 

2.431 

0.799 

2.481 

0 . 400 

1.500 

5.643 

0.515 

3.762 

0.300 

2.333 

14.28 

0.289 

6. 121 

0.200 

4.000 

44.44 

0.1278 

11.11 

0.100 

9.000 

241.6 

0.0315 

26.34 

(iv)  Viscosity  Ratio  = 

50 

1.000 

0 

0 

1.000 

— 

0.975 

0.02564 

0.02009 

2.39 

0.7333 

0.950 

0.05263 

0.04490 

2.64 

0.3532 

0.925 

0.08103 

0.07542 

2 . 66 

0.9302 

0.900 

0.1111 

0.1140 

2.61 

1.026 
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Pressure 


Relative 

Interface  Insitu 

Position  Ratio 

Input 

Ratio 

Gradient 

Reduction 

Factor 

Holdup 

Ratio 

0.875 

0.1429 

0.1630 

2.52 

1. 141 

0.850 

0.1765 

0.2267 

2.41 

1.285 

0.800 

0.2500 

0.4150 

2.17 

1.660 

0.700 

0.4286 

1.2  34 

1.68 

2.880 

0.600 

0.6667 

3.332 

1.232 

4.997 

0.500 

1.000 

8.570 

0.852 

8.570 

0.400 

1.500 

22.10 

0.537 

14.73 

0.300 

2.333 

61.00 

0.300 

26.14 

0.200 

4.000 

203.9 

0.131 

50.96 

0.100 

9.000 

1164 

0.0327 

129.3 

0  Viscosity  Ratio  = 

=  100 

1.000 

0 

0 

1.000 

_ 

0.975 

0.02564 

0.02434 

2.78 

0.9491 

0.950 

0.05263 

0.05344 

2.91 

1.015 

0.925 

0.08108 

0.09194 

2.86 

1. 134 

0.900 

0.1111 

0.1445 

2.76 

1.  300 

0.875 

0.1429 

0.2161 

2.63 

1.513 

0.850 

0.1765 

0.3147 

2.50 

1.  783 

0.800 

0.2500 

0.6263 

2.229 

2.505 

0.700 

0.4286 

2.096 

1.710 

4.891 

0.600 

0.6667 

6.050 

1.249 

9.075 

0.500 

1.0000 

16.17 

0.8591 

16.17 

0.400 

1.500 

42.66 

0.5426 

28.44 

0.300 

2.333 

119.5 

0.3001 

51.20 

0.200 

4.000 

403.1 

0.1308 

100.8 

0.100 

9.000 

2317 

0.03188 

257.4 

(vi)  Viscosity  Ratio  =  500 


1.000 

0 

0 

1.000 

- 

0.9875 

0.01266 

0.02252 

3.216 

1.779 

0.9750 

0.02564 

0.04633 

3.285 

1.307 

0.9625 

0.03396 

0.07434 

3.251 

1.  921 

0.9500 

0.05263 

0.1095 

3.191 

2.080 

0.9375 

0.06667 

0.1560 

3.121 

2 . 340 
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Relative 

Interface 

Position 

Insitu 

Ratio 

Input 

Ratio 

Pressure 

Gradient 

Reduction 

Factor 

Holdup 

Ratio 

0.9250 

0.03103 

0.2126 

3.046 

2.623 

0.9000 

0.1111 

0.3762 

2.391 

3.386 

0.3000 

0.2500 

2 . 304 

2.281 

9.216 

0.7000 

0.4286 

3.931 

1.734 

20.96 

0.6000 

0.6667 

27.79 

1.261 

41.63 

0.5000 

1.000 

76.95 

0.3651 

76.95 

0.4000 

1.500 

207.1 

0.5452 

133.1 

0.3000 

2.333 

537.2 

0.3011 

251.7 

0.2000 

4.000 

1997 

0.1311 

499.3 

0.1000 

9.000 

11540 

0.0319 

1232 

(vii)  Viscosity  Ratio 

=  1000 

1.000 

0 

0 

1.000 

— 

0.9375 

0.01266 

0.03493 

3.370 

2.759 

0.9750 

0.02564 

0.07265 

3.366 

2.333 

0.9625 

0.03896 

0.1194 

3.306 

3.064 

0.9500 

0.05263 

0.1780 

3.232 

3.333 

0.9375 

0.06667 

0.2601 

3.153 

3.901 

0.9250 

0.03103 

0.3617 

3.072 

4.461 

0.9000 

0.1111 

0.6644 

2.909 

5.930 

0.3000 

0.2500 

4.  399 

2.233 

17.60 

0.7000 

0.4236 

17.59 

1.737 

41.03 

0.6000 

0.6667 

54.96 

1.263 

32.44 

0.5000 

1.000 

152.9 

0.3656 

152.9 

0.4000 

1.500 

412.7 

0.5455 

275.1 

0.3000 

2.333 

1172 

0.3012 

502.2 

0.2000 

4.000 

3989 

0.1311 

997.4 

0.1000 

9.000  23070 

0.03192 

2563 

(viii) 
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Viscosity  Ratio 


Maximum  Pressure 
Gradient  Reduction 
Factor 


1 

1 

5.259 

1.61 

10 

1.93 

50 

2.66 

100 

2.91 

500 

3.28 

1000 

Table  A2 (b) 

3.33 

Aspect  Ratio  = 

3:1 

(i)  Viscosity  Ratio  = 

1 

Pressure 

Relative 

Gradient 

Interface 

Insitu 

Input 

Reduction 

Holdup 

Position 

Ratio 

Ratio 

Factor 

Ratio 

1.000 

0 

0 

1.000 

— 

0.900 

0.1111 

0.02962 

0.9716 

0.2666 

0.300 

0.2500 

0.1131 

0.8952 

0.4725 

0.700 

0.4236 

0.2734 

0.7827 

0.6500 

0.600 

0 . 6666 

0.5460 

0.6470 

0.8190 

0.500 

1.000 

1.000 

0.5002 

1.000 

0.400 

1.  500 

1.831 

0.3533 

1.  221 

0.300 

2.333 

3.591 

0.2179 

1.539 

0.200 

4.000 

8.466 

0.1057 

2.117 

0.100 

9.000 

33.76 

0.02873 

3.751 

(ii)  Viscosity  Ratio  =  3 

1.000  0 

0 

1.000 

— 

0.900  0.1111 

0.05383 

1.230 

0.4845 

0.300  0.2500 

0.1773 

1.262 

0.7093 

0.700  0.4236 

0.3371 

1. 113 

0.9032 
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Pressure 

Relative 

Gradient 

Interface 

Insitu 

Input 

Reduction 

Holdup 

Position 

Ratio 

Ratio 

Factor 

Ratio 

0.600 

0.6666 

0.7495 

0.9051 

1.124 

0.500 

1.000 

1.413 

0.6300 

1.413 

0.400 

1. 500 

2.731 

0.4617 

1.354 

0.300 

2.333 

6.043 

0.2715 

2.592 

0.200 

4.000 

16.42 

0.1243 

4.  106 

0.100 

9.000 

79.01 

0.03153 

3.779 

(iii)  Viscosity  Ratio  = 

10 

1.000 

0 

0 

1.000 

— 

0.9875 

0.01266 

0.003332 

1.254 

0.2672 

0.9750 

0.02564 

0.01024 

1.429 

0.3993 

0.9625 

0.03396 

0.01883 

1.  553 

0.4332 

0.950 

0.05263 

0.02862 

1 . 640 

0.5433 

0.3975 

0.06667 

0.03944 

1.701 

0.5916 

0.9250 

0.03108 

0.05124 

1.742 

0.6320 

0.9125 

0.09539 

0.06404 

1.769 

0.6679 

0.9000 

0.1111 

0.07792 

1.733 

0.7013 

0.8000 

0.2500 

0.2402 

1.674 

0.9607 

0.7000 

0.4286 

0.5607 

1.402 

1.  303 

0.6000 

0.6667 

1.214 

1.095 

1.821 

0.5000 

1.000 

2.609 

0.7935 

2.609 

0.4000 

1.500 

5.332 

0.5226 

3.333 

0.3000 

2.333 

14.  39 

0.2934 

6. 167 

0.2000 

4.000 

44.06 

0.1323 

11.01 

0.1000 

9.000 

237.2 

0.03263 

26.35 

(iv)  Viscosity 

Ratio  = 

30 

1.000 

0 

0 

1.000 

— 

0.9375 

0.01265 

0.006034 

1.630 

0.4767 

0.9750 

0.02564 

0.01526 

1.  912 

0.5953 

0.9625 

0.03896 

0.02570 

2.057 

0.6597 

0.950 

0.05263 

0.03724 

2.133 

0.7076 

0.9375 

0.06567 

0.05000 

2.171 

0.7500 

0.9250 

0.03103 

0.06417 

2.135 

0.7914 

0.9125 

0.09539 

0.07997 

2.183 

0.8340 
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Relative 

Interface 

Position 

Insitu 

Ratio 

Input 

Ratio 

Pressure 

Gradient 

Reduction 

Factor 

Holdup 

Ratio 

0.9000 

0.1111 

0.09767 

2.170 

0.3790 

0.8375 

0.1263 

0.1175 

2.149 

0.9272 

0.3000 

0.2500 

0.3463 

1.896 

1.  387 

0.7000 

0.4236 

0.9721 

1.534 

2.263 

0.6000 

0 . 6666 

2.452 

1.  171 

3.673 

0.5000 

1.000 

5.926 

0.3354 

5.926 

0.4000 

1. 500 

14.46 

0.5436 

9.643 

0.3000 

2.333 

33.15 

0.3074 

16.  35 

0.2000 

4.000 

122.9 

0.1354 

30.739 

0.1000 

9.000 

689.0 

0.03300 

76.56 

(v)  Viscosity 

Ratio  =  100 

1.000 

0 

0 

1.000 

— 

0.9375 

0.01265 

0.008572 

2.197 

0.6772 

0.9750 

0.02564 

0.01929 

2.422 

0.7523 

0.9625 

0.03396 

0.03157 

2 . 490 

0.8104 

0.950 

0.05263 

0.04597 

2.501 

0.8734 

0.9375 

0.06667 

0.06313 

2.433 

0.9469 

0.9250 

0.03108 

0.08375 

2.460 

1.033 

0.9125 

0.09539 

0.1036 

2.426 

1. 133 

0.9000 

0.1111 

0.1335 

2.386 

1.247 

0.8875 

0.1263 

0.1744 

2 . 342 

1.  376 

0.3000 

0.2500 

0.6871 

1.996 

2 . 748 

0.7000 

0.4236 

2. 373 

1.587 

5.548 

0.6000 

0 . 6666 

6.750 

1.201 

10.13 

0.5000 

1.0000 

17.50 

0.8514 

17.50 

0.4000 

1.500 

44.65 

0.5516 

29.  77 

0.3000 

2.333 

121.  3 

0.3107 

51.99 

0.2000 

4.000 

399.  1 

0.1364 

99.77 

0.1000 

9.000 

2270 

0.03303 

252.3 

(vi)  Viscosity 

1.000 

0.9075 

0.3750 

0.9625 


Ratio  =  300 
0 

0.01265 

0.02564 

0.03396 


0  1.000 

0.01060  2.596 

0.02356  2.684 

0.04015  2.630 


0.3375 

0.9188 

1.031 
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Relative 

Interface 

Position 

Insitu 

Ratio 

Input 

Ratio 

Pressure 

Gradient 

Reduction 

Factor 

Holdup 

Ratio 

0.9500 

0.05263 

0.06216 

2.650 

1. 131 

0.9375 

0.06567 

0.09153 

2.603 

1.  373 

0.9250 

0.03103 

0.1303 

2.561 

1.603 

0.9125 

0.09539 

0.1309 

2.511 

1.335 

0.9000 

0.1111 

0.2456 

2.459 

2.210 

0.3375 

0.1263 

0.3272 

2.406 

2.581 

0.3000 

0.2500 

1.650 

2.027 

6.593 

0.7000 

0.4236 

6.384 

1.603 

14.90 

0.6000 

0 . 6666 

19.02 

1.210 

23.53 

0.5000 

1.000 

50.56 

0.8569 

50.56 

0.4000 

1.500 

130.9 

0.5533 

37.25 

0.3000 

2.333 

353.9 

0.3116 

153.8 

0.2000 

4.000 

1138 

0.1367 

296.9 

0.1000 

9.000 

6739 

0.03313 

754.3 

(vii)  Viscosity  Ratio  = 

1000 

1.000 

0 

0 

1.000 

— 

0.9375 

0.01265 

0.01477 

2.807 

1.167 

0.9750 

0.02564 

0.03529 

2.799 

1.  373 

0.9625 

0.03896 

0.06675 

2.759 

1.713 

0.9500 

0.05263 

0.1153 

2.703 

2.191 

0.9375 

0.06667 

0.1374 

2.655 

2.811 

0.9250 

0.03103 

0.2399 

2.600 

3.575 

0.9125 

0.09539 

0.4303 

2.543 

4.487 

0.9000 

0.1111 

0.6158 

2.433 

5.551 

0.3375 

0.1268 

0 . 3533 

2.431 

6.771 

0.3000 

0.2500 

5.014 

2.038 

20.06 

0.7000 

0.4286 

20.40 

1.609 

47.61 

0.6000 

0.6666 

61.97 

1.213 

92.95 

0.5000 

1.000 

166.3 

0.3579 

166.3 

0.4-000 

1.500 

432.7 

0.5546 

233.4 

0.3000 

2.333 

1190 

0.3120 

510.1 

0.2000 

4.000 

3940 

0.1363 

937.3 

0. 1000 

9.000 

22600 

0.03312 

2511 

(viii) 
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Viscosity  Ratio 


Maximum  Pressure 
Gradient  Reduction 
Factor 


1 

1.0 

3 

1.29 

10 

1.79 

30 

2.185 

100 

2.50 

300 

2.685 

1000 

2.81 

(i)  Viscosity 

Relative 

Interface 

Ratio  = 

Insitu 

Table  A2  (c) 

Aspect  Ratio  =  1:1 

1 

Input 

Pressure 

Gradient 

Reduction 

Holdup 

Position 

Ratio 

Ratio 

Factor 

Ratio 

1.000 

0 

0 

1.000 

— 

0.950 

0.05263 

0.003419 

0.9919 

0.1599 

0.900 

0.1111 

0.03242 

0.9689 

0.2918 

0.850 

0.1765 

0.07135 

0.9336 

0.4043 

0.300 

0.2500 

0.1258 

0.3834 

0.5033 

0.750 

0.3333 

0.1976 

0.3352 

0.5929 

0.700 

0.4236 

0.2899 

0.7754 

0.6763 

0.650 

0.5385 

0.4072 

0.7108 

0.7562 

0.600 

0.6667 

0.5567 

0 . 6426 

0.8351 

0.550 

0.8182 

0.7490 

0.5719 

0.9154 

0.500 

1.000 

1.000 

0.5001 

1.000 

0.450 

1.222 

1.  335 

0.4283 

1.092 

0.400 

1.500 

1.796 

0.3577 

1.193 

0.350 

1.357 

2.455 

0.2895 

1.322 

0.300 

2.333 

3.450 

0.2248 

1.479 

0.250 

3.000 

5.060 

0.1651 

1.687 

20  3 


Relative 

Interface 

Position 

Insitu 

Ratio 

Input 

Ratio 

Pressure 

Gradient 

Reduction 

Factor 

Holdup 

Ratio 

0.200 

4.000 

7.947 

0.1113 

1.937 

0.150 

5.667 

14.02 

0.06661 

2.47  3 

0.100 

9.000 

30.35 

0.03141 

3.427 

0.050 

19.000 

113.8 

0.008351 

6.252 

(ii)  Viscosity 

Ratio  =  3 

1.000 

0 

0 

1.000 

— 

0.950 

0.05263 

0.01769 

1.114 

0.3361 

0.900 

0.1111 

0.05739 

1.138 

0.5211 

0.350 

0.1765 

0.1173 

1.  121 

0.6649 

0.300 

0.2500 

0.1937 

1.083 

0.7949 

0.750 

0.3333 

0.3070 

1.030 

0.9210 

0.700 

0.4286 

0.4491 

0.9670 

1.048 

0.650 

0.5385 

0.6343 

0.3953 

1.179 

0.600 

0.6667 

0.8777 

0.8179 

1.317 

0.550 

0.3132 

1.193 

0.7351 

1.465 

0.500 

1.000 

1.628 

0 . 6484 

1.628 

0.450 

1.222 

2.216 

0.5594 

1.313 

0.400 

1.500 

3.048 

0.4696 

2.032 

0.350 

1.357 

4. 273 

0.3307 

2.303 

0.300 

2.333 

6.203 

0.2945 

2.658 

0.250 

3.000 

9.461 

0.2141 

3.154 

0.200 

4.000 

15.62 

0.1422 

3.906 

0.150 

5.667 

29.41 

0.08225 

5.189 

(iii)  Viscosity 

•  Ratio  =  10 

1.000 

0 

0 

1.000 

— 

0.950 

0.05263 

0.02956 

1.305 

0.5616 

0.900 

0.1111 

0.03663 

1.327 

0.7301 

0.850 

0.1765 

0.1763 

1.  296 

1.002 

0.800 

0.2500 

0.3130 

1.243 

1.252 

0.750 

0.3333 

0.5112 

1.173 

1.534 

0.700 

0.4236 

0.7912 

1.104 

1.346 

0.650 

0.5335 

1.17  9 

1.02  3 

2.191 

0.600 

0.6667 

1.713 

0.9343 

2.570 
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Pressure 

Relative 

Gradient 

Interface 

Insitu 

Input 

Reduction 

Holdup 

Position 

Ratio 

Ratio 

Factor 

Ratio 

0.550 

0.3132 

2.447 

0.3412 

2.991 

0.500 

1.000 

3.465 

0.7429 

3.465 

0.450 

1.222 

4.905 

0.6411 

4.013 

0.400 

1.500 

7.005 

0.5376 

4.670 

0.350 

1.857 

10.20 

0.4344 

5.493 

0.300 

2.333 

15.36 

0.3345 

6.  533 

0.250 

3.000 

24.  38 

0.2411 

3.126 

0.200 

4.000 

42.01 

0.1532 

10.50 

0.150 

5.667 

32.81 

0.03989 

14.61 

0.100 

9.000 

208.9 

0.03959 

23.21 

0.050 

19.000 

953.7 

0.009563 

50.19 

(iv)  Viscosity 

Ratio  = 

30 

1.000 

0 

0 

1.000 

— 

0.9875 

0.01265 

0.005373 

1.321 

0.4639 

0.9750 

0.02564 

0.01500 

1.404 

0.5351 

0.9625 

0.03396 

0.02602 

1.436 

0.6678 

0.9500 

0.05263 

0.03924 

1.447 

0.7456 

0.9375 

0.06667 

0.05514 

1.449 

0.8271 

0.9250 

0.03103 

0.07425 

1.446 

0.9157 

0.9125 

0.09539 

0.09710 

1.433 

1.013 

0.9000 

0.1111 

0.1244 

1.429 

1.120 

0.8375 

0.1268 

0.1565 

1.417 

1.2  34 

0.350 

0.1765 

0.2377 

1.377 

1.630 

0.800 

0.2500 

0.5636 

1.312 

2.274 

0.750 

0.3333 

1.0129 

1.239 

3.039 

0.700 

0.4286 

1.677 

1.158 

3.913 

0.650 

0.5335 

2.635 

1.072 

4.394 

0.600 

0.6667 

3.990 

0.9790 

5.935 

0.550 

0.3132 

5.896 

0.8303 

7.206 

0.500 

1.000 

3.537 

0.7778 

3.587 

0.450 

1.222 

12.45 

0.6709 

10.19 

0 . 400 

1.500 

13.17 

0.5622 

12.11 

0.350 

1.857 

26.93 

0.4535 

14.  53 

0.300 

2.333 

41.38 

0.3483 

17.73 

0.250 

3.000 

66.35 

0.2503 

22.28 

0.200 

4.000 

117.3 

0.1636 

29.31 

0.150 

5.667 

235.2 

0.09238 

41. 51 

0.100 

9.000 

603.5 

0.04038 

67.07 
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Pressure 

Relative 

Gradient 

Interface 

Insitu 

Input 

Reduction 

Holdup 

Position 

Ratio 

Ratio 

Factor 

Ratio 

(v)  Viscosity 

Ratio  =  100 

1.000 

0 

0 

1.000 

— 

0.9375 

0.01265 

0.007936 

1.483 

0.6270 

0.9750 

0.02564 

0.01911 

1.523 

0.7451 

0.9625 

0.03396 

0.03427 

1.528 

0.3797 

0.9500 

0.05263 

0.05516 

1.523 

1.043 

0.9375 

0.06667 

0.03355 

1.514 

1.253 

0.9250 

0.08108 

0.1212 

1.503 

1.495 

0.9125 

0.09589 

0.1700 

1.439 

1.773 

0.9000 

0.1111 

0.2324 

1.475 

2.091 

0.8375 

0.1268 

0.3090 

1.459 

2.438 

0.350 

0.1765 

0.6478 

1.411 

3.671 

0.800 

0.2500 

1.432 

1.340 

5.728 

0.750 

0.3333 

2.735 

1.263 

3.203 

0.700 

0.4236 

4.740 

1.179 

11.06 

0.650 

0.5385 

7.689 

1.091 

14.23 

0.600 

0.6667 

11.92 

0.9961 

17 . 33 

0.550 

0.8182 

17.92 

0.8959 

21.90 

0.500 

1.000 

26.46 

0.7910 

26.46 

0.450 

1.222 

38.82 

0.6822 

31.76 

0.400 

1.500 

57.19 

0.5713 

33. 12 

0.350 

1.357 

85.65 

0.4607 

46. 12 

0.300 

2.333 

132.4 

0.3536 

56.73 

0.250 

3.000 

215.4 

0.2538 

71.31 

0.200 

4.000 

330.5 

0.1655 

95.12 

0.150 

5.667 

768.7 

0.0933 

135.6 

0.100 

9.000 

1987 

0.0407 

220.3 

0.050 

19.000 

9299 

0.009709 

489.4 

(vi)  Viscosity 

Ratio  =  300 

1.000 

0 

0 

1.000 

— 

0.9375 

0.01265 

0.009553 

1.562 

0.7547 

0.9750 

0.02564 

0.02503 

1.569 

0.9781 

0.9625 

0.03896 

0.05138 

1.551 

1.319 

0.9500 

0.05263 

0.09371 

1.549 

1.780 
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Relative 

Interface 

Insitu 

Input 

Pressure 

Gradient 

Reduction 

Holdup 

Position 

Ratio 

Ratio 

Factor 

Ratio 

0.9375 

0.06557 

0.1574 

1.535 

2.361 

0.9250 

0.03108 

0.2479 

1.520 

3.057 

0.9125 

0.09539 

0.3705 

1.505 

3.364 

0.9000 

0.1111 

0.5327 

1.489 

4.795 

0 . 837  5 

0.1268 

0.7367 

1.472 

5.812 

0.850 

0.1765 

1.663 

1.421 

9.452 

0.800 

0.2500 

3.387 

1 . 34-8 

15.56 

0.750 

0.3333 

7.643 

1.270 

22.93 

0.700 

0.4286 

13.43 

1.136 

31.45 

0.650 

0.5335 

22.12 

1.096 

41.07 

0.600 

0.6667 

34.55 

1.001 

51.33 

0.550 

0.3132 

52.26 

0.9003 

63.33 

0.500 

1.000 

77.53 

0.7950 

77.53 

0.450 

1.222 

114. 1 

0.6854 

93.33 

0.400 

1.500 

163.7 

0.5739 

112.4 

0.350 

1.357 

253.3 

0.4628 

136.4 

0.300 

2.333 

392.4 

0.3551 

163.2 

0.250 

3.000 

639.9 

0.2547 

213.3 

0.200 

4.000 

1132 

0.1661 

233.2 

0.150 

5.667 

2292 

0.09355 

404.4 

0.100 

9.000 

5933 

0.04075 

659.3 

(vii)  Viscosity 

r  Ratio  = 

1000 

1.000 

0 

0 

1.000 

— 

0.9375 

0.01265 

0.01299 

1.595 

1.026 

0.9750 

0.02564 

0.04349 

1.537 

1.696 

0.9625 

0.03396 

0.1086 

1.574 

9  700 

z  •  /  o  o 

0.9500 

0.05263 

0.2259 

1.553 

4.292 

0.9375 

0 . 06666 

0.4131 

1.543 

6.197 

0.9250 

0.03103 

0.6382 

1.527 

8.433 

0.9125 

0.09539 

1.069 

1.511 

11. 15 

0.9000 

0.1111 

1.  531 

1 . 494 

14.23 

0.3075 

0.1263 

2.231 

1.477 

17.59 

(viii) 
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Viscosity  Ratio 


Maximum  Pressure 
Gradient  Reduction 
Factor 


1 

1 

3 

1. 133 

10 

1.329 

30 

1.449 

100 

1.523 

300 

1.569 

1000 

1.595 

Table  A2  (d) 

Aspect  Ratio  = 

1:3 

Pressure 

Relative 

Gradient 

Interface 

Insitu 

Input 

Reduction 

Holdup 

Position 

Ratio 

Ratio 

Factor 

Ratio 

(i)  Viscosity 

Ratio  = 

1 

1.00 

0 

0 

1.000 

— 

0.900 

0.1111 

0.04736 

0.9549 

0.4263 

0.300 

0.2500 

0.1631 

0.3599 

0.6526 

0.700 

0.4286 

0.3413 

0.7456 

0.7964 

0.600 

0.6667 

0.6026 

0.6240 

0.9038 

0.500 

1.000 

1.000 

0.5001 

1.000 

0.400 

1.500 

1.660 

0.3761 

1.106 

0.300 

2.333 

2.930 

0.2545 

1.256 

0.200 

4.000 

6.130 

0.1403 

1.532 

0.100 

9.000 

21. 11 

0.04523 

2.346 

(ii)  Viscosity 

Ratio  =  3 

1.000 

0 

0 

1.00 

— 

0.900 

0.1111 

0.08561 

1.006 

0.7705 

0.800 

0.2500 

0.3240 

0.9181 

1.296 

0.600 

0.6667 

1.372 

0.6372 

2.058 

0.500 

1.000 

2.313 

0.5641 

2.313 
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Relative 

Interface 

Insitu 

Input 

Pressure 

Gradient 

Reduction 

Holdup 

Position 

Ratio 

Ratio 

Factor 

Ratio 

0.400 

1.500 

3.836 

0.4378 

2.557 

0.200 

4.000 

13.35 

0.1805 

3.337 

0.100 

9.000 

45.98 

0.05946 

5.109 

(iii)  Viscosity  Ratio  = 

10 

1.000 

0 

0 

1.000 

— 

0.9875 

0.01265 

0.004504 

1.060 

0.3553 

0.975 

0.02564 

0.01312 

1.072 

0.5115 

0.9625 

0.03896 

0.02516 

1.074 

0.6457 

0.9500 

0.05263 

0.04127 

1.070 

0.7340 

0.9375 

0.06667 

0.06217 

1.065 

0.9325 

0.925 

0.08108 

0.1092 

1.058 

1.092 

0.9125 

0.09589 

0.1211 

1.051 

1.262 

0.900 

0.1111 

0.1602 

1.042 

1 . 442 

0.800 

0.2500 

0.7581 

0.9560 

3.033 

0.700 

0.4286 

1.930 

0.3482 

4.503 

0.600 

0 . 6667 

3.769 

0.7294 

5.653 

0.500 

1.000 

6.530 

0.6050 

6.530 

0.400 

1.500 

10.83 

0.4772 

7.250 

0.300 

2.333 

18.67 

0.3443 

7.999 

0.200 

4.000 

37.32 

0.2056 

9.  331 

0.100 

9.000 

131.3 

0.0675 

14.59 

(iv)  Viscosity 

Ratio  = 

30 

1.000 

0 

0 

1.000 

— 

0.800 

0.2500 

1.934 

0 . 9699 

7.737 

0.700 

0.4236 

5.236 

0.3630 

12.22 

0.600 

0.6667 

10.49 

0.7445 

15.74 

0.500 

1.000 

13.39 

0.6201 

13.  39 

0.400 

1.500 

30 . 74 

0.4914 

20.49 

0.300 

2.333 

- 

0.3577 

- 

0.200 

4.000 

- 

0.2147 

— 

Ml  '  1  ;■/; 


V’K.T  d>?, 
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Pressure 


Relative 

Gradient 

Interface 

Insitu 

Input 

Reduction 

Holdup 

Position 

Ratio 

Ratio 

Factor 

Ratio 

(v)  Viscosity 

Ratio  = 

100 

1.000 

0 

0 

1.000 

— 

0.9375 

0.01265 

0.003913 

1.114 

0.7041 

0.9750 

0.02564 

0.03052 

1.  Ill 

1.190 

0.9625 

0.03896 

0.07492 

1. 105 

1.923 

0.9500 

0.05263 

0.1518 

1.097 

2.883 

0.9375 

0.06667 

0.2699 

1.039 

4.049 

0.925 

0.08108 

0.4374 

1.081 

5.395 

0.9125 

0.09589 

0.6614 

1.072 

5.898 

0.900 

0.111 

0.9483 

1.063 

8.  535 

0.800 

0.2500 

6.012 

0.9770 

24.05 

0.600 

0.6667 

33.97 

0.7502 

50.95 

0.500 

1.000 

59.81 

0.6259 

59.81 

0.400 

1.500 

100.2 

0.4971 

66.79 

0.300 

2.333 

171.9 

0.3627 

73.63 

0.200 

4.000 

343.4 

0.2183 

85.35 

0.100 

9.000 

1225 

0.07140 

136.2 

(vi)  Viscosity 

Ratio  = 

1000 

1.000 

0 

0 

1.000 

— 

0.900 

0.1111 

8.760 

1.065 

78.84 

0.800 

0.2500 

58.54 

0.9776 

234.2 

0.700 

0.4236 

164.9 

0.3709 

334.9 

0.600 

0.6667 

335.7 

0.7525 

503.6 

0.500 

1.000 

592.3 

0.6230 

592.3 

0.400 

1.500 

992.9 

0.4994 

661.9 

0.300 

2.333 

1704 

0.3645 

730.3 

0.200 

4.000 

340  3 

0.2197 

850.9 

0.100 

9.000 

12170 

0.7181 

1352 

(vii) 

Viscosity  Ratio 

Maximum  Pressure 
Gradient  Reduction 

Factor 

1.0 

1.0 

3.0 

(1.025) 

10.0 

1.072 

100.0 

1. 113 

1000 

(1.123) 
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TABLE  A3 

Theoretical  data  for  a  viscosity  ratio  of  5.326  and 
an  aspect  ratio  of  7.95:1. 


(a)  Holdup  ratios  and  pressure  gradient  reduction  factors 


Insitu 

Ratio 

Input 

Ratio 

Volume 

Percent 

Oil  in 

Total  Flow 

Holdup 

Ratio 

Pressure- 

Gradient 

Reduction 

Factor 

0.1111 

0.06688 

93.8 

0.6019 

1.596 

0.2500 

0.2025 

83.2 

0.8101 

1.567 

0.4286 

0.4391 

69.6 

1.024 

1.336 

0.6666 

0.8779 

53.4 

1.317 

1.045 

1.0000 

1.758 

36.4 

1.758 

0.7540 

1.  500 

3.709 

21.4 

2.473 

0 . 4940 

2.333 

3.778 

10.32 

3.762 

0.2815 

4.000 

25.81 

3.77 

6.452 

0.1258 

9.000 

133.3 

0.751 

14.37 

0.03137 

(b)  Interface 

velocities 

Volume 

Percent 

Insitu 

Input 

wi 

Oil  in 

Ratio 

Ratio 

Vrc 

Total  Flow 

0.1111 

0.06683 

2.19 

93.8 

0.251 

0.203 

3.00 

33.0 

0.432 

0.445 

3.20 

69.2 

0.675 

0.890 

3.06 

52.9 

1.014 

1.30 

2.74 

35.7 

1.532 

3.90 

2.31 

20.4 

2.41 

9.40 

1.73 

9.62 

4.  IS 

27.0 

1.208 

3.57 

APPENDIX  B 


PHYSICAL  PROPERTIES  OF  THE  OIL  AND  WATER 

The  density  and  viscosity  of  the  oil  were  determined 
as  a  function  of  temperature  by  two  independent  laboratories  - 
the  Instruments  Laboratory,  Department  of  Chemical  and 
Petroleum  Engineering,  University  of  Alberta  and  the  Gasoline 
and  Oil  Laboratory,  Research  Council  of  Alberta.  Measure¬ 
ments  of  both  quantities  agreed  to  within  0.5  percent  at  any 
temperature.  The  density  had  a  low  temperature  coefficient 
and  was  essentially  constant  at  0.319  grn/cc  during  the  experi¬ 
ments.  The  viscosity  is  shown  as  a  function  of  temperature 
over  a  range  of  72  to  32 °F  in  Figure  43.  At  76.5°F,  the 
average  temperature  of  the  oil  during  the  two-phase  experiment 
the  viscosity  was  4.82  centipoise. 

The  density  and  viscosity  of  water  were  obtained  from 

(43) 

the  Handbook  of  Chemistry  and  Physics  .  The  density  of 
the  water  was  approximately  constant  at  0.997  gm/cc  and  at 
7G°F,  the  average  temperature  of  the  water  during  the  two- 
phase  experiments  the  viscosity  v/as  0.905  centipoise. 

The  viscosity  ratio  is  shown  as  a  function  of  tempera¬ 
ture  in  Figure  43,  for  which  both  liquids  are  assumed  to  be 


OIL  VISCOSITY  (CENT  I PO I  SE  ) 
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FIGURE  43  -  OIL  VISCOSITY  AND  VISCOSITY  RATIO  AS  A  FUNCTION 

OF  TEMPERATURE 
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at  the  same  temperature.  Using  the  average  viscosities  given 
above  the  viscosity  ratio  was  5.326. 

The  oil-water  interfacial  tension  was  determined  by 

(3) 

the  capillary  rise  method  developed  by  Bartell  and  Miller'  # 
the  application  of  which  in  the  determination  of  oil-water 
interfacial  tensions  has  been  described  in  a  previous  investi¬ 
gation^0^.  Using  pure  water  and  pure  benzene  the  diameter 
of  the  capillary  was  found  to  be  0.861  mm.  The  subsequent 
determination  of  the  oil-water  interfacial  tension  at  76.5°F 
yielded  42.2  dynes/cm  using  distilled  water  and  41.9  dynes/cm 
using  water  taken  from  the  separator  tank.  The  interfacial 
tension  was,  therefore,  taken  to  be  42  dynes/cm. 


APPENDIX  C 


CALIBRATION  OF  ROTAMETERS 

The  rotameter  type  flow  meters  were  calibrated  by 
weighing  the  amount  of  fluid  discharged  in  a  known  time  for 
different  rotameter  settings.  For  this  purpose  a  special 
temporary  copper  line  was  installed  which  carried  the  flow 
directly  from  the  rotameters  to  the  separator  without  the 
need  for  it  to  pass  through  the  conduit.  The  outlet  of  a 
side  arm  to  this  copper  line  was  positioned  over  an  open  barrel 
on  a  weigh  scale.  Quick  action  gate  valves  enabled  the  flow 
to  be  diverted  instantaneously  from  the  separator  to  the 
barrel  and  vice  versa. 

The  weigh  scale  itself  was  calibrated  with  known  weights 
and  was  found  to  be  accurate  to  within  0.2  lbs  up  to  250  lbs. 

In  general,  a  sufficient  time  was  allowed  to  elapse  so  that 
at  least  100  lbs.  of  fluid  were  collected.  Two  weight  measure¬ 
ments  were  made  for  each  rotameter  setting  and  were  found  to 
agree  very  closely,  usually  to  within  about  0.2  percent. 

Flow  rate,  in  U.S.  gallons  per  minute,  was  plotted 
against  scale  reading  for  each  rotameter  and  the  resulting 
curves  used  to  determine  the  rotameter  setting  necessary  to 
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provide  a  given  flow  rate.  Flow  rates  obtained  in  this  way 
were  accurate  to  within  1  percent.  In  order  to  present  the 
calibrations  as  concisely  as  possible,  flow  rate  in  U.S. 
gallons  per  minute  divided  by  scale  reading,  which  should  be 
a  constant  for  a  well  designed  rotameter,  is  plotted  against 
scale  reading  in  Figure  44  for  the  three  water  rotameters, 
and  in  Figure  45  for  the  three  oil  rotameters.  The  manufacturer, 
tube  number  and  float  number  are  given  for  each  rotameter. 

Also  included  in  the  figure  for  the  water  rotameters  are  the 
manufacturer ' s  curves . 


FLOW  RATE  (  U  S  G  P  M  )/SC  AL  E  READING 
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FIGURE  44  -  CALIBRATION  CURVES  FOR  WATER  ROTAMETERS 


FLOW  RATE  (U  S  G  P  M  )  /  SCALE  READING 
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FIGURE  45  -  CALIBRATION  CURVES  FOR  OIL  ROTAMETERS 
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APPENDIX  D 

THEORETICAL  ANALYSIS  OF  THE  MOTION  OF  A  SPHERICAL  PARTICLE 
RISING  THROUGH  A  PARABOLIC  VELOCITY  DISTRIBUTION 

Consider  a  small  particle,  such  as  an  air  bubble, 
injected  through  the  bottom  of  a  conduit  containing  a  laminar, 
parabolic  velocity  field.  As  the  particle  rises,  it  will  be 
carried  horizontally  by  the  moving  fluid.  For  the  purposes 
of  the  present  investigation  it  will  be  assumed  that  the 
vertical  velocity  Vn  of  the  particle  is  constant.  Actually, 
of  course,  if  the  particle  were  injected  with  essentially  zero 
vertical  velocity  there  would  be  an  initial  acceleration  period 
until  the  particle  reached  its  terminal  velocity.  For  certain 
combinations  of  particle  and  fluid  this  acceleration  period 
would  be  extremely  short  and  in  any  event  could  be  completely 
eliminated  by  having  an  injection  system  which  would  allow 
the  particle  to  attain  its  terminal  velocity  before  entering 
the  conduit  itself.  The  assumption  of  constant  vertical 
velocity  is,  therefore,  entirely  reasonable. 

If  at  time  t  the  co-ordinates  of  the  position  of 
particle  are  y  and  z  with  the  y-axis  vertical,  the  z-direction 
coincident  with  the  direction  of  flow  of  the  fluid,  and  the 
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particle  moving  in  the  plane  x  =  0,  and  at  time  t  +  At  the 

co-ordinates  are  y  +  Ay  and  z  +  Az  then  the  horizontal  velocity 

of  the  particle  at  (y  +  -^)  is  given  approximately  by  — , 

2  At 

provided  At,  Ay  and  Az  are  small  quantities.  Using  strobos¬ 
copic  or  intermittent  illumination,  Az  could  be  determined 
with  the  aid  of  photography  at  various  points  on  the  particle 
trajectory  and  the  velocity  of  the  particle  established. 

If  the  particle  assumed  the  fluid  velocity  instantaneously 
for  all  values  of  y  then  the  particle  velocity  will  be 
representative  of  the  fluid  velocity.  Obviously  the  particle 
has  inertia  and  this  will  not  be  so.  The  object  of  the 
following  theory  is  to  estimate,  by  making  suitable  assumptions, 
the  magnitude  of  the  difference  of  the  particle  and  fluid 
velocities  as  a  function  of  y. 

It  will  be  assumed  that 

(1)  the  particle  has  a  spherical  shape, 

(2)  the  virtual  mass  of  the  particle  should  be 
used  in  Newton's  law  of  motion  instead  of  the  real  mass  of 
the  particle, 


(3)  the  drag  force  exerted  in  the  horizontal 
direction  on  the  particle  by  the  fluid  is  given  by  3t r£/d  (w^  -  w  ) 
where  jj_  is  the  fluid  viscosity,  d  is  the  particle  diameter, 
wx  is  the  fluid  velocity  and  w^  is  the  particle  velocity  and 
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(4)  any  effects  due  to  the  rotation  of  the 
particle  under  the  velocity  gradient  are  negligible. 

The  inertia  of  a  submerged  particle  is  increased  by 
the  fluid  around  it  and  hence 

m*  =  m  +  m'  (30) 

where  m*  =  virtual  mass 
m  =  real  mass 
m'  =  induced  mass 

It  may  be  shown (7)  that  in  the  movement  of  a  spherical 
particle  through  an  incompressible,  nonviscous  fluid  the 
kinetic  energy  of  the  fluid  is  equal  to  that  which  would  be 
possessed  by  a  particle  moving  with  the  saine  velocity  as  the 
particle  and  whose  mass  m'  is  equal  to  half  the  mass  of  the 
fluid  displaced.  It  will  be  assumed  that  the  virtual  mass  of 
a  submerged  particle  can  be  evaluated  in  the  same  way  for  a 
viscous  fluid. 

The  value  of  the  drag  force  is  that  given  by  Stokes ' s 
law  which  is  strictly  true  only  for  solid  particles  moving 
under  steady- state  conditions  at  particle  Reynolds  numbers  less 
than  0.5  in  an  infinite  expanse  of  fluid.  For  a  particle 
consisting  of  a  second  immiscible  fluid  circulation  within 
the  particle  has  to  be  taken  into  account  and  the  corresponding 
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theoretical  drag  force  is  given  by  3ir/Jd.  (w^  -  w  ) 
where  jjj  is  the  viscosity  of  the  fluid  composing  the  particle, 
If  jj_  is  very  much  greater  than  jj_'  ,  as  in  the  case  of  an  air 


IUL+-ML 


3/Ih  3 fJC 


bubble  moving  relative  to  a  liquid,  the  theoretical  drag 
force  becomes  2irj_ld(w^  -  w  )  .  However  Birkhoff reports 
that  the  experimental  evidence  suggests  that  even  for  a  fluid- 
particle  the  drag  is  given  by  3ir  fJi  &  (w^  -  Wp)  ,  i.e.  a  fluid 
particle  behaves  much  like  a  solid  particle,  and  this  value 
will  be  used  in  the  present  analysis. 

Consider  now  a  particle  of  virtual  mass  m*  moving 
with  horizontal  velocity  w  in  a  fluid  having  a  velocity  w^, 
then  application  of  Newton's  law  of  motion  to  the  z-direction 


gives 


m*  w  =  3ir/J_d(v7^  -  w^) 


(Bl) 


where  the  dot  over  Wp  represents  a  derivative  with  respect 
to  time. 

If  the  origin  of  co-ordinates  is  coincident  with  the 
point  where  the  particle  enters  the  conduit,  then  the  parabolic 
velocity  distribution  can  be  described  as 


w 


m 


?<2 

b 


(32) 
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where  2b  is  the  distance  between  the  top  and  bottom  of  the 
conduit  and  w  is  the  maximum  fluid  velocity  at  y  =  b. 


Since 


y  =  vpt 


(33) 


where  Vp  is  the  vertical  velocity  of  the  particle  and  t 
is  time,  equation  (31)  becomes 


w  +  A  w  =  Bt  -  Ct 
P  P 


(34) 


in  which 


A 


3ir  fjjd  =  13  li 


m* 


d2(0  +Ql  ) 

P  2 


(35) 


B  __  6-rr/jd  WmVp  „  36  Wm  Vp  [X 

m*  b  b  d 2  {D  +  Ql  ) 

'P  2 


(36) 


C  = 


3  it  fX  d  WmVp  _  13  Wm  Vp  p 


m*  b^ 


b2a2  (Pp  +  &  ) 


(37) 


where  and.  f)  are  the  particle  and  fluid  densities 
respectively. 

The  Laplace  transform  of  equation  (84)  is 


s  w  (s)  -  wp(o)  +  A  wp(s) 


B 


2C 
^  3 


(33) 


If  the  particle  has  zero  horizontal  velocity  at  t  =  0 


then  Wp (o)  =  0  and 

wp  (s) 


B 


2C 


s 2  (s  +  A)  s3  (s  +  A) 


(89) 
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l .  e . 


wp(s)  = 


B  +  2C] 

1 

2C 

+ 

■ 

B 

+  2C1 

1 

’  B 

.A2  A3  . 

s  +  A 

As3 

.  A 

A2  . 

s2 

La2 

+ 


2 

A 


(90) 


from  which ,  by  taking  the  inverse  transform 


,  JL 

B  +  IS' 

-At  C  , 2  ,  1 

e  -  —  t  +  — 

B  +  IS' 

B  +  — 

A2 

A  . 

A  A 

A  . 

A2 

A  . 

(91) 


a  solution  which  satisfies  the  original  differential  equation 
and  boundary  condition. 

Alternatively,  using  equation  (83) 


w 

P 


a  y 

r 

1 

b  +  2SL 

e  Vp  _  c  ...r  y2  +  A- 

B  +  ~ 

y  -  —2 

B  + 

A2 

A 

AVp2  AVp 

L  A 

Az 

l  aJ 

(92) 


Lee 


E  =  w.F  -  w 


'f  p 

substitute  for  w^  and  w  ,  and  equate  the  derivative  of  E  with 
respect  to  y  to  zero,  then 

A 


(93) 


e  V 


Ym 


P 


1  - 


AB 


AB  +  2C 


(94) 


i.e.  the  value  of  y  for  which  the  difference  between 

the  fluid  and  particle  velocities  is  a  maximum  is  given  by 


J  in  (. 


AB  +  20 


) 


2C 


(95) 


win 
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FIGURE  4-6  -  DIFFERENCE  BETWEEN  FLUID  AND  PARTICLE 

VELOCITIES  AS  A  FUNCTION  OF  y 

The  general  way  in  which  the  difference  between  the 
fluid  and  particle  velocities  varies  with  y  is  shown  in  Figure 
46.  E  is  initially  positive,  passes  through  a  maximum  and 
then  becomes  negative  as  the  particle  approaches  the  top 
surface  of  the  conduit. 

The  particle  velocity  was  evaluated  as  a  function  of 
y  from  equation  (92)  for  two  particular  cases  of  interest. 

(1)  the  movement  of  an  air  bubble  of  diameter  0.5  mm 
in  water  flowing  with  w  =  lft/sec  and  2b  =  1  inch  and 

(2)  the  movement  of  an  air  bubble  of  diameter  0.5  mm 
in  oil  of  viscosity  5  centipoise  flowing  with  wm  =  1  ft/sec 


and  2b  =  1  inch. 
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For  case  (1)  it  was  assumed  that  Vp  =  0.2  ft/sec, 

P: £  =  1  g/cc,  Pp  =  0.00121  g/cc  and  jj_  =  1  centipoise.  Thus 
A  =  140.97,  B  =  1353.3  and  C  =  3247.9. 

For  case  (2)  it  was  assumed  that  Vp  =  0.07  ft/sec, 

Pf  =  0.8  g/cc,  p  =  0.00121  g/cc  and  (1  =  5  centipoise. 

L  p 

Thus  A  =  398.51,  B  =  3018.9  and  C  =  2535.9. 

In  both  cases  the  vertical  velocities  of  the  particles 
were  estimated  using  the  standard  procedure  involving  drag 
coefficient  and  particle  Reynolds  number.  For  case  (1)  the 
particle  Reynolds  number  for  the  vertical  motion  was  approxi¬ 
mately  30  and  for  case  (2)  it  was  approximately  1.8. 

A  comparison  of  the  fluid  and  bubble  velocities  is 
presented  in  Figure  47.  It  is  apparent  that  the  relatively 
low  vertical  velocity  of  the  bubble  rising  in  oil,  combined 
with  the  relatively  high  drag  exerted  by  the  oil  render  the 
inertial  effect  of  the  bubble  and  its  associated  fluid  completely 
negligible.  By  comparison  the  inertial  effect  for  the  bubble 
rising  in  water  is  by  no  means  negligible  and  the  result  is 
a  distorted  velocity  profile. 

For  the  air-water  system  the  value  of  y  for  which  E  is  a 
maximum  as  calculated  from  equation  (95)  is  0.116b,  while  for 
the  air-oil  system  it  is  only  0.0117b. 
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FIGURE  47  -  COMPARISON  OF  FLUID  AND  BUBBLE  VELOCITIES 
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It  would  appear,  therefore,  that  the  proposed  technique 
would  be  quite  adequate  for  relatively  viscous  liquids  with 
viscosities  greater  than  about  five  centipoise  but  should  be 
used  with  caution  when  water  is  the  flowing  fluid. 
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APPENDIX  E 

APPLICATION  OF  STEREOPHOTOGRAMMETRY  TO  THE  MEASUREMENT  OF 
WAVELENGTH  AND  AMPLITUDE 

Stereophotogrammetry  was  used  to  estimate  the  wavelength 
and  amplitude  of  interfacial  waves.  When  the  stereo  pairs  of 
photographs  obtained  with  the  equipment  described  under 
experimental  procedure  were  viewed  under  the  stereoscope  the 
three-dimensional  character  of  the  interfacial  structure  was 
readily  apparent.  In  this  appendix  the  theory  and  calibration 
procedure  used  in  the  determination  of  wavelength  and  ampli¬ 
tude  are  discussed. 

1.  Determination  of  wavelength 

The  determination  of  wavelength  was  relatively  simple 
and  entailed  the  measurement  of  the  wavelength  as  it  appeared 
on  one  of  the  negatives.  This  wavelength  was  then  converted 
into  an  actual  wavelength  by  the  use  of  a  previously  determined 
scale  factor.  Measurements  were  made  on  the  left  hand,  i.e. 
the  downstream,  negative  with  a  transparent  scale  having  0.2 


mm  divisions. 
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With  the  cameras  in  a  fixed  position  the  scale  factor 
was  a  function  of  water  depth  only,  i.e.  the  scale  factor  was 
not  a  function  of  position  on  the  negative.  The  functional 
relationship  was  determined  by  the  use  of  a  flat  plate  marked 
with  half-inch  squares  immersed  at  various  known  depths  in 
water  inside  the  conduit.  As  may  be  seen  from  Figure  48 
the  scale  factor  was  essentially  a  linear  function  of  water 
depth.  When  used  to  calculate  the  actual  wavelength  the  scale 
factor  was  evaluated  at  the  mean  water  depth,  obtained  by 
taking  the  difference  of  the  conduit  depth  and  the  mean  oil 
depth . 

2.  Determination  of  wave  amplitude 

(i)  Theory 

The  determination  of  amplitude  required  the  use  of  a 
considerably  more  complex  procedure  than  that  necessary  for 
the  determination  of  wavelength.  Consider  the  idealized 
arrangement  depicted  in  Figure  49.  Two  cameras,  with  lenses 
L  and  L'  have  photoplanes  AB  and  A'B'  perfectly  parallel  to 
the  datum  plane.  Rays  of  light  from  the  objects  1  and  2  pass 
through  the  lens  centres  without  refraction.  The  heights 
of  the  objects  above  the  datum  plane  are  h^  and  h^  and  the 
camera  lenses  are  both  distance  I-I  from  the  datum  plane. 


SCALE  FACTOR  ( I N  C  H  E  S /l  N  C  H  ) 
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FIGURE  48  -  HORIZONTAL  SCALE  FACTOR 
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2 


FIGURE  49  -  IDEALIZED  STEREO-MODEL 
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Images  of  object  1  occur  at  1^  and  1^',  while  images  of  object 
2  occur  at  I2  and  I2 ' . 

Also  shown  in  the  figure  are  plan  views  of  the  negatives 
showing  the  positions  of  the  images.  The  lines  joining  the 
fiducial  marks  define  an  x-z  co-ordinate  system  with  origin 
at  0  on  the  left  hand  negative.  The  points  0  and  O'  are  also 
the  points  at  which  the  optical  axes  of  the  cameras  intersect 
the  photo  planes. 

The  parallax  of  object  1  is  given  by 

P  =  Z  -  Z'  (96) 

—  _L  — 

where  Z^  and  Z^'  are  distances  from  the  centres  of  projection 
to  the  images  of  object  1.  Z-,  '  is  assumed  positive  and  Z^ 

negative.  Similarly  the  parallax  of  object  2  is  given  by 

p2  =  z2  "  Z2  '  (97^ 

Thomson ("^  has  shown  that  for  the  situation  depicted 
in  the  figure 


Ah12  =  “ 

^  ap12 

1'  P 
f2 

(93) 

where 

Ah12  =  hl 

h2 

(99) 

and 

AP12  =  P1  - 

P2 

(100) 
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Hence,  by  knowing  the  parallaxes  of  objects  1  and  2 


and  the  height  of  one  of  them,  it  is  possible  to  calculate 


the  height  of  the  other.  If  object  1  is  a  point  in  the 


datum  plane  then  h^  =  0  and  equation  (93)  becomes 


h 


AP12 


(101) 


2 


or 


(102) 


Now,  equation  (102)  has  been  derived  assuming  that 


the  photoplanes  of  the  cameras  are  both  perfectly  parallel  to 
the  datum,  that  no  distortion  is  caused  by  the  camera  lenses, 
and  that  no  medium  other  than  air  lies  between  the  objects  and 
the  camera  lenses.  In  addition,  in  order  to  apply  the  equation 
to  parallaxes  measured  under  the  stereoscope  the  negatives 
must  be  positioned  in  exactly  the  same  relative  positions  with 
respect  to  each  other  as  they  were  in  the  cameras.  Slight 
rotation  of  one  of  the  negatives  with  respect  to  the  other, 
while  not  destroying  the  three-dimensional  effect,  causes 
erroneous  parallax  measurements  to  be  obtained.  While  equation 
(102)  could  no  doubt  be  modified  to  take  into  account  the 
effects  of  tilt  of  the  photoplanes,  the  presence  of  dense  media 
and  relative  orientation  of  the  negatives,  this  would  necessarily 
be  a  very  complex  procedure.  An  empirical  approach  was. 
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therefore,  used  to  develop  adjustments  to  be  applied  to  the 
crude  elevations  calculated  from  equation  (102) . 

The  possible  effects  resulting  from  the  departures  from 
the  idealized  situation  for  which  equation  (102)  applies  will 
now  be  considered.  In  doing  so  reference  will  be  made  to 
Figure  50  where  the  object,  the  oil-water  interface,  is  hori¬ 
zontal  and  parallel  to  the  datum  plane  which  is  the  outside 
surface  of  the  bottom  of  the  conduit.  Elevations  of  points 
(2)  in  the  interface  with  respect  to  a  point  (1)  in  the  datum 
plane  are  assumed  calculated  from  equation  (102) .  The  effect 
of  relative  rotation  of  the  negatives  is  a  linear  deformation 
of  the  parallaxes  and  hence  of  the  computed  elevations.  The 
interface,  therefore,  appears  tilted.  The  effect  of  tilt  of 
the  photoplanes  is  to  add  to  the  tilt  of  the  interface  and, 
in  addition,  to  introduce  a  second  degree  distortion  so  that 
the  interface  acquires  curvature.  The  exact  effect  of  lens 
distortion  is  not  known  and  it  can  either  be  assumed  to  be 
negligible  or  included  in  the  other  effects.  The  presence 
of  a  water  layer,  with  a  refractive  index  of  1.33  and  a 
1-inch  thick  plastic  layer  with  a  refractive  index  of  1.48, 
between  the  interface  and  the  cameras  is  primarily  to  reduce 
the  elevations  considerably  and  make  the  interface  appear 
nearer  to  the  datum  level  than  it  actually  is.  In  addition, 
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FIGURE  50  -  ELEVATION  ADJUSTMENTS  FOR  ACTUAL  SYSTEM 
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the  refraction  effects  may  also  introduce  a  second  degree 
distortion  of  the  interface.  The  sum  of  the  individual 
effects  may  well,  therefore,  cause  the  points  2  to  appear  to 
be  at  2 '  and  a  section  of  the  interface  AA  to  appear  in  the 
position  A' A' . 

For  fixed  camera  positions  the  adjustments  to  be  applied 
for  the  effects  of  tilt  of  the  photoplanes,  relative  rotation 
of  the  negatives  and  lens  distortion  will  be  functious  only 
of  the  position  of  the  point  on  the  negative.  Adjustments 
for  the  effect  of  refraction  through  the  dense  media  besides 
being  a  function  of  position  on  the  negative  are,  of  course, 
also  functions  of  the  water  depth.  It  is,  therefore,  convenient 
to  define  two  adjustments  -  the  basic  adjustment,  Ah^  which 
is  independent  of  the  water  depth  and  the  refraction  adjust¬ 
ment,  Ahr  which  is  a  function  of  the  water  depth.  The  adjust¬ 
ments  Ah-^  and  Ahr  are  shown  in  the  figure  as  well  as  the 
crude  elevation,  Ahc  calculated  from  equation  (102)  .  Obviously 
the  datum  plane  is  also  distorted,  by  an  amount  -Ah^,  relative 
to  point  1.  The  true  datum  plane  BB  appears,  therefore,  to 
be  at  B'B'.  With  the  assumption  that  the  deformations  of  a 
flat  interface  can  be  represented  by  terms  of  second  degree 
in  x  and  z,  the  basic  adjustment  can  then  be  defined  by 


. 
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(103) 


where  the  coefficients  c^  are  constant  only  for  each  viewing 
of  each  stereo  pair  and  Ah  by 


(104) 


Ahr 


where  the  coefficients  a^  are  functions  of  the  water  depth 
only  provided  the  cameras  remain  in  a  fixed  position,  and 
where  x  and  z  are  the  co-ordinates  shown  on  the  left  hand 
negative  in  Figure  49.  The  z-direction  coincided  with  the 
direction  of  fluid  flow  and  the  x-direction  was,  therefore, 
normal  to  it. 


In  an  investigation  of  the  variation  of  wave  amplitude 


across  the  width  of  the  conduit  the  general  forms  of  equation 
(103)  and  (104)  would  be  necessary.  However,  in  order  to 
reduce  considerably  the  amount  of  photographic  interpretation 
and  computation  required  for  the  present  project,  the 
determination  of  the  wave  amplitude  was  limited  to  the  centre 
line  of  the  conduit,  i.e.  to  the  z-axis  of  the  co-ordinate 
system.  With  x  =  0  equations  (103)  and  (104)  reduce  to 


&hb  =  co  +  C2Z  +  C5Z 


2 


(105) 


Ah 


a  +  az  +  a_z 
o  2  5 


2 


r 


and 


(106) 


. 
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The  actual  elevation  Ah  was,  therefore,  calculated 

ci 

from  the  expression 

Ah  =  Ah  +  Ah,  +  Ah  (107) 

a  c  b  r 


It  is  apparent  that  in  order  to  determine  the  coefficients 
in  equations  (105)  and  (106)  at  least  three  values  of  Ah^  and 
Ahr  must  be  evaluated  at  different  values  of  z.  Since  the 
values  of  Ah^,  Ahr  and  z  are  subject  to  experimental  error 
it  is  better  to  use  more  than  the  minimum  number  and  a  least 
squares  approach  to  determine  the  "best"  values  of  the 
coefficients.  Thus  if  N  sets  of  corresponding  values  of  Ah^ 
and  z  are  available  then  cQ,  C2  anc^  c5  are  given  by  the  solution 
of  the  matrix  equation 


N 

V  2 

izi 

co  ' 

lA^bi 

I*i2 

izx3 

>  < 

C2 

>  =< 

ZziAhbi  > 

,Izi2 

Izi3 

izi4 . 

■  C5< 

■  Izi2Ahbi. 

(108) 


where  i  =  1, 2 , 3  .  .  .  N  N  >  3 

Similarly  if  M  sets  of  corresponding  values  of  Ahr 
and  z  are  available  then  a  ,  a^  and  a^  are  given  by  the 
solution  of 


M 

y z  • 
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/ 
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where  i  =  1,2,3  ...  M  M  ^  3 

(ii)  Calibration  procedure 

The  calibration  procedure  consisted  of  the  determination 
of  the  coefficients  c  . c„  and  cC(  and  a  ,  a0  and  ac.  The 
coefficients  c^  must  be  determined  for  each  viewing  of  each 
stereo  pair  of  negatives  while  the  coefficients  a^  are  a 
function  of  the  water  depth  only. 

In  order  to  determine  the  coefficients  c^,  seven  crosses 
were  scribed  on  the  outside  surface  of  the  bottom  of  the 
photographic  section  which  was  levelled  to  within  0.002  inches. 
These  crosses  were  positioned  such  that  they  appeared  more 
or  less  coincident  with  the  z-axis  of  the  negatives  as  is 
apparent  in  the  positive  stereo  pairs  reporduced  in  Figures  40 
and  41.  The  fiducial  marks  are  also  visible  in  these  figures. 
The  left  hand  cross  was  selected  as  the  reference  point  for 
all  elevation  calculations.  The  elevations  of  the  other  six 
points  were  calculated  from  equation  (102)  and  since  these 
elevations  were  not  subject  to  the  refraction  adjustment  they 
were  a  direct  measure  of  the  basic  adjustment  Ah-^,  only  with 
the  sign  changed. 

In  using  equation  (102)  to  calculate  elevations  of 
points  (2)  relative  to  the  reference  point  (1)  it  was  necessary 
to  know  the  distance  of  the  lens  centres  from  the  datum  plane. 
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PI,  the  parallax  of  the  reference  point  P^,  and  the  differences 
in  parallax  AP-j^*  Since  the  cameras  were  always  in  the  same 
positions  H  and  P^  were  constants  and  equal  to  677  mm.  and 
34.5  mm.  respectively.  The  differences  in  parallax  were 
obtained  by  the  use  of  the  parallax  bar,  the  use  of  which  has 
been  described  by  Baker ^  .  The  stereoscope  was  fitted  with 
eyepieces  having  a  magnification  of  eight  times.  With  the 
negative  plates  firmly  attached  to  the  light  table  readings 
were  recorded  to  1  micron,  although  they  were  reproducible 
only  to  within  10  microns.  The  expected  error  in  the 
elevations  calculated  from  equation  (102)  was ,  therefore, 
approximately  0.2  mm. 

Equation  (103)  was  solved  by  using  Jordon's  method  for 
matrix  inversion  and  the  University  of  xAlberta  IBM  1620 

computer.  The  matrix  inversion  sub- routine  was  available  as 
a  library  program  from  the  computing  centre.  The  program  was 
written  in  Forgo.  As  a  check  on  the  calculations,  the  com¬ 
puter  calculated  the  sum  of  the  errors,  error  being  defined 
as  the  difference  between  a  particular  Ah-j_  and  the  value 
calculated  from  equation  (105)  using  the  computed  values  of 
cQ,  C2  an^  c^.  In  all  cases  the  error  sum  was  essentially 
zero,  as  it  should  have  been.  The  errors  were  also  used  to 
calculate  a  standard  deviation  which  varied  between  0.1  and 
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0.3  ram.  Since  the  cameras  were  rigidly  fixed  in  position  and 
because  parallax  measurements  were  confined  to  the  z-axis 
where  the  effect  of  relative  rotation  of  the  negatives  is  a 
minimum  the  values  of  c  ,  c2  and  c^  were  approximately  the 
same  for  each  stereo  pair. 

The  coefficients  aQ,  a2  and  a^  in  the  equation  (106) 
were  determined  in  the  following  way.  A  glass  plate,  painted 
white  and  ruled  with  lines  half  an  inch  apart  was  submerged 
at  various  known  water  depths  in  the  conduit.  Parallax 
readings  were  obtained  for  the  seven  crosses  on  the  bottom 
of  the  conduit  and  for  fifteen  points  on  the  grid  plate.  For 
each  stereo  pair  the  basic  adjustment  was  calculated  in  the 
manner  outlined  above.  This  adjustment  was  used,  together 
with  the  crude  elevations  calculated  from  equation  (102)  for 
points  on  the  plate,  to  determine  the  refraction  adjustment 
from  equation  (107) .  The  coefficients  aQ,  a2  and  a^  were  then 
obtained  by  again  using  Jordon's  method  and  the  IBM  1620 
computer  to  solve  the  matrix  equation  (109)  . 

The  coefficients  aQ/  a2  and  a^  are  shown  as  functions 
of  the  water  depth  in  Figure  51.  The  coefficient  aQ  was  a 
linear  function  of  water  depth,  the  slope  being  0.295  mm/mm. 
The  coefficients  a2  and  a^  scattered  somewhat,  probably  due 
to  errors  in  the  parallaxes.  However,  the  coefficients  were 
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FIGURE  51  -  REFRACTION  ADJUSTMENT  COEFFICIENTS 
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relatively  small  and  differed  little  from  zero.  Consequently 
both  a2  and  a^  were  assumed  to  be  zero  for  the  subsequent 
calculations  of  wave  amplitude. 

In  order  to  check  the  calibration  procedure  and  its 
use  in  the  determination  of  wave  amplitude  a  plastic  model  of 
a  square  wave  having  a  wavelength  of  one  inch  and  an  average 
trough  to  crest  height  of  2  mm  was  submerged  in  water  in  the 
conduit.  Parallax  bar  readings  were  obtained  for  points  on 
the  crests  and  troughs  of  the  model.  The  actual  and  computed 
elevations  or  water  depths  are  shown  in  Figure  52.  The 
agreement  is  very  reasonable,  considering  the  error  involved 
in  the  parallax  bar  readings.  The  standard  deviation  of  the 
computed  from  the  actual  elevations  is  0.35  mm.  In  computing 
the  refraction  adjustment  the  water  depth  used  to  evaluate  aQ 
was  the  actual  water  depth.  In  general,  however,  the  water 
depth  is  not  known  and  a  trial  and  error  procedure  is  necessary. 
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FIGURE  52  -  COMPARISON  OF  COMPUTED  AND  ACTUAL  ELEVATIONS 
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APPENDIX  F 

EXPERIMENTAL  DATA 


Data  for  the  flow  of  oil  alone  are  given  in  Table  Fl, 
data  for  the  flow  of  water  alone  in  Table  F2  ,  and  data  for  the 
two- phase  system  in  Table  F3.  Table  F3  is  subdivided  according 
to  the  superficial  oil  Reynolds  number,  ReQ,  which  was  calcu¬ 
lated  using  a  viscosity  of  4.82  centipoise,  i.e.  the  viscosity 
corresponding  to  the  average  oil  temperature  of  76.5°F.  All 
other  quantities  involving  the  oil  viscosity  were  calculated 
using  the  actual  oil  viscosity  prevailing  in  each  experiment. 

The  oil  depth  readings  reported  are  based  on  the  average  conduit 
depth  of  1.007  inches  and  were  used  in  the  computation  of  the 
actual  water  velocity.  Since  the  water  viscosity  varied  only 
slightly  during  the  two-phase  experiments  the  average  value 
of  0.905  centipoise,  corresponding  to  a  temperature  of  76°F, 
was  used  in  the  calculation  of  the  water  Reynolds  number. 

A  set  of  sample  calculations  follow  each  table  to 
illustrate  how  each  quantity  was  computed. 
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Sample  calculations  for  Table  Fl 


For  experiment  number  101  oil  flow  rate  -  13.8  USGPM 


Average  velocity  = 


flow  rate 


conduit  cross-sectional  area 


(13.8)  (0.1336)  (144) 
(60)  (8.01)  (1.007) 


0.548  ft/sec 


From  equation  (14)  Reynolds  number  = 


Pv  De 

fl 


(0.819)  (62.4)  (0.548)  (1.786)  (3600) 
(4.95)  (2.42)  (12) 


=  1250 


From  equation  (13)  Fanning  fraction  factor  = 


1  djp  9c 

2  dz  p Vz 


(-0.104)  (32.2)  (1.786) 

(2)  (0.819)  (62.4)  (12)  (0.548)  2 


0.0162 
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Sample  calculations  for  Table  F2 


For  experiment  number  115  water  flow  rate  =  3.99  USGPM 

flow  rate 


Average  velocity  =  - - 

conduit  cross-sectional  area 


=  (3.99)„C0^1.3,36I.,(144)  =  0.1585  ft/sec 

(60)  (8.01)  (1.007) 


From  equation  (14)  Reynolds  number  = 


Pv  De 


(0.997)  (62.4)  (0.1585)  (1.786)  (3600) 

(0.953)  (2.42)  (12) 


=  2300 


From  equation  (13)  Fanning  friction  factor  = 


1  dp 

2  bz 


gc  De 

P  V2 


(-0.00605)  (32.2)  (1.786) 

(2)  (0.997)  (62.4)  (12)  (0.1585)2 


0.00923 
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Sample  calculations  for  Table  F3 

For  experiment  number  215 


oil  flow  rate  =  1.20  USGPM 
Superficial  average  oil  velocity,  VQ  = 

_  _ oil  flow  rate _ 

conduit  cross-sectional  area 


(1.20)  (0.1336)  (144) 

(60)  (3.01)  (1.007) 


=  0.0476  ft/sec 


Superficial  oil  Reynolds  number 


Po^pDe 

=  Po 

=  (0.319)  (62.4)  (0.0476)  (1.736)  (3600) 

(4.32)  (2.42)  (12) 


=  111.5 


Water  flow  rate  =  7.00  USGPM 
Superficial  average  water  velocity,  Vw 

=  (7.00)  (0.1336)  (144) 

(60)  (8.01)  (1.007) 


=  0.278  ft/sec 


. 
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Superficial  water  Reynolds  number 

PwVwDe 

=  (0.997)  (52. 4)  (0.278)  (1.786)  (3600) 

(0.905)  (2.42)  (12) 


=  4240 


Input  ratio  =  —  =  5.34 

1.20 


Input  oil  fraction  =  - — — -  =  0.1465 

1.20  Jr  7.00 


.  water  depth  1.007  -  0.320  ^ 

Insitu  ratio  =  - 1 -  -  -  =  2. 


oil  depth 


0.320 


15 


Equation  (11)  reduces  to 


^  =  1026  jl  Q 


for  the  particular  conduit 
the  pressure  gradient  for 


used  for  the  experiments 
the  flow  of  oil  alone, 


dz 


Hence 


o 


=  (1026)  (4.35)  (2.42)  (1.20)  (0.1336) 

(3600)  (60) 

=  0.00394  lb/ ft2/ ft 
Pressure  gradient  reduction  factor 


QjlQQJIM  =  0.203 
0.0429 


1026  M0Q0 
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From  a  plot  of  the  data  presented  in  Table  A3 (a)  the 
theoretical  pressure  gradient  reduction  factor  for  an  input 
ratio  of  5.84  is  0.365. 

Hence  the  theoretical  two-phase  pressure  gradient 
assuming  laminar  flow  in  both  phases 


=  ,0- ,00 3 94  =  0.0245  lb/ft2/ft 
0.365 

v  =  _i.  tee]  — i — 

rc  4ab  \dz/  1026 fl 

_ (14-4)  (0.0245)  (3600) _ 

(4)  (4.005)  (0.5035)  (1026)  (4.85)  (2.42) 


=  0.130  ft/sec 


Velocity  ratio, 


0.285 

0.130 


2.19 


Actual  average  water  velocity 

=  (7.00)  (0.1336)  (144) 

(60)  (8.001)  (0.687) 


0.407  ft/sec 


